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EXECUTIVE SUMMARY 


The Center for Advanced Space Propulsion 
(CASP) is one of sixteen NASA sponsored Centers es¬ 
tablished to facilitate the Commercial Development of 
Space. It is the only Commercialization Center with 
a technology concentration in propulsion. CASP has 
completed the second year of operation with a signif¬ 
icant list of accomplishments: 

The number of industrial partners has increased 
to nine (9), almost double the original number. 

The number of active projects has increased to 
fourteen (14), almost three times the original 
number. Fifteen (15) other projects are tn vari¬ 
ous stages of development. 

Four projects leading to flight experiments are 
being actively pursued. 

The selection of a contractor for the CASP of¬ 
fice/laboratory building is underway. Occupancy 
is expected in the Fall of 1990. 

Fifteen Technical papers were published or pre¬ 
sented in national conferences by CASP principal 
investigators distributing results of their research 

effort. 


The first CASP Annual Technical Symposium 
was held with over 70 people in attendance. Posi¬ 
tive feedback on the one-day of presentations was 
received and the symposium will probably be ex¬ 
panded to two days next year. 

During the past year we have continued to solid¬ 
ify and add projects in five focus areas: 

• Advanced Chemical Propulsion 

• Electric Propulsion 


• AI/Expert Systems 


• Microgravity Fluid Management 
a Propulsion Materials Processing 


It has become evident that CASP projects fall 
into two general categories; 1) technology develop¬ 
ment usually supported by the larger Aerospace com¬ 
panies and the NASA field centers, and 2) product 
development and flight experiments usually found in 
teaming arrangements with smaller companies. Fl¬ 
ares I and 2 depict the relationships of out current 
projects with their anticipated role in CASP s de 
opment. 
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Figure 2. 


During the first two years CASP has developed 
the capability to interface effectively with NASA and 
industrial partners. These communications have ma¬ 
tured into a concentrated project workload in the mi¬ 
crogravity and chemical propulsion areas. 

In order to meet CASP long term objectives, em¬ 
phasis will be placed on the following: 

• Develop Electric Propulsion 

• Develop and Process Flight Experiments 

• Form Partnerships with Small Business 

To adequately respond to the increasing interest 
by both industry and government in electric propul¬ 
sion, CASP, CAR, and UTSI elected to jointly fund 
the refurbishment of a 10 x 28 ft. space chamber. 
Upon completion, the chamber will include provisions 
for state-of-the-art laser diagnostics to support elec¬ 
tric thruster testing and subsequent space flight ex¬ 
periments. 

during the year Dr. George Garrison assumed 
the position of Director and Dr. F. Speer became 
CASP’s Associate Director. CASP continued to as¬ 
semble a uniquely qualified team of investigators and 
two additional faculty were added as Principal Inves¬ 
tigators. At the close of the period, requisitions for 
three propulsion specialists were being advertised by 
UTSI. 


Personnel census at CASP s interim facility in 
the UTSI Research Park increased this year. Twenty 
seven investigators are now involved in fourteen active 
projects. CASP’s microcomputer center was com¬ 
pleted with the addition of one Compaq and four 
CompuAdd 386 microcomputers and a SUN 386i 
workstation. 

To bring together the growing number of inves¬ 
tigators and assistants working on CASP projects, 
a 10,000 ft . 2 office and laboratory facility is being 
funded by the State of Tennessee and the Appalachian 
Regional*Commission. Site selection, design and floor 
plan layout of the new building have been completed, 
and bids for construction are being reviewed. Ground 
breaking is scheduled for early January, 1990. 

NASA code C organized and hosted quarterly 
meetings of the CCDS directors, NASA field center 
representatives, and Office of Commercial Program 
(OCP) support contractors. There is strong support 
for the CCDS program at NASA Headquarters and, 
in spite of budget reductions across all of NASA, 
grants for the commercialization centers remained un¬ 
changed. 

It has become very evident that the CCDS s are 
expected to fly experiments and that OCP will do 
everything possible to facilitate flights and shorten, 
where possible, the time from concept to flight. CASP 
continues to evaluate the skill and experience required 
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to meet future (light experiment needs and identify 
the specific role the Center will take in the devel¬ 
opment of a follow-on to the current CCDb sound¬ 
ing rocket program. The matrix, show in figure 3, 
displays the tradeoffs and implied management chal¬ 
lenges in developing future flight vehicles. 

CASP ended the current year slightly under the 
financial goals set in the original 5-year plan, but fin¬ 
ished the year with income from grants and cost shar¬ 
ing being slightly larger than expenses Growth in 
project funding was slower than expected and ll ^in 
contributions closely tracked the plan. Future CASP 
growth will demand greater income. The total was 
far below that necessary to support a flight program. 

Along with the continuing challenge of daily op¬ 
erations, CASP’s management emphasis now shifts to 
the next major set of design, fabrication, testing an 

flight of a commercially viable propulsion experiment. 
Of CASP’s focus areas microgravity fluids manage¬ 
ment and electric propulsion related systems are ma¬ 
terializing as prime candidates for space flight exper 
iments. 
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INTRODUCTION & OVERVIEW 


The Center for Advanced Space Propulsion is 
part of The University of Tennessee-Calspan Center 
for Aerospace Research (CAR). Formed in 1985 to 
take advantage of the extensive research faculty and 
staff of The University of Tennessee and Calspan Cor¬ 
poration, CAR is a not-for-profit organization estab¬ 
lished to conduct aeronautical and aerospace research 
and development. The broad technical knowledge and 
assistance provided by CAR make it a strong par¬ 
ent organization for the Center for Advanced Space 
Propulsion. 

On October 31, 1989 CASP marked it’s second 
successful vear as a NASA Center for the Commercial 
Development of Space (CCDS). In carrying out it's 
mission to strengthen U.S. competitiveness m space 
propulsion technology, CASP added new industry 
partners to the Management Advisory Council and 
new projects to it’s list of active pursuits. Based on 
investigators’ qualifications in propulsion system de¬ 
velopment. and matching industry strong intent the 
Center focused its efforts in the following technical 
areas: 


• Advanced Chemical Propulsion 

• Electric Propulsion 

• AI/Expert Systems 

• Fluids Management in Microgravity 

• Propulsion Materials Processing 

Throughout the year, the Center continued to 
integrate academic excellence with commercial space 
requirements and increased the number of active 
projects to 15. These are supported by 10 Principal 
Investigators and 19 Graduate Research Assistants 
and six staff. Each project has a defined commercial 
application and either enables or enhances the related 
technologies with respect to space commercialization. 

CASP marked it’s second anniversary in the Cen¬ 
ter's interim facility at the UTSI Research Park. De¬ 
sign aiul funding of CASP’s new 10,000 ft.' office 
building was completed during the year. Construc¬ 
tion is scheduled to begin in late November, 1989. 

Valuable inputs from both CASP’s Management 
Advisory Council and NASA's Office of Commercial 
Programs formed a framework for CASP’s perfor- 
mance goals. Figures l and '2 of the Executive Sum¬ 
mary depict the relationship of the various projects 


accomplished by CASP during the year, and relates 
each to it’s proposed objective. Critical for reach¬ 
ing these goals are the ground based facilities which 
are used by the investigative teams in performing the 
project work. Supercomputer support continued to 
be provided by MSFC and was the empowering tool 
for our advances in Computational Fluid Dynamics. 
Installation at UTSI of an Alliant computer to compli¬ 
ment the existing VAX allowed some programs which 
had been run on the NASA Marshall CRAY to be 
moved to UTSI. To supplement these computational 
systems, CASP completed outfitting it’s interim com¬ 
puter laboratory with three additional 386 PC’s and 
a Sun 386i. With these existing computers, CASP 
offers a full spectrum of computational hardware for 
it’s researchers. 


Earlier emphasis on marketing paid dividends 
this year with addition of two projects to the ac¬ 
tive list and three more industrial organizations which 
have expressed an interest in joining CASP in future 
commercial space ventures. Figure 3 tabulates the 
projects which have a potential for being funded dur¬ 
ing the next vear. Those most closely matching the 
commercial objectives of CASP will be activated. 


CASP’s second vear was marked by increased fi¬ 
nancial stability, but there is a continuing problem 
with the extended time to consummate project agree¬ 
ments. Project costs continued to vary somewhat 
from their budgets, however total CASP spending 
ended somewhat less than earlier projections. Cash 
flow problems caused by delayed project funding con¬ 
tinued to restrict programs in some areas. To accom¬ 
modate the shortfall, management action was taken 
to contain costs. CASP ended the year with expen- 

/• . i e J _ ,1 nkla 




NASA gave impetus for increasing the CCDS 
presence in space by providing a number of GAS can 
and mid deck experiment opportunities. As briefed 
bv NASA, flight openings are becoming available, on 
Shuttle, Spacehab, and the Space Station. Such re¬ 
served” space allows CASP’s industrial partners to 
take immediate advantage of space experimentation 
without encountering a lengthy queue. In prepara¬ 
tion for possible space flight, CASP selected experi¬ 
ments in Microgravity Fluids Management and Elec¬ 
tric Propulsion for flight in the ’92-'93 time frame. 


A precursor micro-gravity fluids management ex¬ 
periment. planned to be flown on a KC-135 during the 
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summer of 1990, was being formulated with MSFC as 
the year ended. Knowledge gained from this exper¬ 
iment will flow directly into Boeing’s SOFTE space 
flight experiment design. 

As the year ended, CASP’s goals remained un¬ 
changed and designs of several space flight experi¬ 
ments were in progress. The industrial commitments 
and project prospects appear to be solid and should 
result in growth in each of the CASP focus areas. 

In CASP’s annual review with NASA Headquar¬ 
ters, Office of Commercial Programs, Jim Rose ex¬ 
pressed approval of the path chosen and the progress 
made by the Center. He noted the need for elements 
of infrastructure which must be developed to service 
the continuing needs of all users of commercial space. 
CASP has the objective of being a significant player 
in development of one of the critical infrastructure 
elements - transportation. 



ADMINISTRATIVE / BUSINESS STATUS 


Funding / Bus iness Goals 

CASP is one of two operating elements of the 
Center for Aerospace Research, a not for profit orga¬ 
nization founded by The University of Tennessee and 
Arvin/Calspan Corporation. CASP’s primary source 
of funding is an Annual Grant from NASA code C 
with industrial contracts and agreements providing 
additional sources of revenue. 

Figure II.1 describes the major elements of CASP 
funding and reflects the business goals set by CASP 
during the past two years. One significant philo¬ 
sophical change emerged during the year: NASA has 
agreed that a meaningful flight experiment cannot 
generally be formulated, tested, and flown within the 
five year limit imposed by the initial grant. As a re¬ 
sult NASA has accepted the need for Center funding 
beyond the five year limit. CASP’s flight plans mirror 
the requirement for extended funding of core needs. 
It is our belief that a usable microgravity or electric 
propulsion experiment cannot be accomplished before 
FY 93. To reach the point of a flight experiment will 
require a great deal of nurturing (cost sharing) with 
both a longer period of baseline and augmented fund¬ 
ing required to be successful. 

In FY ’89, the Center was not able to meet it’s 
overall goal for funding. While in-kind contributions 
met the planned level, cash from contracts and grants 
was somewhat lower. A carry-over of industrial con¬ 
tributions from the prior year allowed project work to 
proceed at a level which was slightly above the year’s 
income. In general, there is still a shortfall in funds 
caused by delays in the receipt of funds from NASA 
Centers and industrial partners. The effect ot the de¬ 
lays can be seen in figure II.2 where project workload 
was restricted so it remained below the level of funds. 
Note that receipt of the CASP’s NASA Headquar¬ 
ters Grant typically funds all work early in the year 
while the Center finalizes terms with our industrial 
partners. 



Figure II.1 


Personnel 

On February 1, 1989, Dr. George Garrison as¬ 
sumed the position of Director with Dr. F. Speer 
named CASP’s Associate Director. CASP continued 
to assemble a uniquely qualified team of investigators 
and additional faculty were added as Principal Inves¬ 
tigators for the Ion Propulsion and Inconel Welding 
projects. At the close of the period, CASP had filled 
one of the three open requisitions for propulsion spe¬ 
cialists. Mr. Verlin Friedly accepted a position as 
a Research Engineer and will begin working primar¬ 
ily in electric propulsion. He will continue the work 
started by Dr. Y.K. Chen who left late in the year to 
take a position with Ames Research Center. 

Initiation of the Helical Coil project with Sund- 
strand Inc. marked the first opportunity to tap ex¬ 
pertise at the UT Knoxville campus. Dr. Edward 
Keshock, of the Engineering Department and two as¬ 
sistants are investigating two phase flow in helical 
coils for this project. 
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Figure II.2 

Administrative support remained constant dur- 
ing the year and workload peaks were accomplished 
by using temporary hires. 

Facilities 

An expanded project workload and initiation of 
two projects requiring laboratory space to conduct 
ground testing pushed CASP’s interim facility to ca¬ 
pacity. Design of the new CASP Building was com¬ 
pleted and offers the most immediate prospect for re¬ 
lief of the crowded conditions being experienced in 
the interim facility. 

Figure II.3 shows the floor plan of the proposed 
building. A Fiber-optics link will connect the build¬ 
ing with the UTSI Computer Center which is the 
port through which researchers access the CRAY 
at MSFC, and VAX and Alliant at UTSI. Included 
within the 10.000+ ft. 2 building is a small experi¬ 
ment i>uild-up area and laboratory. These areas have 
gained additional importance as CASP researchers 
have begun (light experiment fabrication. 



Figure II.3 

Renovation of the 10’ x 28’ Space Chamber at 
UTSI was initiated during the year. This chamber 
will support electric propulsion system and compo¬ 
nent testing and will be initially configured to perform 
a detailed experimental investigation of ion thrust 
grid erosion. Testing of arc jets and resisto jets can be 
accommodated, but depending on the flow rate, ad¬ 
ditional pumping capacity may be required. The sys¬ 
tem is complete with a LHi pump and liner system, 
roughing pumps, and control system. Test hardware 
includes an X-Y positioning system and laser diag¬ 
nostics equipment. Refurbishment of the chamber is 
scheduled for completion before June, 1990. Funding 
for activation is provided by CAR, UTSI, and CASP 
non-grant resources. 

Addition of a Sun Microsystems workstation and 
three 386 microprocessors completed hardware acqui¬ 
sitions for the interim CASP computer center. Two 
PC’s plus two Visual terminals are connected to the 
UTSI Computer Center switch through Codex 9600 
baud modems. From the switch, investigators link to 
NASA k DOD Supercomputers and VAX and Alliant 
computers within the University system. 


Seminar Pu blications 

In October 1989, CASP hosted it’s first annual 
Technical Seminar. Lasting one day, the seminar al¬ 
lowed each investigator to present the salient aspects 
of his efforts to date and field questions from the m- 
dustrial and Government attendees. Feedback from 
die attendees was favorable and their comments will 
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be included in our planning for the next annual sem- 
inar. 

Publications which have been either presented or 
published by CASP investigators are listed below. 


FY 89 CASP PUBLICATIONS 


AUTHOR 

UTCHFORD, RON 


M. H. McCAY and 7. 0. McCAY, 
eL aL 

M. H. McCAY and T. D. McCAY, 
et at 

a GARRISON 


W.E. DIETZ, Ei_ KJECH and 
MOCNIS AU 


W.E DETZ, HJ. FER8ER and 
MOONIS AU 


LUC GUPTA and MOONtS ALI 


W.E. DIETZ, EL KIECH and 
MOONIS AU 

C. F. LO 
J. ORR 


TITLE 

[SONS-, Master's Thesis. The University of Tennessee. Novemoer ,989 
PULSE LASER WELDING OF NCO 718’, Presented Octooer 1S-20. 1989, ICALEO ,989. 
Mando, Florida 

EFFECT OF PULSE DUTY CYCLE ON INCONEL 7,8 LASER WELDS’. SAE 
-ECHNICAL PAPER SERIES #892299. The Use of Lasers in 
Manufacturing, Septemoer 25-28, 1989 

STATUS OF THE CENTER FOR ADVANCED SPACE PROPULSION’ 
nternatlonal Conference on Engineering Management. Septemoer 1989 

CLASSIFICATION OF DATA PATTERNS USING AN AUTOASSOCIA7TVE N0JRAL 
siETWCRK TOPOLOGY* The Second IntemaUonal Conference on 
ndustnal & Engineering Applications of Artificial Intelligence 
t Expert Systems, Tuilahoma, TN, June 6-9, 1989 

■SPACE SHUTTLE MAIN-ENGINE COMPONENT ASSEMBLY, ASSIGNMENT AND 
SCHEDULING EXPERT SYSTEM*. The Second International Conference 
on Industrial 1 Engineering Applications of Artificial htelSgence 
«, Expert Systems, Tuilahoma, TN, Juna 6-9. 1989 

•LEAOER - AN INTEGRATED AUTOMATIC LEARNING ANO FAULT DIAGNOSTIC 
^PQ,T SYSTEM FOR SPACE SHUTTLE MAN ENGNE (SSME). The Second 
International Conference on Industrial & Engineering Applications 
of Artificial tntelSgence & Expert Systems, Tuflahoma, 

June 6-9, 1989 

-NEURAL NETWORK MODELS APPUED TO REAL-TIME FAULT DIAGNOSIS’ 

Journal of Neural Network Computing. VoL 1, No. 1, 1989. 

INTELLIGENT HYPERTEXT SYSTEMS FOR AEROSPACE ENGNEEFUNG 
TABES '69 Technical and Businesa Exhibition and Symposum. May 16 17. 1989 


RC CODE APPLICATION TO HIGH AREA RATIO NOZZLES* 


, a rnnfamnrii Anfll 20. 1989 


FORMAT 

MASTER'S THESIS 


PRESENTATION 


PRESENTATION 

PUBLICATION 


PRESENTATION 


PRESENTATION 

PROCEEDING 


PRESENTATION 

PROCEEDING 


PRESENTATION 

PROCEEDING 


PUBLICATION 


PRESENTATION 

PRESENTATION 
MASTER'S THESIS 
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STATUS REPORTS FOR ACTIVE PROJECTS 


SECTIONS 

Advanced Chemical Propulsion 
A.I. Propulsion Applications 
Microgravity Fluid Management 
Electric Propulsion 
Laser Material Processing 
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advanced chemical propulsion 


• OMV Variable Thrust Engine Analysis 

• Spray Combustion Stability 

• Advanced High Area Ratio Nozzles 
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Project Title: 

Principal Investigators: 


OMV/VTE Variable Thrust Engine Analysis* 
L. Larosillere, R. Litchford and S. M. Jeng 


Project Report 


SUMMARY 

The objective of the present work is to develop 
a predictive CFD based analytical tool for the vari¬ 
able thrust engine (VTE) in the orbital maneuvering 
vehicle (OMV). This objective is being accomplished 
within the framework of the Los Alamos KIVA com¬ 
puter code for chemically reactive flows with sprays. 
For the OMV application, the main structure of KIVA 
is to be retained while reformulating many of the phe¬ 
nomenological submodels, enhancing some of the nu- 
merics, and adding more features. 

The analytical model consists of the general con¬ 
servation equations for two-phase reactive flows and 
of submodels for turbulence, chemical reactions, and 
bipropellant sprays. Tailoring this model to the OMV 
engine brings about the added complexities of com¬ 
bustion and flow processes that occur in a liquid hy- 
pergolic propellant rocket chamber. 

This report exposes the foundation upon which 
the analytical tool is being constructed and devel¬ 
oped. Results from a cursory computational exercise 
involving the simulation of the flow and combustion 
processes in a hypothetical N 2 H 4 /N 0 O 4 rocket engine 
thrust chamber is presented and discussed. 

INTRODUCTION 

TRW is currently developing a hypergolic bipro¬ 
pellant variable thrust engine to be qualified specif¬ 
ically for OMV. The available analytical prediction 
methodology is incapable of providing physical and 
phenomenological insight about the complex spray 
combastion and flow processes occurring in the en¬ 
gine tftrust chamber. Thus, an adequate CFD model 
would significantly contribute towards the identifica¬ 
tion of existing design flows and the evaluation of al¬ 
ternative new designs. 

Government/Industry Partners - MSFC/TRW 


A complete description of the combustion and 
flow processes occurring in a liquid hypergolic rocket 
chamber does not exist. It’s known that a wide range 
of interrelated and coupled physical phenomena man¬ 
ifest themselves in the combustor. For example, botli 
homogeneous and heterogeneous chemical reactions 
may occur; the spray droplet trajectory and life his¬ 
tory is determined by the chamber gas environment: 
all gas fluid in the chamber is either from the liq¬ 
uid evaporation process or the reaction process in the 
mixed liquids. Many physical processes occurring m 
the VTE are similar to other complex flow phenom¬ 
ena and have been modelled and applied to the gas 
turbine combustor, diesel engine, cryogenic liquid pro¬ 
pellant rocket; however, the advanced CFD model For 
the hypergolic propellant rocket has never been es¬ 
tablished. 

T his report is a current status update of an on¬ 
going effort to develop a predictive tool for the com¬ 
bustion chamber of the VTE. A framework for con¬ 
structing the model is presented in particular, results 
from a cursory computational exercise are presented 
in support of the proposed framework. 

p flYfiT.nPVfENT OF S P Ml COMBUSTION 
^Np FT.OW MODEL 

The case is considered in which liquid sprays of 
fuel and oxidizer are injected into the combustion 
chamber of a rocket engine. Upon entering the cham¬ 
ber volume, the two continuous liquid phases will en¬ 
counter forces leading to the formation of fuel and 
oxidizer drops with a spectrum of properties. Com¬ 
mencing at the point of drop formation, the aim is to 
track the drop dynamics, the gas dynamics, and the 
ensuing interactions. Hence, the realm of two-pluise 
chemically reacting flows is entered. A separate-phase 
flow model is actually used: gas phase is described us¬ 
ing Eulerian coordinates and dispersed-phase liquid 
droplets are described via Lagrangian coordinates. 
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MODEL DESCRIPTION 


Dispersed-Liquid Phase 


Gas-Phase 

The unsteady conservation equations for a tur¬ 
bulent mixture of chemically reactive ideal gases are 
formulated. Specific equations will not be listed here 
because they are well known and are available else¬ 
where. Suffice to say that they are a system of cou¬ 
pled partial differential equations having the following 
structure: 


+ div (pity) ~ "h (1) 

where the dependent variable $ may represent 
bulk velocity, internal energy, gas composition, turbu¬ 
lence kinetic energy, and rate of decay of turbulence 
kinetic energy. The variable T* is an effective diffu¬ 
sion coefficient for <*,, and S<p is an appropriate source 
term representing effects such as droplet/gas interac¬ 
tions, turbulence production, and chemical reactions. 
The primary 

limitations of this formulation are the modeling 
of the source terms and the numerical solution tech¬ 
nique. 

Two models are available to represent the effects 
of turbulence. A standard version of the k-E turbu¬ 
lence model, modified to include volumetric expansion 
effects, and spray/turbulence interactions has given 
results of sufficient accuracy for engineering purposes 
in a variety of applications. The subgrid scale turbu¬ 
lence model 1 is also available. This model reduces to 
the k-E model near walls where all turbulence length 
scales are too small to be resolved by the computa¬ 
tional mesh. 

The confinement of the gas is taken into account 
by the boundary conditions. Velocity boundary con¬ 
ditions on rigid walls can be free slip, no slip, or tur¬ 
bulent law-of-the-wall. For turbulent law-of-the-wall 
conditions, the tangential velocity components are de¬ 
termined by matching to appropriate wall function. 
Temperature boundary conditions on rigid walls are 
introduced by specifying either the wall temperature 
or the wall heat flux. The boundary conditions for the 
turbulence kinetic energy and its decay rate are zero 
gradient of the kinetic energy at the walls and equi¬ 
librium between generation and dissipation of turbu¬ 
lence. Boundary conditions for the individual species 
are zero gradients of the mass fractions at the walls, 
corresponding to chemically-neutral surface and no 
fluxes through the wall. 


Each spray (fuel or oxidizer) is followed by em¬ 
ploying a discrete droplet model which represents the 
spray by a statistically meaningful number of indi¬ 
vidual droplets. A computational droplet represents 
a parcel of droplets all having the same characteristics 
such as position and size. This formulation is taken 
from Sukorvicz 2 who stresses a statistical rather than 
a deterministic approach to the modeling of sprays. 

A droplet’s trajectory is tracked using the bal¬ 
listic equations as derived via Newton’s second law. 
Hence, a Lagrangian technique is implied. The turbu¬ 
lence effects on droplet motion are simulated statis¬ 
tically by superimposing upon the gas field turbulent 
eddies, each having a length, lifetime, and velocity 
fluctuations. This model is based on a Gaussian dis¬ 
tribution of the local gas velocity perturbations. 

Boundary conditions are needed in order to con¬ 
fine the sprays within the chamber walls. A ques¬ 
tion which must be answered is what happens when 
a droplet impinges on a hot rigid wall? Some sketchy 
results obtained by Watchers and Westerling 3 with 
water show that the observed phenomena are strongly 
linked to the Weber number. The technique employed 
in the current model is to set the droplet velocity 
equal to the wall velocity and not to allow any heat 
transfer between the droplet and wall. This is a pro¬ 
visional model for spray/wall interactions and needs 
to be critically reevaluated. Conditions at the loca¬ 
tion of drop formation need to be specified also. This 
is presently achieved by specifying a droplet mass 
flow rate and a distribution of droplet sizes, veloci¬ 
ties, temperatures, and other relevant parameters. 

PHYSICAL SUBMODELS AK D 

MATERIAL P ROPERTIES 

Chemistry 

Both equilibrium and kinetic chemical effects 
need to be accounted for. Unfortunately, chemical 
kinetic mechanisms and rates are not readily a\ail- 
able for the proposed bipropellant combustion. The 
following nine step reaction is being enlisted Tor the 
NiH-i/NiO-i system: 


2N 2 Ih —♦ 2NHz + Hi + Nl ( 1 ) 

N 2 O 4 —- 2N0i (' 2) 

2 NOi - 2NO + O 2 

ff 2 + WO2 - H2O + NO ( 1 ) 
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2 Hi + Oi — 2 HiO ( 5 ) 

2 Hi + Oi + 2 NO —* 2HiO + 2NO (6) 

2NH 3 + 3 NOi — + 3 HiO + 3 NO + Ni (7) 
4 NH 3 + 3 Oi —- 6H 2 0 + 2N 2 (8) 

2 NH 3 + Hi + 20i —” *HiO + Ni (9) 


The rate constants are assumed to be of a gener¬ 
alize Arrhenius form. Reaction rates are determined 
from the law of mass action or any other experimen¬ 
tally derived rule. Once the chemical scheme has 
been established, the chemical source term in the 
species continuity equation and the chemical heat re¬ 
lease term in the energy equation are obtained. 

pmpjgt Evaporation 

The hypergolic propellants of primary interest 
for liquid rocket applications are hydrazine and its 
derivatives - MMH (monomethyl hydrazine) and 
UDMH (unsvmmetrical dimethyl hydrazine). For 
high ambient temperatures, as occurs in rocket cham¬ 
bers, conditions are conducive for formation of an 
exothermic decomposition flame enveloping the fuel 
droplet near the surface as depicted in Fipire 1. This 
process of pure monopropellant combustion can sub¬ 
stantially enhance the gasification rate over the non¬ 
reactive case. If, in addition, an oxidation propellant 
is introduced, as is the case for hypergolic bipropel¬ 
lant systems, the decomposition envelope flame may 
be further surrounded by an oxidation flame where 
the decomposition products react with the oxidizing 
agent. This process of hybrid droplet combustion is 
depicted in Figure 2. The additional presence of the 
oxidation flame would further enhance the gasifica¬ 
tion rate over the pure decomposition case. 

Since an MMH/N 2 0 4 hypergolic bipropellant 
system has been designated for the OMV-VTE, the 
preceding phenomena will occur in the combustor and 
should be properly modeled as part of the analysis 
effort. Depending upon the local flow conditions the 
fuel will gasify in one of three possible modes: (1) pure 
evaporation, (2) pure monopropellant droplet com¬ 
bustion, or (3) hybrid monopropellant droplet com- 
busticih (dual flame model). Since the KIVA code is 
equipped with a pure evaporation model, no modifi¬ 
cations are required in this case; however, theoretical 
models must be obtained and coded for the single and 
dual flame cases. 

Much work has been done in the past on mono¬ 
propellant combustion at steady-state condition, but 
little effort has been devoted to the development of 


transient engineering submodels which include the 
heat-up period and are applicable for time depen¬ 
dent CFD applications. An exception is some work 
by Faeth 4 on pure monopropellant droplet combus¬ 
tion at high ambient temperatures. He considered 
quasi-steady, spherically-symmetric, gas-phase trans¬ 
port coupled to a time dependent constant property 
liquid phase with one step global decomposition reac¬ 
tion rate kinetics of the Arrhenius form. The theory 
also allowed for corrections in transfer rates due to 
forced convection from empirical correlations for the 
film thickness. The gas-phase solution results in a 
system of nonlinear algebraic equations easily solved 
with an iterative algorithm, and the liquid-phase is 
governed by two ordinary differential equations. 

Faeth performed calculations for EN (ethyl ni¬ 
trate) and PGDN (propylene glycol dinitrate) and fit¬ 
ted the results to experimental data by selecting a re¬ 
action order and activation energy and adjusting the 
Arrhenius constant. The procedure correlated exper¬ 
imental data reasonably well and appears justified as 
a good, efficient engineering submodel for pure mono- 
propellant combustion. 

Because Faeth’s theory is good and readily avail¬ 
able it was selected as the basis submodel for imple¬ 
mentation in the KIVA code. Before incorporation 
in the code, however, the model was studied for ap¬ 
plicability to hydrazine type fuels. Calculations were 
made for hydrazine, MMH, and UDMH monopropel¬ 
lant combustion in an atmosphere of inerts anti com¬ 
pared with the experimental results of Allison and 
Faeth 5 . The computations were made such that the 
process reached steady-state and the droplet diameter 
attained the size corresponding to experiment. The 
kinetic parameters were adjusted to fit the data as 
prescribed by Faeth with good correlation obtained. 
The comparison of theory with experiment is depicted 
in Figures 3, 4, and 5 for hydrazine, and MMH, and 
UDMH, respectively. With the activation energies 
chosen for hydrazine and MMH, the burning rate con¬ 
stant falls off sharply at decreasing temperatures. If 
a smaller activation energy is chosen, the theory may 
be extended to lower temperatures. Also note that, 
a reaction order of unity was specified in the cur¬ 
rent calculations while a second order reaction may 
be necessary to correlate data at higher pressuies as 
stated by Faeth. These issues will be examined and 
addressed in future study. Calculations for pure e\up 
oration are also shown as dashed lines on the figures 
to indicate enhancement of the gasification rate. 

A transient model for hybrid droplet combus- 


19 



tion is not available at this time although Allison and 
Faeth 5 have formulated and correlated a steady-state 
hybrid combustion model for hydrazine type fuels. 
At this time, two alternatives are being considered. 
They are to modify Allison and Faeth’s theory to a 
quasi-steady approach or make ad hoc modifications 
to the pure monopropellant combustion model. Fu¬ 
ture study will determine the approached deemed ap¬ 
propriate. 

Because it is unwise to expend effort modify¬ 
ing the KIVA code until complete droplet gasification 
models have been established and experimentally val¬ 
idated, the original KIVA pure evaporation model will 
be retained during modifications to other parts of the 
code structure. Artificial amplification of the gasifi¬ 
cation rate may be used to reflect actual combustion 
processes for code development purposes until model 
development is completed. 

p r np/Dron Interactions 

Dense sprays exist in a region close to the in¬ 
spector. Hence, considerations of drop collisions are 
warranted. There are two types of drop/drop interac¬ 
tions: collisions between drops of the same kind and 
collisions between two different kinds of drops (i.e. 
fuel and oxidizer). A literature search and review is 
being conducted in order to formulate an appropriate 
model. 

Thermochemical Data Base 

Thermophysical properties for hydrazine and its 
derivatives have been collected and placed in a block 
data subprogram. Properties for N0O4 are also avail¬ 
able. 

CALCULATION TECHNIQUE 
Crid Generation 

A rapid grid operation technique has been con¬ 
structed specifically for the VTE geometry. Further 
developments are proceeding as the program evolves. 
The grid consists of arbitrary hexahedrons and is able 
to be arbitrarily packed. 

Numerics 

Presently, the numerics are taken exactly from 
KIVA 1 . KIVA employs a fully 3-D finite volume ap¬ 
proach over a staggered mesh. Time marching is per¬ 
formed by use of operator splitting and the solution 
is advanced in a series of phases. Variable implicit¬ 


ness is employed to circumvent the current stability 
requirements. Enhancements to the inflow/oiitllow 
boundary condition treatment are contemplated. 

EESHUS AND DISCUSSION 

A cursory calculation is presented for an assess¬ 
ment of the proposed modeling framework. It is per¬ 
formed using a 15 x 40 grid (Fig. 6), covering a cham¬ 
ber radius of 1.83 cm and a length of 12 cm. A small 
amount of mesh packing is done near the walls. At 
time t=0, pressure, temperature, and density are set 
uniform in the chamber. The gas composition at t=0 
is virtually CO2 with initial temperature of 1600 °k 
and pressure of 4.5 atm. A constant back pressure 
equal to the initial pressure and adiabatic wall condi¬ 
tions are applied for all times. 

Fuel (N 2 H 4 ) and oxidizer (N 2 O 4 ) droplets are en- 
jected approximately 2 cm downstream of the back 
wall and 0.3 cm from the center line. The injection 
parameters are: 

Sauter mean diameters - SMD = 50 \i m 

Injection velocity - V[nj = 15 5 m/s for NjH.| 

V[ nj = 13.3 m/s for N 2 O 4 

Injector orientation - tilt — 30°, cone angle — 25 

Total injected propellant flow rate - 197. g/s 

Fuel oxidizer ratio - F/O = 0.83 

A A-squared distribution about the SMD is used 
for the sizes of injected droplets. Hence, the drop sizes 
are statistically sampled at the beginning of injection. 

Figure 7 presents the velocity field inside the 
chamber at two different times. It can be seen that the 
velocity field is mostly one directional away from the 
injection point. Figure 8 is a temporal evolution of the 
N 0 H 4 droplet parcel ballistics within the chamber. A 
high percentage of N 2 H 4 droplets can be seen imping¬ 
ing on the chamber wall and adhesing to it. Figure 9 
shows similar trends for the N 1 O 4 droplets. An imme¬ 
diate deduction from these results is that the droplet 
vaporization rates are grossly under-determined using 
the pure droplet evaporation (i.e. non-combusting) 
model. Figures 10 and 11 show the temporal evo¬ 
lution of the droplet parcel velocities for NMl., and 
N 2 O 4 , respectively. 

The temporal evolution of the temperature lickl 
is presented in Figure 12. A uniform region or high 
temperature (II) can be seen propogating downstream 
from the injection point where the temperature is rol- 
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atively low (L). The maximum temperature within 
the chamber is well below the adrabatic flame tem¬ 
perature due to incomplete droplet gasification. 

Figure 13 shows the fuel oxidizer ratio distribu¬ 
tion within the chamber for different times. Again, 
this points to the low droplet gasification rate. Mass 
fraction contours of nitrogen, which is the primary 
combustion product in this system, are presented in 
Figure 14 for different times. The N 2 concentration is 
seen to be maximum near the chamber exit tending 
towards an equilibrium value as time evolves. 


FUTURE D EVELOPMENTS 

The plan is to fully develop all of the aforemen¬ 
tioned submodels and efficiently incorporate them 
into the overall computational scheme. Some basic 
research will be undertaken specifically in the areas 
of drop/drop interactions and hypergolic bipropellant 
gasification modes. It is believed that the perfor¬ 
mance of this hypergolic system is strongly controlled 
more so by the gasification process than the gas phase 
chemical kinetic rates. 
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Decomposers 

Decomoos 1 tion Producis 



Figure 1. Schematic of pure monopropellant droplet combustion 



Figure 2. Schematic of hybrid monopropellant droplet combustion 
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Figure 4. 
various 


Theoretical and experimental burning rate constants with 
ambient temperatures for MMH monopropellant droplets 



Figure 5. Theoretical and experimental burning rate constants with 
various ambient temperatures for UDMH monopropellant droplets 


t =* 3.9 ms 




Figure 7. Gas velocity fields 



Figure 8. N 2 H, droplet trajectories 



Figure 6. The computing mesh 
(flow is from right to left) 


l. 0.74 ms 



t - 3.9 ms 



Figure 9. N 2 0 4 droplet trajectories 
































































Project Title: Spray Combustion Stability 

Principal Investigators: Dr. Pak-Yan Liang and Dr. S-M Jen; 
Other Investigators: Mr. Ronald Litchford 


Project Report 


PROJECT OBJEC TIVES 


The central purpose of this project is the im¬ 
provement of liquid-fueled rocket motor design tech¬ 
nology in order to assist the establishment of econom¬ 
ical commercial access to space through the develop¬ 
ment of engines with enhanced performance and relia¬ 
bility. Specific research effort in the project is focused 
on spray physics and associated combustion instabil¬ 
ity phenomena. Results garnered from this work wil 
contribute to the development of new computational 
took for design of stable liquid propellant rocket en- 

gines. 

The specific objectives of the research effort in¬ 
clude identifying and evaluating physical submodels 
which pertain to spray combustion stability with th 
idea of enhancing or refining existing submodels with 
a more comprehensive approach. In particular, any 
refinements to the spray combustion physical su - 
models which are achieved during the project will be 
channeled back to Rocketdyne for incorporation in 
their ARICC liquid rocket combustor code as secon- 
generation improvements. Also, as the ARIC 
forms the basis of future CFD development, some ef¬ 
fort is devoted to an evaluation of the code s capabi- 
ity for modeling oscillating pressure waves within the 

combustor. 

Over the past year, personnel associated with 
this project have changed significantly. Dr. Jeng as¬ 
sumed the role of principal investigator m May, an 
Mr. Litchford joined the project m June. The major 
research activities concerned with this project during 
1989-piay be summarized in the following three areas. 
(1) equation of state-of-the-art submodels for liquid 
rocket combustor instability study, (2) improvemen 
of droplet combustion model at elevated pressure, 3 
development of new, computationally efficient drop e 
tracking modek. 


^rnMPUSHMENTS 

Existing physical submodek have been evaluated 
with special emphasis on their capability for spray 
combustion instability study. The ARJCC code has 
ako been tested for transient spray and instability cal¬ 
culations. From these evaluations, submodel develop¬ 
ment work crucial to the establkhment of a predictive 
liquid rocket combustion instability code have been 
identified. It was realized that submodel improve¬ 
ments must be made before any CFD code can accu¬ 
rately calculate combustion instability, and a decision 
was made that initial work be devoted to physica 
submodel development. Three principal tasks were 
derived as the focus of a three year research effort. 

I - Hig h Bisssme Droplet 
fiasificn-tion Model 

The droplet vaporization rate k believed to be 
one of the most important parameters in determin¬ 
ing liquid rocket combustion instability. Furthermore, 
currently available models for rocket combustor cal¬ 
culations are based on derivations for low ambien 
pressure, whereas practical combustors operate at su¬ 
percritical chamber pressures. Because high-pressure 
thermodynamics equilibria calculations have shown a 
significant deviation in equilibria concentrations and 
heat of vaporization from low pressure assumptions, 
a high-pressure vaporization model is required to ac- 
predict the evaporation rate. Work toward 
development otan engineering model ..th the dcsned 
high-pressure characteristics has comprised much of 
the research effort during 1989 and b “" “‘' C “ E 
trated on LOX droplet vaporization within the SSi 

environment. 

- rimnlet T» rh»lent Dispersion Model 

Lagrangian tracking for droplet positions and as- 
sociatef properties has been recommended by rocket 
researchers as the best method for spray combusuo 
modeling. However, due to limited computer spe 
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and memory, a “computational droplet in a CFD 
code generally represents 10“* - 10 6 physical droplets. 

In addition, it is uneconomical and impractical in the 
engineering design practice to have a massive, statis¬ 
tically significant number of computational droplets 
in a CFD code. The reverse problem is that an insuf¬ 
ficient number of computational droplets will gener¬ 
ate computational noise and forbid spray combustion 
instability calculations. During the past year, initi¬ 
ation of a new model requiring a small number of 
computational droplets and capable of averting nu¬ 
merical difficulties has been made. The basic concept 
of the model has been established, and calculations 
have been compared against simple spray flows. 

Task III - Snrav Atomization Model 

Spray angle and droplet size distributions calcu¬ 
lated from the existing spray atomization models are 
based upon semiempirical correlations. The pressure 
and flow fluctuations upon the atomization process 
have not been taken into consideration, and this lim¬ 
its code capability for instability study. It has been 
recognized that a fully predictive method for spray 
atomization processes will not be mature within the 
next decade, and the intention in this study is to de¬ 
velop a transient spray atomization process based on 
a semi-empirical method. This task will commence 
during 1990. 

High Pressure Droplet Ciasification 

Research effort during 1989 has been concen¬ 
trated on evaluation of existing and newly formu¬ 
lated high-pressure droplet gasification models. A 
full report of this study can be found in Reference 
[1], Such attention is warranted by the role of pro¬ 
pellant vaporization as a rate limiting participant in 
instability mechanisms and by the fact that practi¬ 
cal liquid rocket combustors operate at supercritical 
mean chamber pressures. Furthermore, pressure sen¬ 
sitivity of the propellant vaporization rate must be 
correctly reflected by the mathematical submodel in 
response to acoustic pressure perturbations. Because 
the current droplet evaporation submodel in the AR- 
ICC code is a heuristic extension of conventional low- 
pressiire theory, reformulation of the model has been 
undertaken from first principles, following an exten¬ 
sive review of the literature to assimilate a basis for 
submodel development. 

A thorough study of previous work in this area 
has indicated that in approaching the modeling of 
droplet gasification at elevated pressures, the fun¬ 
damental physical mechanisms which distinguish the 


high-pressure case from the low-pressure case must, bo 
recognized and taken into consideration. The pecu¬ 
liarities encountered at elevated pressures derive pri¬ 
marily from thermodynamic considerations. Most no¬ 
tably, with increasing pressure absorption of ambient 
gas into the liquid occurs and real gas effects become 
important in determining thermophysical properties 
and thermodynamic vapor-liquid equilibria. There¬ 
fore, multi-component thermodynamics based on a 
realistic equation of state must be considered in the 
model formulation. Inclusion of realistic thermo¬ 
dynamics further implies the attainment of higher 
droplet temperatures which in turn demand that the 
surface regression effect and bulk thermal expansion 
be included in the transport model. The preceding 
considerations formed the basis for a high-pressure 
droplet vaporization model formulation. 

The fundamental approach used in the formu¬ 
lation is based on a quasi-steady gas phase coupled 
to a constant property liquid phase in time. Trans¬ 
port is modeled as a spherically symmetric, diffusion 
controlled process which can be corrected for convec¬ 
tive effects using steady-state heat and mass trans¬ 
fer correlations. To account for the high-pressure ef¬ 
fects, surface regression and bulk thermal expansion 
of the drop are considered, and high-pressure multi¬ 
component thermodynamics is employed to couple 
the liquid-phase and gas-phase solutions at the in¬ 
terface. 

The transport processes during gasification of a 
propellant droplet are depicted in Figure 1. The de¬ 
rived governing system of equations reduce to a set of 
ordering differential equations in time where the heat 
and mass transport coefficients are obtained from the 
simultaneous solution of two analytic equations [Ref¬ 
erence 1]. 

Because the droplet surface, under supercriti¬ 
cal ambient pressures, can be heated to the critical 
state where the surface tension and heat of vaporiza¬ 
tion tend to zero, finite rate supercritical gasification 
can occur. Furthermore, since the preceding derived 
transport model is conceptually valid only to the crit¬ 
ical state and due to the absence of definitive theory, 
the stripping vaporization criteria as proposed by Dr. 
Liang is tentatively retained. In this way, finite rate 
vaporization may be heuristically considered. 

Because the conventional low-pressure thermo¬ 
dynamic model is invalid at moderate and in¬ 
creasingly higher pressures, proper evaluation o 
vapor-liquid equilibria and enthalpy of vaporization 
demands rather rigorous formulation from mult.i- 
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component thermodynamics. From fundamental the- 
ory, the vapor-liquid equilibria constraints may be 
enforced such that the temperature, pressure, and 
component fugacities are identical across the phase 
boundary. If an equation of state is specified, the gen¬ 
eral thermodynamic relations for the fugacity coeffi¬ 
cients may be integrated to provide analytic equations 
for the equilibria concentrations in each phase. An 
application of the Redlich-Kwong equation of state 
with the mixing rules of Prausnitz and Chueh, [Re er- 
ences 2 and 3] known to yield good, verifiable results, 
is the conventional approach for droplet analyses at 
high pressure and has been effected in this study. In 
addition, the heat of vaporization may be determined 
from the rigorous thermodynamics relation between 
the partial molar enthalpies and the fugacity coeffi¬ 
cients for each phase. The resulting thermodynamic 
model may be partially simplified in special cases to 
neglect ambient gas absorption while retaining the 
real gas effects. 

The equilibria model was numerically solved for 
several hydrocarbon/inert gas mixtures for evaluation 
with respect to experimental data. Figure 2 shows the 
calculated results for a n-hexane/nitrogen binary sys¬ 
tem in comparison with experiment. Excellent Agree¬ 
ment with theory is obtained. Prediction of the low- 
pressure ideal-gas model is depicted by the dashed line 
to illustrate the associated error with increasing pres¬ 
sure. The calculated enthalpy of vaporization for the 
n-pentane/nitrogen binary system is shown in Figure 
3 indicating substantial deviation from the latent heat 
of vaporization with increasing pressure. 

Upon validation of the equilibria model, the for¬ 
mulated transport/thermodynamic governing equa¬ 
tions were combined and numerically solved in a sin¬ 
gle droplet vaporization code. The general numerical 
procedure adopted for the study was to evaluate ther¬ 
modynamic and thermophysical properties, heat and 
mass transfer coefficients, and vapor-liquid equilibria 
explicitly and assume these values hold over the nu¬ 
merical time step used in the solution of the governing 
ODE’s. Calculations were conducted for comparison 
with experiment and some results are shown in Fig¬ 
ures 4 and 5. Reasonable correlation with the data 
was obtained without adjusting any empirical con¬ 
stants. 

With the basic validity of the theory established, 
attention was turned toward modification °f the 
model to include the cryogenic propellants (0 2 /H 2 ) 
of primary interest. Since hydrogen is a quantum 
gas, the mixing rules were modified in the manner 


proposed by Chueh and Prausnitz to accommodate 
the desired propellant system. Equilibria calculations 
for the H 2 /0 2 system using the quantum gas mixing 
rules are given in Figures 6 and 7. Expected quali¬ 
tative trends are obtained. As an investigative pre¬ 
cursor to evaluation in a CFD code, the droplet va¬ 
porization code was then modified to a simplified 1- 
D rocket combustor code for the H 2 /0 2 propellant 
system. LOX droplets were assumed to be axially in¬ 
jected as discrete, identically sized particles in a cylin¬ 
drical chamber. The axial gas velocity was assumed 
to increase in proportion to the amount of mass evap¬ 
orated from the droplets, and the temperature and 
pressure were taken to be constant throughout the 
combustor. Calculations with the highpressure equi¬ 
libria model are compared with calculations using the 
current thermodynamic model in the ARICC code in 
Figures 8 and 9. with the principal result being that 
the high-pressure model predicts a much faster rate of 
development and a droplet lifetime roughly half that 
predicted by the currently used model. 

The high-pressure gasification model formulated 
for implementation into already computationally in¬ 
tensive CFD codes is obviously a rather sophisti¬ 
cated submodel with computational demands exceed¬ 
ing those encountered for the low-pressure models 
currently used. However, the required additional 
computing time is demonstrably not excessive and 
the ability to encompass the physics of high-pressure 
gasification is crucial in instability analyses where the 
vaporization response to pressure perturbations acts 
as the rate-controlling Rayleigh mechanism for acous- 
tic wave reinforcement. 

Direction concerning future refinement of the 
currently proposed high-pressure gasification model is 
focused on extension for supercritical gasification and 
will depend heavily upon experimental results avail¬ 
able for validation purposes. Currently, our thin ing 
is to incorporate an integral equation formulation for 
convective liquid-phase energy transport to obtain the 
internal surface temperature gradient so that finite- 
rate vaporization may be considered once the surface 
reaches criticality. Other possible improvements in¬ 
clude the addition of an envelope flame for conjugate 
supercritical combustion of LOX droplets and inclu¬ 
sion of H 2 0 in the equilibria calculation. 

p r ,->pW Turbu lent. Dispersion Modej 

The spatial distribution of propellant chemical 
energy release within a liquid rocket engine is primar 
ily controlled by the droplet life history - droplet tra¬ 
jectory and vaporization rate. For spray combustion 



modeling, ihe Lagrangian droplet tracking method 
has been most recommended- The stochastic sepa¬ 
rated flow (SSF) models which utilize the Lagrangian 
tracking method consider both slip and fluctuations 
relative to the gas flow and have been well calibrated 
for simple flows. 

In the SSF model, droplets are assumed to in¬ 
teract with a succession of turbulent eddies using 
random-walk calculations. For example, the droplet 
acceleration is calculated by 




where u^u^p^p^D are the droplet velocity, gas 
velocity, droplet diameter, droplet density, gas den¬ 
sity, and drag coefficient. The gas velocity is calcu¬ 
lated from its mean velocity plus random sampling 
from the fluctuating velocity. 

In extension of the SSF model for rocket combus¬ 
tion instability study, however, two major difficulties 
are encountered. First, if an SSF model is used for 
time-dependent, three dimensional spray calculations, 
a rather sufficient number of computational droplets 
are required to cover the droplet probability distribu¬ 
tion over nine independent variables - three spatial co¬ 
ordinates, three velocity components, droplet temper¬ 
ature, droplet size, and time. The random walk sam¬ 
pling over nine independent variables for each droplet 
property is still computationally too expensive or im- 
possible even with the aid of supercomputers. Sec¬ 
ond, since each droplet in the SSF model represents a 
Delta function over nine independent variables, “com¬ 
putational shot noise” will be generated if an insuf¬ 
ficient number of computational droplets are used in 
the code. 

Schematic representation of a typical spray into 
stagnant air is shown in Figure 10. If a large num¬ 
ber of computational droplets are used (number of 
droplets/number of computational cells > 100), the 
distribution of droplets per computational cell will be 
very smooth. However, if a small number of droplets 
are u$ed (number of droplets/number of computa¬ 
tional cells < 10) the simulated droplet distribution 
will be significantly altered from real flows. This zig¬ 
zag distribution represents spatially non-uniform en¬ 
ergy release and hinders successful combustion insta¬ 
bility calculations. 

Because of this inadequacy, it is desirable to 
modify the SSF model for time dependent spray 


combustion within rocket chambers. The modified 
model should require a small number of computa¬ 
tional droplets and should generate less computa¬ 
tional noise. Effort on this task began in August, and, 
thus far, the preliminary model has been established 
and calibrated against nonevaporating spray flows. 

The proposed model considers that a computa¬ 
tional droplet has a probability density function of 
temperature, size, number of physical droplets, and 
locations at different times. The shape of these dis¬ 
tributions over each variable will not be changed, al¬ 
though the variance will be altered. In the original 
SSF model, the random sampling over whole PDF 
distributions are needed. For the proposed method 
only the mean and variance of the PDF are calculated. 
This results in a significant savings in cpu time and 
memory requirements. The model has been tested 
against simple non-evaporating spray flows as shown 
in Figure 10. 

The calculations for the gas phase were per¬ 
formed using a modified version of GENMIX. The 
droplet trajectory was computed using a second order 
finite difference algorithm. Details of the SSF model 
may be found in Reference 4. 

Several test cases using from 100 to 8000 compu¬ 
tational droplets have been performed with 50 gas- 
phase computational cells radically. Figures 11-13 
show the radial distribution of the number of droplets 
per gas-phase cell calculated from the original SSF 
model and the modified SSF model. Results for the 
original SSF model using 200 and 600 droplets show 
significantly irregular distributions, whereas the mod¬ 
ified SSF model greatly reduced the noise. Both mod¬ 
els predict similar distributions when a large number 
of computational droplets (N p = 8000) were used. 
Clearly, the modified SSF model can predict a better 
spatial distribution when a small number of compu¬ 
tational droplets are used. 

Figure 14 shows the predictions of three different 
dispersion models. The original deterministic sepa¬ 
rated flow (DSF) model assumes that the gas-phase 
fluctuation velocity has no effect on the droplet trajec¬ 
tory and that each computational droplet represents 
a point source for a group of physical droplets. In 
the modified DSF model, droplet dispersion is con¬ 
sidered through the turbulence effect on the position 
variance. The results from both SSF models using 
8000 computational droplets to simulate the spray 
are also plotted in the figure. Both DSF models 
simulated the spray using 40 computational droplets 
which represents a significant reduction in the number 
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of droplets considered compared to the SSF model. 
Clearly, the DSF model without dispersion consid¬ 
erations underpredicted the distribution width. The 
DSF model with dispersion rate calculations can pre¬ 
dict the global radial distribution of droplets fairly 
well, although it did a poor job in calculating the lo 
cal structure. 
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pWQ.1F.CT O BJECTIVES 

The objective is to develop computational tech¬ 
niques for the design of high-area-ratio nozzles and 
to validate these models by comparison with experi¬ 
ments and computations using other codes. 

FVSQ AC COMPLISHMENTS 

Progress was made in two areas during the past 
year. First, performance computations were added to 
the PARC2D code and the performance of the Soi l 
nozzle was computed for inviscid, laminar and turbu¬ 
lent flow assuming a perfect gas with 7 - 
ond the PARC2D code was modified in a non-^AS 
project to compute equilibrium flow about hypersonic 
blunt bodies. Progress has been made toward mod¬ 
ifying this code to compute equilibrium H 2 /0 2 How 
through the SSME and related nozzles. 

First consider the perfect gas computations m 
the SSME nozzle. Last year it was reported that the 
PARC2D code in the Euler mode gave excellent agree¬ 
ment with inviscid solutions to the same flow com¬ 
puted by the Rockwell USA code and by the JAN- 
NAF method of characteristics code. This year lami¬ 
nar and turbulent solutions to the flow field were ob¬ 
tained and are reported in the thesis of Joseph L. 
Orr Jr. Several computational grids were used in¬ 
cluding one that clustered grids in the location of the 
embedded shock wave plus a tight clustering along 
the wall. The other grids had various wall pack¬ 
ings The results of the computation were found to 
be sofhewhat dependent upon the grid, especially for 
the boundary layer temperature distribution and the 
predicted adiabatic wall temperature. After this work 
had been completed it was discovered that other in¬ 
vestigators found that the first grid point must have 
v+ < 0 2 for the solution to be independent of the 
grid (J.S. Ryan, J. Flores and C.Y. Chow »Develop¬ 
ment and Validation of CNS {Compress bde Navier- 
Stokes) for Hypersonic Applications, AIAA-8J- Loaa, 


presented at AIAA 20th Fluid Dynam.cs Pl^rna Dy¬ 
namics and Laser Conference, June 12-14, 1989]. All 
future computations will use this criteria to ensure 
grid-independent solutions. Typical velocity and tem¬ 
perature boundary layer profiles are shown in Figures 

1 and 2. 

The solutions used a new artificial viscosity 
model for the PAR.C2D code, developed by person¬ 
nel at Sverdrup Technology, Inc., AEDC Group. This 
model decouples the artificial viscosity in the ‘ wo co¬ 
ordinate directions and reduces the fourth-order dis¬ 
sipation (based on the local cell Reynolds number) in 
the boundary layer to zero at the wall. Th.s new ar¬ 
tificial viscosity model is not nearly as robust as the 
original model and great care had to be taken to ge a 
converged solution. The turbulent solution was even 
more difficult to obtain and L2 decreases of only 
to 10 s could be achieved. 

Performance was added to the code by integrat¬ 
ing across the nozzle exit plane to get the mass flux 
and the thrust. Integration across the throat gave 
exactly the same value of mass flux. Because o the 
thin boundary layer in the SSME nozzle the pred.cted 
mass flux and specific impulse were very independent 
of the computational grid. The various grids are de¬ 
scribed in Table 1 and the performance results are 
given in Table 2. The specific impulse decrement Tor 
the laminar boundary layer was -1.4 sec an ■ 0 sec 
for the turbulent boundary layer, compared to the 
inviscid performance. 

The specific impulse for a perfect gas depends 
upon 7 M > the molecular weight. Examples o 

the sensitivity to these parameters are given 111 _ 

3 The predictions given in Table 2 assumed O/ _ 

8 so that only H 2 0 would appear alter combusuon. 
Thus M was assumed to be equal to 18. ° uev 
the actual SSME operates with 0/F = 6, ^mg an 
average molecular weight of 14 after complete co 

bustion. 1, Edition. 7 '«*' ?!“' 

The first part of Table 2 compares the change. . 
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pected in specific impulse if'/ — 1.2, M — 16 while the 
second assumes 7 = 1.18, M = 16. In both cases the 
specific impulse increased compared to the computed 
results for M = 18, 7 = 1.2 given in Table 2. Only 
an equilibrium code can predict the actual specific 
impulse. 

During the past two years the PARC2D code was 
modified to compute the flow field about blunt bodies 
moving at hypersonic speeds through the upper alti¬ 
tudes of the earth’s atmosphere. The thermodynamic 
properties are computed from statistical mechanics, 
incorporating the latest spectroscopic values to com¬ 
pute the partition function and including a significant 
number of excited electronic states. The mixture is 
assumed to be that of perfect gases because of the rel¬ 
atively low density, even after the bow shock, at the 
high altitude of interest. The composition at a given 
temperature and pressure is computed from the law 
of mass action, using statistical mechanics to compute 
the equilibrium constants. Diffusion is ignored except 
in the computation of the coefficient of thermal con¬ 
ductivity. 

The transport properties are obtained using the 
latest approximation to the exact kinetic theory re¬ 
lations and the latest computations of the collision 
integrals. The latter were obtained using accurate 
representations of the interaction potentials. There 
are three contributions to the coefficient of thermal 
conductivity of the mixture, from the transport of 
translational energy, the transport of internal energy 
and a reactive contribution caused by the diffusion 
that occurs in regions of dissociation, for example. 
This last contribution can be an order of magnitude 
greater than the other contributions to the thermal 
conductivity and can have a significant influence on 
the flow field. 

The equilibrium code was applied to the flow field 
about the Aeroassist Flight Experiment (AFE) vehi¬ 
cle which skips along the outer edge of the atmosphere 
at a Mach number of 33. The perfect gas tempera- 
ture was over 40,000 K after the shock wave but the 
predicted equilibrium temperature was around 8,000 
K. The stagnation point heat transfer was found to 
be significantly influenced by the real transport prop¬ 
erties, compared to the usual Sutherland approxima¬ 
tion. 


7,000 R, as compared to the JANNAF tables. Com¬ 
puted values of the specific heat are shown in Figure 
3. The composition for various O/F ratios, pressures 
and temperatures as well as the mixture properties 
can be accurately computed, when compared to the 
results tabulated in NASA SP-3011. The adiabatic 
flame temperature can also be accurately computed. 
As examples, the temperature and pressure computed 
along the centerline of the SSME nozzle by the per¬ 
fect gas code were used to compute the composition 
and ratio of specific heats for various locations in the 
nozzle. These results are shown in Figures 4 and 5. 
Although these results will not be repeated by the 
equilibrium code, they do indicate how the composi¬ 
tion and gamma will later be found to vary through 
the nozzle flowfield. 

Final modifications to the equilibrium version 
of PARC2D are presently underway and this new 
code will be used for various projects during the next 
year. These will include the SSME nozzle as well as 
the enlarged-throat SSME nozzle that was tested at 
NASA MSFC. The computed results for the latter 
nozzle will be compared with the measured exit flow- 
field. 

PUBLICATION 

Joseph L. Orr, "SSME Nozzle Flow Field Calcu¬ 
lations Using PARC CFD Code”, M.S. Thesis, The 
University of Tennessee, Knoxville, December 1989. 


This equilibrium version of PARC2D is presently 
being modified for H 2 /O 2 • At the present time the 
thermodynamic properties of the molecules present 
in H 2 /O 2 combustion, namely, H 2 , O 2 , H, O, OH, 
H 2 0 and H0 2 , can be accurately computed up to 
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A.I. PROPULSION APPLICATIONS 


• Identification and Interpretation of Patterns m Rocket Engine Data. Artific 
Intelligence and Neural Network Approaches 


• SSME Propellant Path Leak Detection 

. Intelligent Hypertext Manual Development for the Space Shuttle Hazardous Gas 
Detection System 

• Classification of Data Patterns Using an Autoassociative Neural Network Topol 

ogy 


SSME Component Assembly and Life Management Expert System 

Leader - An Integrated Engine Behavior and Desip Analyses Based on Real-Time 
Fault Diagnostic Expert System for Space Shuttle Mam Engine 
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INTRODUCTION. 


The Space Shuttle Main Engine (SSME) is a 
complex system whose health is monitored by a set 
of sensors placed at various points within the engine. 
While serious anomalies in engine behavior will ul¬ 
timately cause large fluctuations in the sensor data, 
these large fluctuations may not occur until it is too 
late to shut down the engine or to take other correc¬ 
tive actions. The purpose of the work reported here is 
to automatically detect anomalies in engine behavior 
as earlv as possible, when their indications in the sen¬ 
sor data are much subtler. Early detection of these 
anomalies is challenging because of the large amount 
of noise present in the data. In the presence of this 
noise, early indications of anomalies become very diffi¬ 
cult to distinguish from fluctuations in normal steady 
state operation due to the noise. 


Detection of anomalies is also challenging be¬ 
cause there is an extremely large set of possible fault 
conditions. Characteristics of only a fraction of these 
possible fault conditions are known, and only a very 
small fraction have ever been observed. 


This paper describes an expert system which is 
designed to perform automatic data analysis, iden¬ 
tify anomalous events and determine the characteris¬ 
tic features of these events. We have employed both 
artificial intelligence and neural net approaches in the 
design of this expert system. The artificial intelli¬ 
gence, approach is useful because it provides. 

1. The use of human experts’ knowledge of sensor 
behavior and faulty engine conditions in inter¬ 
preting data. 

2. The use of engine design knowledge and physi¬ 
cal sensor locations in establishing relationships 
among the events of multiple sensors. 


3 . The use of stored analysis of past data of faulty 
engine conditions. 

4 The use of knowledge-based reasoning in distin¬ 
guishing sensor failure from actual faults. 


The neural network approach appears promising 
for four reasons: 

1. Neural networks can be trained on extremely 
noisy data, and produce classifications which are 
more robust under noisy conditions than other 
classification techniques. 

2 . The neural network technique thus avoids the ne¬ 
cessity of noise removal by digital filtering, and 
therefore avoids the need to make assumptions 
about frequency bands or other signal character¬ 
istics of anomalous behavior. 

3. Neural networks can, in effect, generate their own 
feature detectors based on the characteristics of 
the sensor data used in training. This is in con¬ 
trast to other techniques which are based on the 
analysis of predetermined features of the sensor 
data. If anomalous engine behavior is indicated 
by some unanticipated feature present in the sen¬ 
sor data, a neural network is more likely than 
other techniques to detect this unanticipated fea¬ 
ture. 

4. Neural network algorithms are inherently paral¬ 
lel, and therefore are potentially implementable 
in special-purpose parallel hardware. Neural net¬ 
work algorithms therefore have the potential ot 
being useful for the detection of anomalous en¬ 
gine behavior quickly enough for suitable correc¬ 
tive action to be taken. 
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Architecture A 

The architecture of the Data Analysis Expert 
System (DAES) is shown in Fig. 1. The high fre¬ 
quencies of the noisy data are first filtered through a 
two-band digital filter so that the significant features 
of data corresponding to anomalous engine conditions 
are not lost. The upper and lower thresholds of the 
filter were determined by minimizing the number of 
segments in the filtered data. A segment is defined 
as the time duration between two consequent zero 
crossings in the first derivative of the filtered data. 
The Segmentation Module, shown in Fig. 2, identi¬ 
fies the segments. These segments are input to the 
Significant Event Selection Module (SESM) shown in 
Fig. 3. SESM compares the features of each input 
segment with the features of the most deviant seg¬ 
ment belonging to normal steady state sensor data 
of the corresponding sensor stored in the data base. 
The input segments which are more deviant than the 
segments of normal data are given as the significant 
events at the output of SESM. The event identifier, 
shown in Fig. 1, separates abnormal significant events 
from scheduled events by employing the knowledge 
about the scheduled events. This knowledge includes 
a planned event schedule and computer commands 
profile as well as specific features of scheduled events. 


Results of Prototype A 

A prototype called ‘PROTOTYPE A’ has been 
developed on the basis of the above described ar¬ 
chitecture A, shown in Fig. 1. !n order to illus¬ 
trate the results, we selected the analysis of Hr r 1 
(High Pressure Fuel Turbo Pump) failure in SSME 
test #901-364, specifically, hot gas intrusion to ro¬ 
tor cooling. The probable cause of the failure was 
a new HPFTP thermal shield retainer nut assembly 
used for the first time on this test. The geometry of 
the nut allowed a direct leak path through the heat 
shield for the high temperature gas which produced 
two jets impinging directly on the turbine end cap 
(Kaiser helmet) and reducing material properties in 
the impingement zone. Main Combustion Chamber 
Hot-Gas Injection Pressure (MCCHGP) sensor data 
from 320 to 400 seconds after the start of the test 
is shown in Fig. 4. The occurrence of the fault is 
noticeable from the global pattern which is super¬ 
imposed by the high frequencies. The filtered data 
is shown in Fig. 5 in which the global pattern is 
more visible. SESM identifies three significant events 
which are then further analyzed and their features 
are determined. These features are shown in Fig. 6. 
The identification of the significant events by DAES 


is compared with the NASA and Rocketdyne human 
expeVts’ analysis. DAES findings were identical to 
human experts’ 


Architecture B 


In Architecture A, digital filtering, segmentation, 
and significant event selection are performed by three 
separate modules, as shown in Figure 1. This method 
of identifying significant events is expensive computa¬ 
tionally, and requires filtering out high frequencies in 
order to extract a reasonable number of segments for 
further analvsis. The alternative is to replace these 
three modules with a neural network. Because of the 
noise resistance of neural networks, it is unnecessary 
to remove high frequencies by digital filtering. In¬ 
stead, the neural network accepts data in the form of 
a complete power spectrum for each successive win¬ 
dow of data. The purpose of the neural network is to 
identify significant events on the basis of these spec¬ 
tra, without further filtering or segmentation. Fig¬ 
ure 7 depicts the neural network architecture we are 
currently investigating as an alternative to the digi¬ 
tal filter*, segmentation, and significant event selection 
modules of Figure 1. Fast Fourier Transform (FFT, 
preprocessing is employed to produce data wlncn is 
time-shift invariant within each FFT window, but 
without filtering out the high frequencies. The result¬ 
ing spectral data is compared to the average power 
spectrum of normal steady-state data, and the differ¬ 
ence in power at each frequency is fed into the neu¬ 
ral network. The neural network is trained with this 
spectral data using Kohonen’s [1988] self-organizing 




Current tapes of SSME tests typically include not 
only sensor data, but also control information speci¬ 
fying the overall commands (such as thrust changes) 
given to the engine and the resultant commanded 
value settings. The overall sequence of steps in train¬ 
ing the network and classifying the data proceeds as 

follows: 

I If control information for the test is available, 
remove scheduled transients from the new data. 
If control information is not available, analyze 
the entire data stream. 

2. Perform Fast Fourier Transform of sequential 
windows of the time series test data. 

3 . If normal steady-state data for the same sensor 
is in the database, subtract the normal power 
spectrum from the power spectrum of the new 
data, yielding a delta (difference from normal) 
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power spectral density (P. S- D.) 

4 . Set the weights in the neural network to random 
initial values. 

5. Using Kohonen’s unsupervised learning algo¬ 
rithm to adjust weights within the network, se¬ 
quentially feed in the delta P. S. D.’s of successive 
windows of test data. 

6 . Turn off the learning algorithm. The weights in 
the network now represent the result of the learn- 
ing algorithm. 

7. Again feed in the delta P. S. D.’s of successive 
windows of test data. For each window of data, 
the output of the network reflects the neural net¬ 
work’s classification of that data. 

8 . Apply a threshold to the output activation level 
based on the observed outputs of normal steady- 
state data. The amount by which activity ex¬ 
ceeds this threshold represents the level of confi¬ 
dence that an event has occurred. 


Preliminary R e sults of Prototype £ 

Our initial implementation of the architecture 
Lbove has been operational for a few weeks, dur- 
ng which time we have evaluated its performance 
in SSME tests, for which both raw sensor data and 
inal NASA Investigative Board Reports were aval - 
ible. Our initial analysis was confined to each sen¬ 
ior meeting both of the following conditions: (i) The 
\jASA final report indicated an anomaly in the data 
for that sensor which occurred significantly before 
engine cutoff; and (ii) availability of data for that 
sensor. Results for the three sensors meeting these 
conditions are shown in Table 1. The original time- 
series sensor data (Figure 8), and plots of neural net- 
work activation obtained from FFT window sizes of 64 
data points (Figure 9) and 32 data points (Figure 10) 
are shown for one of the sensors, Main Combustion 
Chamber (MCC) Hot-Gas Injection Pressure The 
neural network activation plotted is the result of Step 
8 above, and represents the level of confidence that 
an event has occurred, as computed by the network 
baser? on the data with which it has been trained in 
each case. 


As shown in Table 1, and illustrated in Figures 
9 -10, the neural network detects events in agreement 
with the NASA Investigative Board analysis. Obvi¬ 
ously, a smaller FFT window size yields more precise 
timing of the events detected. On the other hand, on 
the basis of this pilot data, it appears that decreasing 


the FFT window size makes the anomaly detection 
process less resistant to steady-state noise. (Com¬ 
pare Figures 9 and 10). There would thus appear 
to be a trade-off between precision and reliability in 
the anomaly detection process. Our next step will 
be to investigate the effect on this trade-off of vari¬ 
ous alternatives in preprocessing steps 1-3 above, such 
as window shape, window overlap, and the metlio 
of comparing new data with known normal data in 
Step 3. While the initial results of tins technique ap¬ 
pear promising, much more work is needed before this 
neural network approach can be adequately compared 
with other approaches to identifying and classifying 
anomalies in SSME sensor data. 


rnNCLUSION 

The preliminary results obtained from both pro¬ 
totypes are in agreement with analyses done by N AS/ 
human experts. While much more work is needed .or 
a complete evaluation, it appears that these architec¬ 
tures are a promising first step in automating the ex¬ 
traction of features and subsequent diagnosis of rocket 
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Project Report 


PROJECT O BJECTIVES 


The complicated high-pressure cycle of SSME 
propellant path provides many opportunities for ex¬ 
ternal propellant path leaks while the engine is run¬ 
ning. This mode of engine failure may be detected 
and analyzed with sufficient speed to save critical en¬ 
gine test hardware from destruction. The leaks indi¬ 
cate hardware failures which will damage or destroy 
an engine if undetected, therefore, detection of both 
cryogenic and hot gas leaks is the objective of this 
investigation. 

The primary objective of this phase of the inves- 
ligation is the experimental validation of techniques 
for detecting and analyzing propellant path externa 
leaks which have a high probability of occurring on 
the SSME. The selection of candidate detection meth¬ 
ods requires a good analytic model for leak plumes 
which would develop from external leaks and an un¬ 
derstanding of radiation transfer through the leak 
plume. 

One advanced propellant path leak detection 
technique is obtained by using state-of-art technology 
of infrared (IR) thermal imaging systems combined 
with computer, digital image processing and expert 
systems for the engine protection. The feasibility of 
the IR leak plume detection will be evaluated on sub¬ 
scale simulated laboratory plumes to determine sen¬ 
sitivity, signal to noise, and general suitability for the 
application. The theoretical analysis was undertaken 
with the objective of developing and testing simple, 
easy-to-use models to predict the amount of radia¬ 
tion coming from a radiation source, background plate 
(BP), which can be absorbed, emitted and scattered 
bv the gas leaks. 

THEORECTICA l ACCOMPLISH MENT 

During this period the investigation concentrated 
on developing theoretical and numerical models to 
predict radiation transport in absorbing, emitting and 


scattering media. 

The transport equation describing the spectral 
intensity field is: 

^ ( 1 ) 

do 

This equation together with boundary condi¬ 
tion constitutes the basic mathematical problem to 
be solved for the description of intensity field. The 
source function I\{k\, w) is defined as. 


= (1 ~ ftoA K\fc(k\) 



w,wi)dwi 


( 2 ) 


and the albedo for scattering fio.\. defined as the ra¬ 
tio of the scattering coefficient, ct,a, to the extinction 
coefficient, K\ — 0jA + ®Ai 


Qoa = 


<r,\ 

K\ 


gVA 
0j,\ 4* 


where a\ is absorption coefficient. 

Equation (3) is the integrated form of the equation of 
transfer 


i\(ks) = *a(0) exp (-fc\) 

+ J w) exp [~(k\ - k]]dli\ (3) 


where k\ is a dummy variable of integration. 

Equation (3) is interpreted physically as the in¬ 
tensity being composed of two terms at optical depth 
k\. The first term is the attenuated incident radiation 
arriving at k\ (including, however, the contribution of 
induced emission along the path). The second term is 
the intensity at k\ resulting from spontaneous emis¬ 
sion and incoming scattering in the S direction y 
all thickness elements along the path, reduced by ex¬ 
ponential attenuation between each point of emission 
and incoming scattering k\ and the location k\. 
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The radiation transport calculations have been 
done for a non-gray and non-isothermal gas ( HnO ). 
Water vapor was selected as the molecule to concen¬ 
trate detection on due to the water vapor present 
in the hot hydrogen propellant lines. The radiation 
transport equation in the absorbing and emitting me¬ 
dia with a non-gray and non-isothermal gas is much 
more complicated in cylindrical geometry than paral¬ 
lel plate case which can be found in literature. These 
models will be used to evaluate the intensity of back¬ 
ground and leak plume radiation reaching the IR de¬ 
tector. The current radiation transport model for a 
simulated leak plume in form of cylindrical geometry 
has been developed (Figure 1). Numerical calculation 
of directional intensity passing through this cylindri¬ 
cal plume of water vapor at each y location (see Fig¬ 
ure 2) has been completed. Figure 3 and 4 are the 
typical numerical variations of intensity at different 
y locations. Figure 3 indicates that the intensity has 
maximum values at y = 0 and minimum values at 
y = R, which has the same pattern for each value 
of wave length but the peak values appear at differ¬ 
ent intensity due to the absorption coefficient at that 
particular wavelength. Figure 4 shows the minimum 
intensity at y = 0 and maximum value at y = R, for 
the case of hot background plate and cold gas plume. 
Most of the incoming intensity will be absorbed by 
cold gas at the center line, where y = 0 and 5 is 
maximum. When y values increases, the S value de¬ 
creases, therefore, there is not enough distance for the 
incoming energy to be absorbed. Also, closer to the 
outer edge of the plume the temperature decreases 
with increasing r. 

The variation of intensity with wavelength is 
shown in Figure 5 and 6. For case of hot plate with 
cold plume Figure 5 indicates that most of the incom¬ 
ing energy has been absorbed between 2.4 and 2.9 /x 
for different y values, but for case of cold plate with 
hot gas, Figure 6, shows that most of energy has been 
emitted by the hot gas. 

EXPERIMENTA L RESULTS 

The Thermovision 870 system operating in the 
2-5 micron band was selected for the initial experi¬ 
ments' The system can measure 10,000 point array 
at a rate of 25 times a second. The optical system 
with an extensive choice of lenses and filters enables 
measurement of a large range of temperatures -20 C 
to 1500°C over various distances. The availability and 
successful utilization of several Thermovision 870 sys¬ 
tems at the U.S. Air Force Arnold Engineering Devel¬ 
opment Center, aided the selection process. A formal 


loan request from CASP to AEDC lias been approved, 
and the IR instrumentation was made available in 
July 1989. 

The feasibility experiments have been completed 
using the Thermovision system to obtain images of 
both hot and cold plumes in front of a warm back¬ 
ground plate (350°F). The cold plume was provided 
by a room temperature mixture of air, water vapor 
and 2fi water droplets. Figure 7 shows a distinct 
plume against the warm grey background plate as ob¬ 
tained by the IR camera and enhanced by the VICOM 
digital image processing system. These results have 
confirmed the sensitivity of the Thermovision system 
to low intensity signals. 

P lanning is being conducted for more quantita¬ 
tive experiments with leak plumes which closely sim¬ 
ulate SSME propellant path properties. These exper¬ 
imental results will be used to validate the radiation 
transport model. 



F1fur« I. SchtMtic Mlutiritlon of Up*rt«ef»ul S*t-up 



Figure 2 . Plu*e Cross Section 
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for the Space Shuttle Hazardous Gas 
Detection System 
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Project Report 


This research is designed to utilize artificial in¬ 
telligence (AI) technology to increase the efficiency of 
personnel involved with monitoring the Space Shut¬ 
tle hazardous gas detection systems at the Marshall 
Space Flight Center. The objective is to create a com¬ 
puterized service manual in the form of a hypertext 
and expert system which stores experts’ knowledge 
and experience. The resulting Intelligent Manual will 
assist the user in interpreting data timely, in identi¬ 
fying possible faults, in locating the applicable doc¬ 
umentation efficiently, in training inexperienced per¬ 
sonnel effectively and updating the manual frequently 
as required. 


INTRODUCTION 

The complexity of space shuttle launch opera- 
tions is well recognized. Monitoring the shuttle main 
propulsion system for abnormal gas leakage during 
the launch phase is one method used to detect sub¬ 
system faults. In general, Operation and Main ten Mice 
Requirements and Specification Document (OMRSD) 
information and Operation, Maintenance and Instruc¬ 
tion (OMI) procedures have been developed to at¬ 
tempt, on a very limited basis, to isolate a leak prior 
to scrubbing a launch. Once a launch is scrubbed, 
more detailed isolation procedures are in place to iso¬ 
late a leak prior to detanking. 

The proposed Hazardous Gas Intelligent Hyper¬ 
text Manual (HIHM) will provide a significantly im¬ 
proved mechanism for the Hazardous Gas Detection 
System (HGDS). The proposed system will furnish 
“quick and easy” access to appropriate information 
for fault detection, data interpretation and action op¬ 
tions for the Huntsville Operations Support Center 
(HOSC). It may display detailed Mass Spectrome- 
ter information and Main Propulsion System (MPS) 
information for the Console Engineer and Manager. 
This information might include the hardware con¬ 


figuration, the sensor port locations, the gas species 
“real” concentration, etc. 

Additional Application Scenarios. The Intelli¬ 
gent Hypertext Manual (IHM) offers a system for use 
in training personnel through the large volume of ma¬ 
terial available for a specific area of space shuttle. The 
management of documentation for the HGDS " i a so 
be greatly simplified through easy maintenance and 
update capabilities. 

Benefits. The Intelligent Hypertext Manual will 
provide accurate service information in a rapid and 
efficient manner. It will eliminate most human er¬ 
rors and reduce the required resources including man¬ 
power and material cost to achieve the same speci¬ 
fied goal. In short, the system will assist to improve 
productivity and quality of service for fault detection 
from the HGDS. 


intf.t.ugf .nt HYPERTEXT MAHIIAL 
rnw-nmiOAT. DESCRIPTION 

A H vnszl&xt Manual 


The traditional printed service manual restricts 
laterial to a linearl arrangement. In searching for 
simple procedure, one must constantly make re- 
•ences to other chapters or sections, index entries, 
lotnotes, figures, tables, etc. Often, a text invites 
ie user to skip a section if he is at a different tech- 
ical competence level. This reference trail among 
aper documents requires users to exert substantia 
hysical effort and time to gather the complete pro 
•dure information. The present proposal offers a 
Dlution by converting the linear printed text into 

“hypertext/hypermedia” system supported by cur- 
ent computer technology. Hypertext (Ref. 1 )ms an 
dvanced document management system in which a 
ser can manipulate a database via graphical objects 
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on a computer screen that represent objects inside the 
computer memory. 

The basic feature of the computer-based hyper¬ 
text'is that windows on the screen are associated with 
objects in a database, and links are provided between 
these objects graphically as well as in the database. 
Usually, screen objects are manipulated with a mouse, 
and several blocks of documentation can be displayed 
simultaneously in different windows as shown in Fig. 
1. For example, specific documents of the Hazardous 
Gas Detection System for the Space Shuttle could be 
displayed in a window and at the same time the de¬ 
tailed schematic for the MPS flight measurement lo¬ 
cations could be displayed in another window simulta¬ 
neously. The on-line customized document can be or¬ 
ganized to meet users’ current requirements. In short, 
the hypertext manual has the advantages of tracing 
references easily, of creating new references, of creat¬ 
ing customized information display, and modularity 
of information. 

Fortunately, it is not required to design and pro¬ 
gram a hypertext system from the basic AI program¬ 
ming environment, since several commercial systems 
are readily available (Ref.l), e.g. KMS (Knowledge 
Management system Ref. 2) which runs on the SUN 
workstation and enables users to create a hypertext 
of a knowledge base. The software package Guide 2 
(Ref.3) runs on the IBM Personal Computer or com¬ 
patible or the Macintosh computer and has a com¬ 
plete development and management system for large 
volumes of documentation. The Hypertext technol¬ 
ogy has been given much attention recently. In a 
special issue of the ACM Journal(Ref. 4), authors 
presented six papers at HYPERTEXT ’87 (Univer¬ 
sity of North Carolina). Also, an in-depth review of 
Hypertext appeared in BYTE (Ref. 5). In addition, it 
should be possible to extend hypertext to hypermedia 
which will enable the user to capture knowledge in a 
specific domain with graphics and other media. Both 
Guide 2 and KMS offer graphic display capability. 


system will provide the appropriate recommendation, 
comment, and advice in response to the user s answers 
via a search algorithm in the system. For example, 
the user inputs to the computer the gas sample in¬ 
formation from the mass spectrometer. The expert 
system will recommend a corrective action and sup¬ 
porting information if the data indicate a problem, 
otherwise the expert system will indicate no action is 
required but may display a comparison with similar 
data observed in the past. Recommendations for cor¬ 
rective actions will be based on experts’ stored opin- 
ions and historical data. 

The organization and representation of the 
knowledge base for the expert system will be an¬ 
other key effort in the project. A model will 
be selected from the available paradigms which in¬ 
clude frame-rule representations and objective ori¬ 
ented techniques. Interleaving of the expert advi¬ 
sory knowledge base and the computerized manual 
will create an Intelligent Hypertext Manual. 

Q Intelligent Rvpert.ext Manual 

To couple an expert system with a computerized 
manual is an idea recently proposed independently by 
Lo Sc Steinle (Ref 6) and Williamson (Ref <). This 
system is contrasted with the traditional expert sys¬ 
tem which usually does not include a large volume 
of information. The proposed system will place in 
memory a large volume of text /graphics documenta¬ 
tion, such as, schematics, optional actions, and service 
procedures from the hypertext database, the archi¬ 
tecture for Intelligent hypertext manual is shown in 
Fig. 2 with the selected softwares. 

The effort under this project concentrates specif¬ 
ically on the engineering domain of the subsystems. 
However, the underlying software architecture and 
the other computer technologies will be applicable to 
other domains, such as engineering management or 
intelligent tutoring. 


fi T Rule-Rased Ex pert Systems 

lit order to more effectively utilize the comput¬ 
erized manual a rule-based expert system will be ^cre¬ 
ated. For selected systems of the shuttle, experts’ ex¬ 
perience and knowledge will be stored in the form of 
rules. Expert systems can typically explain their rea¬ 
soning and offer a user-friendly interface. A typical 
consultation session consists of a series of dialogues 
between the expert system and the user. The user 
answers questions asked by the system. The expert 


SOFTWARE AN D HARDWARE 

rnNSinF,RATIONS 

A Hypertext S oftware Evaluation 

The examination of commercial Hypertext pro¬ 
grams resulted in the purchase of two programs ^ 
'valuation purposes. KMS operates on* he J“ n f 38 ®' 
workstation. The hypertext system, Guide 2 (Ret. J). 
runs on an IBM personal computer. 
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A demonstration of hypertext systems using 
KMS and Guide 2 was prepared using the current 
HGDS IHM project plan as its data base. The demon¬ 
stration package takes advantage of all the hypertext 
database features incorporated in KMS and Guide 2, 
including text and graphics. 

KMS offers a programming language for au¬ 
tomation of repetitive actions, which is equivalent to 
macros. 


Three features of Guide 2 are judged superior 
to that of KMS. Guide 2 runs on IBM PC com¬ 
patible computers which makes hardware acquisition, 
software installation and portability better than a 
workstation-based system. Guide 2 is currently in¬ 
stalled on a COMPAQ portable computer which is 
for demonstrations away from the CASP office. Also, 
Guide 2 allows more than two windows, albeit smaller 
than a single KMS window, of information to exist 
on the computer screen simultaneously. Guide 2 pro¬ 
vides three types of links; notes, expansion buttons 
and references expansion to the database as oppose 
to the single type of KMS link frames. The three links 
in Guide 2 offer the user greater flexibility than the 
KMS frames. 


KMS offers more powerful features than Guide 
2 for formatting and displaying information on the 
monitor. KMS displays two full pages of information 
simultaneously as opposed to Guide’s several partia 
pages of information. The KMS screen is a WYSI¬ 
WYG (what you see is what you get) display tha 
offers better resolution. 


Both KMS and Guide 2 offer limited graphic im¬ 
age editing capabilities. With Guide 2 the user can 
resize images in two dimensions. Guide 2 also allows 
moving images around the document. The user cuts 
images to the Microsoft Windows Clipboard for past¬ 
ing into new locations. KMS gives the hypertext user 
the same graphics editing features as Guide 2. Addi¬ 
tionally, KMS contains drawing tools to include lines, 
rectangles, polygons, ellipses and arcs in a document. 
Completely closed figures such as a circle can also be 
filled with varying levels of shading. 


KMS and Guide 2 each provide only one method 
for importing graphics into a hypertext document. 
Guide 2 uses the Windows Clipboard program for 
pasting images from external programs into a doc¬ 
ument. Therefore the computer system used to build 
a hypertext database must have programs capable o 
cutting or pasting to the Clipboard. The Clipboard 
limits the amount of an image that can be pasted into 


Guide. A sample A4-page-size image had to lie cut 
to the Clipboard in five individual pieces for pasting 
into Guide 2. KMS imports images using a single 
command that is attached to a point in the docu¬ 
ment. The command inserts the entire image file into 
the KMS document. The size of the file is not limited 
to the display. One image tested occupied two full A t 
pages. 

Although, KMS offers only one type of link, its 
links can be shared over directories and even net¬ 
works. The access performance of KMS on the Sun 
386i, about one second to fill a screen with new data, 
exceeds that of Guide 2 running on a 386 PC. 


KMS offers the user powerful indexing capability 
for rapid navigation through the hypertext database. 
The index feature gives the hypertext manual a def¬ 
inite advantage over the IBM-PC based Guide 2 hy¬ 
pertext. However KMS does not offer pop-up note 
capability like Guide 2. Advantages of KMS is the 
ability of the SUN to run several programs simulta- 
neouslv. 


Guide 2 only runs under the Microsoft Windows 
Environment using DOS 2.0 or higher as an operating 
jvstem which is subjected to the RAM 640k limit. 
On the other hand, KMS runs under standard UNIX 
aperation system directly and has no limitations as 


p Expert System So ft ware Evaluatio n 

VP-Expert and EXSYS are DOS compatible pro¬ 
grams that offer rule based artificial intelligence ca¬ 
pabilities. The inference engine in VP-Expert and 
EXSYS consists of both backward and forward chain- 
ing with confidence factors assigned to each knowl¬ 
edge based production rule. The shell is also compat¬ 
ible with spreadsheet, database and ASCII text exter¬ 
nal programs and is able to exchange data with these 
files. The executable programs such as filename.exe, 
filename.com, or filename.bat files can be accesse 
from the shell. The intelligent knowledge bases are 
constructed from these external programs in addition 
to knowledge-based rules acquired from domain ex- 




r Integration of Hy p ertext and Expert gvstern s 

Calls to/from the hypertext document must be 
made to/from an expert system for full implemen¬ 
tation of the intelligent hypertext manual. Guue - 
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offers a "launch" capability in conjunction with com¬ 
mand buttons to interface external programs. How¬ 
ever, when Guide 2 launches an external program, 
that program must also be compatible with the Mi¬ 
crosoft Windows environment. In one application a 
reverse process worked successfully using VP Expert 
and Guide 2. Instead of Guide 2 making an external 
reference, VP Expert was used to call Guide 2 when 
the expert system user wanted to search the hypertext 
database. 

KMS, when operated under the SUN View en¬ 
vironment can easily interface with a DOS window 
that contains an expert system. To switch between 
DOS and KMS the user establishes two windows on 
the screen and simply points the mouse cursor to the 
window of choice. Unix, which is the operatmg sys¬ 
tem on the SUN, also provides a powerful “piping" 
capability for passing information from one applica¬ 
tion to another. The evaluation of the piping feature 
is now beginning. 


D Database Mass Storage 

In anticipation of the large volumes of documen¬ 
tation that will have to be assimilated, the mass- 
storage technology for text as well as for graphic 
images is being investigated. Text requires approx¬ 
imately 1 Kbyte of storage per page. For graphics at 
a low pixel density, the typical file of compresse lm 
ages requires 15 to 30 Kbytes. Magnetic disk storage 
may satisfy the requirement for the prototype intel¬ 
ligent Manual. However, it is a key objective of the 
proposed research to evaluate mass-storage technol¬ 
ogy, such as Erasable Optical,CD-ROM and WORM, 
for the future extension of the project. This tech¬ 
nology is also required to implement the system on 
an information server computer. The CD-ROM is a 
technology similar to digital audio compact disk (CD) 
used for the storage and retrieval of digitally encoded 
(machine-readable) information. For mass informs 
tion storage, CD-ROM is a reliable and permanent 
medium. The memory capacity of one CD-ROM is 
600 Mbytes which is equivalent to 1500 (5.25-inch) 
floppy disks. One disadvantage is that CD-ROM has 
longer, access time and a slower data transfer rate 
compared to a magnetic disk. WORM is another 
technology which is just like CD-ROM but allows the 
user to “Write Once and Read Many” times. The 
CD-ROM/WORM technology will be evaluated for 
storing the large database required for an intelligent 
hypertext manual. 

The new erasable optical disks offer another al¬ 


ternative to storage of large databases. Current first 
veneration systems store 300 megabytes of data per 
side per disk. Costs range from $4700 to $6500 per 
unit This technology is new and each vendor offers 
unique formats so data is not interchangeable between 
machines. (Ref 13). Due to its rewrite capability, 
the erasable optical drive offers significant advantages 
over both the WORM and the CD-ROM disks. 


Disk storage requirements for text and graphical 
files were established. Text data in ASCII formats 
occupies 1 bit per character. The different hypertext 
formats were also evaluated and it was determined 
that text storage requirements will not be a problem 
for the Intelligent Hypertext Manual. Three types 
of graphical files were evaluated for size requirements 
and included the SUN graphics image format, PC- 
Paintbrush format and Microsoft Windows format. 
These files can occupy up to 1 Megabyte of storage 
for an eight by ten inch figure. The UNIX operat¬ 
ing system on the SUN offers a compress program 
which gave significant reductions in graphic file stor¬ 
age needs. A compress program was found for the 


F Ont.ical Scanners 


Requirements. An Intelligent Hypertext Manual 
will consist of large volumes of textual and graphical 
information which must be available, on-line, to the 
user at the HOSC. To meet the goal of demonstrating 
a hypertext manual of the HIHM, the assimilation of 
the many volumes of information linked to the HGDb 
arrnmolished rapidly. 


The use of electronic, desktop scanners for input 
of information to the Space Shuttle HIHM will expe¬ 
dite the task of making required systems information 
available to the expert system. Scanners alleviate the 
need for information to be retyped, redrawn or repho 
tographed. They also assure greater accuracy than 
manual input through a keyboard and/or mouse. 


Selection Criteria. Scanners perform two func- 
>ns, graphical input and optical character recogm- 

>n (OCR). The scanner for use with the HIM mu 
•rform both functions with a high 
yr and accuracy. The manuals for the IIGDS contain 
xt, figures and photographs which the user will need 

the IHM. 

Although the hardware for scanners is similar the 
uality of scans and the capability of OCR determine 
ie usefulness of particular systems. Scanners accom 
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plish both tasks with software, and usually it is the 
software that determines the capabilities of a given 
piece of hardware. OCR uses software that reads type 
and stores it in ASCII, or other, character formats. 
Graphics are stored using bit-mapped formats that 
can be manipulated using paint or similar programs. 
Most scanners require separate scans for each type o 
input if a document contains both text and grap - 
ics. The scanner must possess the capability to read 
multiple styles of type for input into the hypertext 
database. 

The scanner may need to handle the bulky and 
bound books of the Operation and Maintenance Re¬ 
quirements and Specification Document (OMRiiD) 
and Operation, Maintenance and Instruction (OMI) 
procedures volumes. The sheet-fed scanners require 
that bound documents be taken apart before input. 
The compromise between sheet-fed scanners with au¬ 
tomatic document feed and scanners that handle 
bulky documents will become one of time and avail¬ 
ability. If extra copies of the documents are available 
for breaking down into single sheets, a high capacity 
ADF machine may be most efficient. On the other 
hand, if the documents must be rebound a flatbed or 
overhead scanner may be the most efficient, depend- 
ing upon the speed of the scanner. 


S T 1 MMi* RV nF ACCQMPTiISffMEHl 

The objective of the present Phase I is to estab¬ 
lish the feasibility and merit of the proposed system. 
An intelligent hypertext manual will be created to 
demonstrate the feasibility of a selected subsystem 
with a limited knowledge base. This hypertext man¬ 
ual will be demonstrated to NASA personnel at the 
completion of Phase I. NASA’s comments and recom¬ 
mendations will be incorporated into the final design 
of the hypertext manual. The capability of the expert 
system’s inference mechanism will be illustrated by a 
few production rules and facts. 

As the first step, it was initiated to evaluate hy¬ 
pertext and expert systems software, mass-storage 
technologies, scanners and host computers. Results 
of th*e evaluations are included in this report. A 
demonstration System was constructed using both 
KMS and Guide-2 software. The integration of the 
expert system shell EXSYS into the hypertext sys- 
terns is in progress. 


FUTURE WORK 

In the future work, the development of the expert 
systems for a selected knowledge base will be com¬ 
pleted aided by NASA’s domain experts. The expert 
system alone will be demonstrated separately before 
integrating with the hypertext manual. The major 
task that will be accomplished during this phase will 
be the integration of the hypertext manual and expert 
system into the prototype Intelligent Hypertext Man¬ 
ual A demonstration of the system to NASA person¬ 
nel at this stage will be crucial for the implementation 
of all improvements recommended by NASA. The ac¬ 
quisition of the knowledge base for the HIHM will 
continue in the selected subsystems. A critical eval¬ 
uation of HIHM will be performed by NASA in the 
later part of the work, after assembling and incor¬ 
porating the improvements recommended by NAS/ 
personnel. 
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Abstract 

An effort is underway at the University of Ten¬ 
nessee Space Institute to develop diagnostic expert 
system methodologies based on the analysis of pat¬ 
terns of behavior of physical mechanisms. In this 
approach, fault diagnosis is conceptualized as the 
mapping or association of patterns of sensor a 
patterns representing fault conditions. Neural net- 
works are being investigated as a means of storing 
and retrieving fault scenarios. Neural networks 
offer several powerful features useful in fault diag¬ 
nosis, including 1) general pattern mapping 
capabilities, 2) resistance to noisy input data, 3) the 
ability to be trained by example, and 4) the potential 
for implementation on parallel computer architec- 
tures. 


Results of a prototype neural-network-based 
diagnostic system applied to jet and rocket en &™ 3 
are documented in Refs. 1 and 2. Although the 
performance of the earlier diagnostic system was 
quite good, some desirable improvements were iden¬ 
tified. The most important improvement deals with 
how the networks respond to fault classes for which 
the networks have not been specifically trained. 
(These fault classes are referred to as unknown 
classes in this paper). The networks used in Reft. I 
and 2 were trained to categorize inputs representing 
known fault conditions. However, there was no 
method whereby a network could be trained to 
recognize when an input pattern corresponded to a 
previously unobserved, i.e. unknown, fault condi¬ 
tion. When presented with sensor patterns repre¬ 
senting an unknown behavioral scenario, the 
networks would sometimes incorrectly categorize 
the patterns as a known fault condition. As a result, 
the diagnostic system rarely produced false negative 
results (i.e., reporting no fault condition when given 
a set of patterns indicative of a known fault condi¬ 
tion), but would sometimes produce false positive 
results when given input patterns representing a 
behavioral scenario for which the networks were not 
specifically trained. 


A neural network topology is presented which 
effectively addresses these problems. The network 
topology is autoassociative, i.e. the network inpu 
and output is identical when properly trained, and 
hence learning is unsupervised. Classification of 
input patterns is accomplished by examining the 
activations of the nodes in the hidden layer. A range 
of training curves indicative of a particular fault 
class will result in hidden node activations which 
map out a surface in N-dimensional space, where N 
ft the number of hidden nodes. An input curve 
presented to the network for classification will in¬ 
voke hidden node responses which can be concep¬ 
tualized as a point in N-space. If the hidden node 
activation lies near a surface corresponding to a 
particular fault, the input pattern will be associated 
with that fault If, on the other hand, the resultant 
hidden node activation lies away from any surfaces 
defined by the training set, the input will be clas- 
sified as an unknown pattern. Therefore, 
pattern produces a hidden node response which is 
sufficient to either 1) associate the input with a 
known fault scenario, or 2) reject the input - 
together as not representative of any known fault 

condition forwhich the netwcnrk has been trained. H 

an input is associated with a fault condition, the 
hidden node activations will also serve to identify 
the severity and duration (i.e., rapidity of fault 
onset) of the fault. 

Another improvement to the system reported in 
Refs. 1 and 2 concerns the number of networks 
required for diagnosis. The earlier system required 
five networks for each sensor, one to diagnose the 
fault type, and two additional networks to identify 
fault severity and duration for each fault. Two main 
faults were considered; as a result, the entire system 
required 20 networks. Adding another fault type to 
the system's diagnostic capability required the 
training of twelve additional networks. 

With the new approach, a single networkcanbe 
trained to respond accurately to input curves from 
all sensors; as a result, a single network performs 
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the functions of all 20 networks required in the 
earlier system. In the present approach, a single 
networkis trained to identify and categorize all fault 
conditions. As a result, adding the capability of 
identifying additional faults requires the training of 
only a single network. 

This paper describes 1) the network topology 
used in the present implementation of the diagnos¬ 
tic system, 2) the training regimen used, and 3) the 
response of the system to inputs representing both 
previously observed and unknown fault scenarios. 
The effects of noise on the integrity of the diagnosis 
are also evaluated. 

Introduction 

The University of Tennessee Space Institute 
(UTSI) is currently involved in the development of 
advanced methodologies for the detection and diag¬ 
nosis of faults in aircraft and space systems. Faults 
in these domains often develop rapidly, requiring 
diagnostic systems to respond in near real-time. 
Also, these problem domains require diagnostic sys¬ 
tems to exhibit resiliency in the event of noisy data 
and sensor degradation or failure. 



Figure 2. Mapping Of Sensor Patterns To Fault Conditions 


Sensor data typical of a jet engine fault are 
depicted in Fig. 1. Four sensors, exhaust gas 
temperature (EGT), combustion temperature (CT), 
and the rotational speeds of the low and high pres¬ 
sure compressors (N1 and N2, respectively) are 
monitored during engine operation. The data are 
normalized with steady-state values for con¬ 
venience. In the present effort, it is assumed that 
the engine is operating at some steady-state cona¬ 
tion when the fault condition occurs. The fault will 
be manifested as a deviation of sensor values from 
the steady-state condition. In the present effort, a 
time window containing 40 sensor readings span¬ 
ning four seconds is used. 

It is not sufficient for the purposes of diagnosis 
to simply detect when and whether a deviation from 
steady-state conditions has occurred; how the devia¬ 
tion is manifested over time is also important. For 
instance, an observation that a particular sensor 
parameter is decreasing linearly will likely result in 
a different diagnosis than that obtained from ob¬ 
serving an asymptotic decrease (Ref 3). Therefore, 
to be effective, a diagnostic system must be respon¬ 
sive to the qualitative (as well as the quantitative) 
behavior of the engine. The diagnostic process must 
also exhibit resilience to noise. A noise-corrupted 
version of the fault is depicted in Fig. 1 for com¬ 
parison. The 6% noise level means that the standard 
deviation for a Gaussian distribution of perturba¬ 
tions about the noise-free curves is 0.05. A 5% noise 
level would be considered excessive for most in¬ 
strumentation; however, satisfactory operation at 
this noise level is used as a goal in the present effort. 


Many fault diagnosis problems may be concep¬ 
tualized as the mapping or association of temporal 
patterns representing sensor data to patterns rep¬ 
resenting fault conditions. This process is il¬ 
lustrated in Fig. 2, which depicts a mapping process 


52 












from patterns of incoming sensor data to a pattern 
of real numbers which can be associated with a 
particular fault condition. Each sensor inputs a pat¬ 
tern into a classifying system, which attempts to 
determine which fault condition is indicated by that 
sensor. In the example of Fig. 2, sensor SI has been 
classified as indicative of either normal operation or 
a bearing failure. This ambiguity may occur if the 
bearing failure is of a low severity, or has just 
started to manifest itself in the sensor patterns. 
Also, sensor SI may not be highly responsive to 
bearing failures. Sensor S2 has more definitely in¬ 
dicated a bearing failure, while sensor S3 is unable 
to m alf« a positive diagnosis of any known failure. 
The ambiguities are resolved by an arbitrator, 
which con tains domain knowledge about the diag¬ 
nostic process, such as the sensitivity of certain 
sensors to different faults. Using this knowledge, 
the final diagnosis of bearing failure is performed 
correctly. 

The conceptualization of the diagnostic process 
as a pattern association or mapping process leads 
naturally to the investigation of pattern mapping 
and classifying methodologies which are applicable 
to real-time fault diagnosis. Neural networks ex¬ 
hibit several features useful in diagnostic systems, 
such as 1) general mapping capabilities, 2) resis¬ 
tance to noisy input data, 3) the ability to be trained 
by example, and 4) the potential for implementation 
on parallel computer architectures. 

A system based on neural networks to diagnose 
the development of fault conditions in jet engines 
was designed and developed at UTSI (Refs. 1 and 2). 
The system utilized 20 neural networks arranged in 
a hierarchical fashion which 1) input unprocessed 
data from four sensors, 2) identified the fault condi¬ 
tion, and 3) estimated the severity and duration (i.e., 
whether the fault was developing slowly or rapidly) 
of the identified fault. The networks were trained 
using data from a jet engine simulation program 
(Refs. 4 and 5). The system exhibited the ability to 
render accurate and reliable diagnoses at high noise 
levels, and was able to generalize to correctly iden¬ 
tify fault patterns which were not included in the 
training sets. 

Neural network models were also employed in 
the development of diagnostic techniques applied to 
the rocket engine domain. The requirement for real¬ 
time response in the rocket engine domain is even 
greater than that of the jet engine domain. The 
rocket engine domain is characterized by the pos¬ 
sibility of very rapid development of fault condi- 



Figure 3. Autoassoclative Neural Network Topology 


tions, with potentially catastrophic results if 
remedial action is not taken quickly. 

The neural network models used in the 
prototype engine diagnostic systems were conven¬ 
tional feedforward networks trained using well- 
known back-propagation algorithms (Ref. 6). 
Training was accomplished by clamping the input 
and output nodes to patterns representing Bensor 
data and fault conditions, respectively, and then 
applying the back-propagation algorithm until the 
network weights effected a proper mapping between 
input and output patterns. Although these 
methodologies were effective in diagnosing faults 
when given samples of data corresponding to a fault 
for which the networks were trained, several 
deficiencies of the conventional feedforward model 
were noted: 

1. Networks can be trained to associate a set of 
patterns with a set of fault conditions quite readily. 
However, when presented with an input pattern 
qualitatively different from those included in the 
training examples (e.g., a pattern representing a 
previously unobserved fault condition), the net¬ 
works can produce spurious categorizations in an 
unpredictable manner and produce false positive 
identifications. There is no method to train a con¬ 
ventional feedforward network to, for instance, ac¬ 
tivate an output node if the input pattern is not 
similar to those patterns for which the network is 
trained. Instead, the diagnostic systems had a ten¬ 
dency to produce a positive categorization cor¬ 
responding to the fault in the training set which, in 
some sense, was "closest" to the input. Referring 
back to Fig. 2, the output node marked "Unknown 
cannot be implemented with conventional neural 
network topologies. 
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2. Conventional feedforward networks are 
limited in the number and type of training examples 
which can be used. Some of the networks in the jet 
engine diagnostic system required long training 
times, which were a function not only of the number 
of training examples required, but also how similar 
training examples representing different faults 
were to each other. 

These deficiencies can adversely effect the ap¬ 
plication of neural networks to the classification of 
fault patterns* The present work attempts to formu¬ 
late a neural network topology and training regimen 
which will successfully extend the capabilities of 
feedforward network models. 


The Autoassodatlve Back-Propagation Model 


The network topology employed in the present 
work is autoas soaative , with at least one layer of 
semi-linear hidden nodes, as depicted in Fig. 3. Node 


activations are continuous variables. The autoas- 
sociative network contains the same number of 


input and output nodes and is trained to map or 
associate input patterns with identical patterns 
clamped at the output nodes. As a result, training is 
unsupervised; it is left to the weights between the 
input, hidden, and output layers of the network to 
organize in a fashion such that the mapping is 
performed correctly. In the present effort, three 
layers of nodes are used. The input and output 
layers contain 40 nodes; the hidden layer contains 
three nodes. Only the hidden layer nodes are semi- 
linear, the input and output nodes are linear. There¬ 
fore, the activations of nodes in the input and output 
layer are not restricted to specified ranges. 


Training The Network 

In the present effort, training the autoassocia- 
tive model is accomplished in two distinct steps. 
Suppose, for instance, it is desired to train a network 
to classify a large number of temporal patterns from 
engine sensors. In the first step a small repre¬ 
sentative subset of the total set of curves is selected. 
This subset of curves is used to obtain a system of 
weights which will perform the autoassoaative 
mapping, using a back-propagation algorithm. 

In th^ second step, the entire set of available 
data is input sequentially to the network. Each 
pattern will invoke a specific hidden node response. 
When patterns are related, i.e. exhibit similar 
qualitative features, the hidden node responses will 


*A variation of the autoaaaociative network topology ha* been 
previously applied to 2-D image analysis and image compression, 
as reported in Ref. 7. 


map out a surface. If N hidden nodes are used, the 
hidden node activations will define a surface in 
N-dimensional space. For example. Fig. 4 illustrates 
surfaces generated by training an autoassodative 
network with three hidden nodes on exhaust gas 
temperature data from a jet engine. The larger 
surface is generated from hidden node activations 
which correspond to a fuel interruption; the smaller 
surface corresponds to a bearing failure. Both sur¬ 
faces represent fault conditions over a range of 
severities and durations. A single point in three 
dimensional space represents the hidden node ac¬ 
tivation resulting from normal (non-transient) 
operation. 

The first step of training, using the back- 
propagation algorithm, requires most of the total 
computational effort in the two-step training 
process. However, the amount of computational ef¬ 
fort required using the subset of ex a m ples is much 
less than would be required if all available patterns 
were used in the training set. For instance, the first 
step of training for the case in Fig. 4 used a repre¬ 
sentative set of 24 patterns out of an available 450. 
In the present effort, failure scenarios of moderate 
severity were used. It is known that some subsets of 
data produce better classification performance than 
others; methods of selecting an optimal subset are 
presently under investigation. 

In the second step of training, the only computa¬ 
tional burden is to obtain hidden node activations 
for every available curve. The second step of training 
for the case in Fig. 4 utilized all 450 examples, but 
required a relatively insignificant amount of com¬ 
putational effort compared to the first step. The end 
result is a network that, in effect, is trained on all 
460 patterns, with a total training time far less than 
would be required to conventionally train a network 
on all 460 examples using back-propagation alone. 

It should be stressed that the training process 
does not in any way require that only three hidden 
nodes be used, but is fully extendible to higher 
numbers of hidden nodes. If higher numbers of 
hidden nodes are used, higher-dimensional surfaces 
(hypersurfaces) will be generated. 

Classification Of Input Patterns 

Once a network is trained and the hypersurfaces 
generated, it may be used for classification of input 
patterns. An input pattern will be mapped into a 
point in multi-dimensional space. If the point lies on 
or near a surface associated with a specific fault, the 
pattern is classified as an example of an input cor¬ 
responding to that fault. Fig. 4 depicts the point in 
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Figure 4. EGT Surfaces 


three-dimensional space representing a fuel inter¬ 
ruption scenario. If the point lies in a region off the 
surface (beyond a preset tolerance), the input is 
classified as corresponding to a type for which the 
network was not trained, as illustrated in Pig. 4. In 
this sense a network can determine if it has been 
properly trained to recognize an input; if the input 
is not recognized as belonging to a class correspond¬ 
ing to the training set, it can be identified as an 
instance of another class and processed further, if 
desired, by other (e.g., knowledge-based) methods. 

The position of a point on the surface also may 
be used to further quantify aspects of an identified 
input pattern. In the case of the jet engine diagnostic 
system, each fault pattern is associated with a 
severity and duration. As a result, each point on a 
hypersurface is associated not only with a particular 
fault, but is also associated with a unique seventy 
and duration. Therefore, an input pattern which 
results in hidden node activations which fall on a 
surface can be associated with a fault severity and 


duration. In the example in Fig. 4, the hidden node 
activation corresponding to a fuel interruption 
scenario falls on a point on the fuel interruption 
hypersurface corresponding to a severity and dura¬ 
tion of 75% and 1.5 seconds, respectively. 

Diagnostic Process 

Diagnosis proceeds in two steps. First, it is 
determined whether a fault condition is occurring, 
and, if so, which fault scenario is most likely. This 
step is fault identification. Secondly, the seventy 
and duration of the fault is determined. 

iHantification qf F aillt Scenarios. 

From the standpoint of fault identification, it is 
desirable if hypersurfaces corresponding to dif¬ 
ferent faults are widely separated. The hypersur¬ 
faces depicted in Fig. 4 are widely separated and 
therefore can be used reliably to distinguish be¬ 
tween fuel interruptions and bearing failures. How¬ 
ever this fortunate circumstance will not always 
occur. The position of the hypersurfaces is related to 
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Figure 5. N1 Surfaces 


the q uali tative and quantitative features of a fault 
scenario. Different fault scenarios may produce 
ni-milnr responses in a given sensor. For instance, 
Fig. 6 depicts N1 hypersurfeces for fuel interruption 
and bearing failure scenarios. The proximity of the 
two surfaces indicates that, in some cases, the sen¬ 
sor responses to the two faults will be similar. In 
addition, since the two surfaces intersect, some com¬ 
binations of fault severity and duration will invoke 
identical sensor responses. Hence, the usefulness of 
the N1 seftsor in discriminating between bearing 
failures and low severity fuel interruptions is 
limited. However, the use of other sensors will al¬ 
most always remove any ambiguities introduced by 
the N1 sensor (see Fault Arbitration, below). 
Moreover, the response of N1 to higher severity fuel 
interruptions may be used successfully, since t e 
corresponding region of the hypersurface is isolated 
enough to avoid ambiguities. 


Although failure scenarios are represented as 
hypersurfaces, normal (noise-free) engine behavior 
is represented as a point in hyperspace. The 
desirable spatial relationships between the normal 
point and the hypersurfeces are identical to those 
between hypersurfeces, i.e., it is desirable for the 
point representing normal behavior to be separated 
from the hypersurfeces. In Fig. 4, it is dear that this 
desirable condition is not met to a sufficient degree. 
The point representing normal steady-state be¬ 
havior is quite close to the hypersurface repre¬ 
senting bearing failures. This is because bearing 
failures do not cause severe deviations from the 
steady-state values for EGT. As a result, EGT is not 
a reliable sensor to use for discriminating between 
normal and bearing failure scenarios, especially in 
high noise environments. 

The situation is much better in the case of Nl, 
however, as can be seen in Fig. 5. In this case the 
point representing normal behavior is almost equi- 
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distant from both the fuel interruption and bearing 
failure hypersurfaces, indicating that N1 is a reli¬ 
able means of discriminating between normal be¬ 
havior and fault conditions. Therefore, EGT is best 
used to discriminate between bearing failure an 
fuel interruption, and N1 is best used to dis¬ 
criminate between normal behavior and either fault 
condition. Knowledge of the diagnostic strengths 
and limitations of each sensor can be used very 
effectively to resolve situations where, for instance, 
a particular sensor detects an incoming pattern that 
has characteristics of more than one faul t scenario. 
This process is performed by the arbitration system 
(see Fig. 2). 


Fault Arbitration. 

Fault arbitration is essentially a process 
whereby final fault identifications are made by 
weighting some sensors more than others, based on 
the ambiguous fault identifications initially ob¬ 
tained. For example, consider a low-severity bearing 
failure. The hidden node activations for EGT and Nl 
will be near the point representing normal engine 
behavior. In the case of EGT, the hidden node ac¬ 
tivation will also be close to the hypersurface cor¬ 
responding to a bearingfailure. As a result, EGT will 
identify the incoming signal as being characteristic 
of either a bearing failure scenario or normal be¬ 
havior. The N1 sensor will unambiguously identify 
the fault as a bearing failure. Overall, the diagnostic 
system at this point has identified either a bearing 
failure or normal behavior, a "tie" exists. However, 
it is known that Nl is quite good at breaking ties 
between bearingfailure and normal behavior; there¬ 
fore, the results of N1 are weighted more heavily 
than those of EGT in this case. The final result is an 
unambiguous identification of a bearing failure. 


The weighting would be different in the case of 
a tie between a fuel interruption and bearing failure 
scenario. In this instance, it is known that EGT is 
superior in discriminating between fuel interrup¬ 
tions and bearing failures; as a result theweight 
associated with the results from EGT would be lngh 
and the weight associated with N1 would be low. The 
final result would be an unambiguous identification 
based ijn the results of EGT. 


In the present effort, four sensors in the jet 
engine are monitored. The weights assigned to each 
sensor in the event of ambiguous initial diagnoses 
are presented in Table 1^. As stated earlier, nor- 

1 In the present effort, binary weightings are used; however , 
fractional weightings are entirely plausible and may be necessary 
for optimal performance. Optimal weighting assignments are 
currently being investigated. 


mal/bearing-failure ties are broken by Nl. No sen¬ 
sors are suppressed in the event of a tie between fuel 
interruption and normal behavior scenarios, since 
all sensors can adequately discriminate between 
those scenarios. Nl, as indicated previously, cannot 
reliably distinguish between fuel interruptions and 
bearing failures; in the event of a tie of that nature, 
all sensors but Nl are used for a final diagnosis. 


Table I 


Sensor 

Dual Fault Indication 

Nl 

or 

N2 

EGT 

1.0 

1.0 

1.0 

1.0 

Normal & Fuel 

1.0 

0.0 

0.0 

0.0 

Normal & Bearing 

0.0 

1.0 

1.0 

1.0 

Fuel & Beanne 


Determinatio n Qt Optimal Radius Of InflUgflCO . 

As stated previously, fault identification is ac¬ 
complished by determining the proximity of a hid¬ 
den node activation produced by an input pattern to 
a hypersurface. In the present system, a preset 
tolerance is used to determine how close a hidden 
node activation should be to a hypersurface to be 
considered "on" the hypersurface. In effect, a hidden 
node activation has an associated "radius of in¬ 
fluence"; if the distance to a hypersurface is less 
than this radius, the hidden node activation is con¬ 
sidered to be dose enough to be identified with the 
attributes of the hypersurface. 


The ideal radius of influence for the present 
lystem was determined by experiment. A large 
lumber of fault scenarios were presented to the 
system, with an associated noise level of 5%. The 
radius of influence was varied from 0.0 to 2 . 0 , and 

the resultant percentages of correct fault identifica¬ 
tions were obtained, as depicted in Fig. 6. For small 
radii (e.g., less than 0.1), all input fault patterns 
were dasaified as unknown fault classes. This is 
expected, as it is unlikely that any noisy input 
pattern will invoke a hidden node activation close 
enough to a hypersurface to be dassified proper y. 
As the radius of influence increases, the percentage 
of correct fault identifications increases Correct 
identifications of fuel interruptions and bearing 
failures increase to 100% and remain at that level 
for radius values of greater than - 0.6. However, the 
percent of correct classifications for nolB y 
input patterns peaks at a radius of about 0.9. This 
is because, as the radius increases beyond 0.9, some 
of the normal input patterns are incorrectly das¬ 
sified as faults. This peak indicates dearly that an 
optimal radius of influence exists, and that the 
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Figure 6. Diagnostic Accuracy vs Radius Of Influence 

diagnostic system can in effect be "tuned" to maxi¬ 
mize performance for a given noise level. 

Results 

Results of diagnoses for a radius of influence of 
0.9 are presented in Fig. 7 for noise levels from 0% 
to 20%. A large number of normal and failure (fuel 
interruption, bearing failure, and the unknown 
bleed pipe rupture) scenarios were processed by the 
diagnostic system. The illustration depicts the per¬ 
centage of correctly identified faults before and after 
arbitration, as a function of noise level. The effects 
of the arbitration process on the overall fault iden¬ 
tification process is pronounced, particularly in the 
identification of normal input patterns. Initially, all 
normal scenarios are identified as either normal or 
bearing failure scenarios. After the arbitration 
process, normal scenarios are correctly identified up 
to fairly high noise levels. Bearing failure scenarios 
initially confused with normal scenarios are, after 
arbitration, correctly identified at a level of 100% up 
to nearly 5% noise levels. As expected, fuel interrup¬ 
tion identifications are only very slightly affected by 
noise levels up to 10%. At extremely high noise 
levels, the hidden node activations of normal and 
fault scenarios are usually far from any hypersur- 
face, and are therefore classified as unknown pat¬ 
terns. Identification of a fault condition for which 



_BEFORE ARBITRATION 

_ AFTER ARBITRATION 



Figure 7. Diagnostic Accuracy vs Noise Level 


the network was not trained (bleed pipe rupture) is 
constant at nearly 100% for all noise levels This is 
also expected, since the addition of noise will usual y 
cause the associated hidden node activations to fall 
even farther away from the hypersurfaces when 
compared to the no-noise case. 


It is interesting to note that, before arbitration, 
the bearing failure and normal scenarios are actual¬ 
ly correctly identified more often under noisy condi¬ 
tions. This is because random noise fluctuations 
may produce hidden node activations which are less 
ambiguous than those produced by the noise-free 
conditions. For instance, a noisy normal behavioral 
pattern may produce a hidden node activation which 
is farther away from the bearing failure hypersur¬ 
face than that produced in a noise-free environment. 
In this case, there will be no confusion between 
normal behavior and a bearing failure scenario. A 
similar process occurs with true bearing failure 
scenarios; noise may move a hidden node activation 
away from the point designating normal behavior, 
thereby removing any ambiguity. However, arbitra¬ 
tion significantly improves results both in noisy and 
noise-free environments. 
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(b) Duration Error 


Figure 8. Error Distribution For Fuel Interruption Severity 
(a) and Duration (b) 

Determinati o n of Fault Severity find Duration . 

Once a fault has been identified, the hypersur¬ 
faces may be used to estimate the severity and 
duration of the fault. Each point on a hypersurface 
is associated with a unique fault severity and dura¬ 
tion. During the fault identification process, a hid¬ 
den node activation corresponding to the input 
pattern will be generated. To estimate fault severity 
and duration, the closest point on the hypersurface 
to the point representing the hidden node activation 
may be'found; the closest point will correspond to a 
unique fault severity and duration. 

Fault identification is mostly affected by the 
separation of the hypersurfaces from each other and 
from the point representing normal steady-state 
behavior. Estimates of fault severity and duration 
are affected primarily by the size of the hypersurface 
representing the identified fault. This is because, in 


Figure 9. Error Distribution For Bearing Failure Severity 
(a) and Duration (b) 


noisy environments, the perturbations of the hidden 
node activations will be relatively more severe for 
small hypersurfaces than large ones. FigB. 8 a) and 
b) depict histogramB of fuel interruption seventies 
(expressed as decimals) and durations, respectively, 
for a 6% noise level. The standard deviation (c) for 
severities is about 0.03, which means that ap¬ 
proximately 99% of estimated severities will be 
within 0.06 (2a) of the true values. The correspond¬ 
ing standard deviation for fuel interruption dura¬ 
tions is 0.3 seconds. As can be seen in Figs. 4 and 5, 
the hypersurfaces for fuel interruption scenarios are 
quite large. As a result, good estimates of fuel inter¬ 
ruption severity and duration may be obtained. 
However, the hypersurfaces associated with bearing 
failures are relatively small. This has a detrimental 
impact on the ability to accurately estimate the 
severity and duration of bearing failures. As can be 
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seen in Figs. 9 a) and b), the 2 a limits for bearing 
failure severity and duration are 0.18 and 0.9 
seconds, respectively. 

Conclusions 

A diagnostic system based on neural networks 
has been developed and applied to the real-time 
diagnosis of jet and rocket engines. The approach 
used in the present work addresses deficiencies in¬ 
herent in many feedforward neural network models. 
The present system greatly reduces the number of 
false positive diagnoses that can be produced when 
applying standard feedforward networks to pattern 
classification. In addition, the present approach 
greatly reduces the number of networks necessary 
to 1) identify the existence of a fault condition, and 
2) estimate the duration and severity of the iden¬ 
tified fault. The present method also incorporates a 
two-stage training approach which greatly reduces 
the overall training time. 
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Abstract 

The Space Shuttle utilizes several rocket engine 
systems, all of which must function with a high 
degree of reliability for successful mission comple¬ 
tion. The Space Shuttle Main Engine (SSME) is by 
far the most complex of the rocket engine systems. 
In earlier spacecraft, rocket systems (and, in fact, 
the entire spacecraft) were designed for use on only 
a single mission and were discarded after use. In a 
major departure from earlier practices, almost all 
systems on the Space Shuttle, including the rocket 
systems, are designed to be reusable; only the exter¬ 
nal tank is discarded during each mission. 

The reusability of spacecraft systems introduces 
many problems related to testing, reliability, and 
logistics. Components must be assembled from 
parts inventories in a manner which will most effec¬ 
tively utilize the available parts. Assembly must be 
scheduled to efficiently utilize available assembly 
benches while still maintaining flight schedules. 
Assembled components must be assigned to as many 
contiguous flights as possible, to minimize com¬ 
ponent changes. Each component must undergo a 
rigorous testing program prior to flight. In addition, 
testing and assembly of flight engines and com¬ 
ponents must be done in conjunction with the as¬ 
sembly and testing of developmental engines and 
components. 

The development, testing, manufacture, and 
flight assignments of the engine fleet involves the 
satisfaction of many logistical and operational re¬ 
quirements, subject to many constraints. The pur¬ 
pose of the SSME Component Assembly and Life 
Management Expert System (CALMES) is to assist 
the engine assembly and scheduling process, and to 
insure that these activities utilize available resour¬ 
ces as efficiently as possible. 

PROBLEM DESCRIPTION 

Fig. 1 illustrates the general problem CALMES 
is designed to address. Typically, a priori knowledge 
consists only of parts inventories and a proposed 
flight schedule. Using this input, CALMES per¬ 


forms four basic functions: 1) determining how com¬ 
ponents may be optimally assembled from the parts 
inventory, 2) assigning components to flights, 3) 
scheduling component testing, 4) scheduling com¬ 
ponent assembly. The first function deals exclusive¬ 
ly with component assembly, while the remaining 
functions are considered to be scheduling opera¬ 
tions. The order in which these functions are per¬ 
formed is indicated by the direction of the arrows in 
Pig l The four major functions are explained in 
more detail below: 

Part-Comoon ftn* Assignment 
For the purposes of the present discussion, the 
fflMF. is considered to be composed of 1) major 
components, and 2) minor components (or parts). 
Major components include the four turbopumps, i.e., 
the high pressure fuel and oxidizer turbopumps 
(HPFTP and HPOTP, respectively) and the low 
pressure fuel and oxidizer turbopumps (LPFTP and 
LPOTP, respectively). Minor components or parts 
comprise the constituent parts of the major com¬ 
ponents, and cannot be subdivided further. Minor 
components or parts indude turbine blades, im¬ 
pellers, impeller shafts, housings, and turbine discs. 
Both major and minor components are rated accord¬ 
ing to remaining allowable life, and remaining al- 


ASSEMBLY 


PARTS 
INVENTORY^ 


PART-COMPONENT 

ASSIGNMENT 
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COMPONENT-ASSEMBLY 
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Figure 1. CALMES Input And Main Functions 
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lowable starts. If. for example, the HPOTP has 800 
seconds of remaining allowable life, and two remain¬ 
ing allowable starts, it may be used once for 800 
seconds, or twice for a total of 800 seconds. 

It is relatively easy to calculate allowable 
lifetimes and number of starts for minor com¬ 
ponents. Each minor component is assigned an al¬ 
lowable life and number of starts at the time of 
manufacture; these estimates are based on. en¬ 
gineering design considerations and extensive test¬ 
ing. During the life of a part, the length of time the 
part is used and the associated number of starts are 
simply subtracted from the initial values each time 
the part is used. 

The allowable life and number of starts of a 
major component, however, are determined by the 
smallest allowable lifetime or number of starts of its 
parts. Therefore, it is critical that major components 
be assembled from parts with sufficient remaining 
life and number of starts to accommodate an¬ 
ticipated test and flight requirements. Generally, 
there will be many ways in which available parts 
may be combined into major components. A single 
turbopump consisting of five parts, where each part 
comes from an inventory of 60 parts, may be as¬ 
sembled in almost a billion different ways. However, 
many of the possible component configurations will 
not satisfy the required number of starts or allow¬ 
able lifetimes, based on anticipated flight schedules 
and durations. It is the task of the present system 
to assemble a set of major components which effi¬ 
ciently utilizes the available resources from a large 
parts inventory. 

< ^rpponent-Fl |qht Assignment. 

Each flight of the shuttle requires three engines 
which have remaining allowable lifetimes greater 
than the longest anticipated shuttle flight (i.e., 770 
seconds, which would be required in an abort situa¬ 
tion). Of course, each engine, and therefore each 
f ngini» component, must have at least one allowable 
start. As in the case of major components, the 
remaining life and number of starts of engines are 
determined by the smallest remaining life or num¬ 
ber of starts of its components. Generally, it is 
desirable to assign an engine or component to 
several contiguous flights to minimize engine chan¬ 
ges. An engine spanning more than one flight re¬ 
quires longer allowable lifetimes and multiple 
allowable starts. The task of CAIMES is to assign 
engines and major components to flights in an effi¬ 
cient manner. Generally, this means that 1) engines 
and components will be assigned to as many con¬ 
tiguous flights as possible, and 2) the remaining life 


of an engine or component after its use on flights will 
be smaller than that of any flight remaining on the 
schedule. 


r^ m ponpnt-Test Stand Scheduling. 

There are presently several types of engine 
ground tests, including 1) development, 2) certifica¬ 
tion, 3) flight certification extension, 4) acceptance, 
and 5) component "green runs". Development tests 
are performed on hardware which has not been 
certified forflight. Certification tests are performed 
on hardware which has matured to a point where 
use on a flight is desired. Flight certification exten¬ 
sion tests are performed on flight certified hardware 
to extend the certified hardware lifetime or number 
of allowable starts. Acceptance tests are performed 
on entire engines prior to installation on a shuttle 
orbiter. Green run tests are designed to identify any 
problems in a component in a controlled environ¬ 
ment before it is installed in an engine which will be 
used on a manned flight Usually, green runs consist 
of three engine firings of 1.6, 250, and 520 seconds 
over the course of one week. Other types of tests may 
range from one to several weeks, and may involve a 
wide range of engine firing times and starts. 

CALMES is designed to assign test stands to 
components or engines to fulfill test requirements. 
The assignment process attempts to 1) schedule an 
available test stand within a specified time range, 
and 2) to preferentially assign a particular test 
stand usually used for particular types of tests. The 
former requirement is more important than the 
latter, and therefore takes precedence. 


fyynponent- A ssnmbtv Bench Schedu l ing . 

The final function of CALMES is to assign as¬ 
sembly benches to components which are required 
for ground tests. Presently, there are six available 
assembly benches. Components may take several 
weeks to assemble, depending on the fype 
ponent and the condition of its parts. CALMES is 
designed to assign assembly benches to «tber 1) 
utilize as many assembly benches as possible, an 
hence minimize the overall time required for the 
assembly of all components, or 2) use as few assemb¬ 
ly benches as possible, but assembling the com- 
nnnents as late as possible. 


Simplifying A ssumptions. 

The approach used in the prototype has been 
designed to be able to accommodate the general 
problem outlined above. However, the following 
simplifying assumptions have been made forthe 
purpose of verifying the approach used in CALMES. 
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1) A single type of major component (high-pres¬ 
sure oxidizer turbopumps, or HPOTPs) is con¬ 
sidered. More major components will be considered 
in the operational program. 

2) HPOTPs are assumed to consist of five life- 
limiting parts, i.e. blades, impeller, housing, im¬ 
peller shaft, and turbine disc. Additional parts may 
be considered in the operational program. 

3) A single launch vehicle is assumed. The 
operational program will consider three launch 
vehicles. 

4) Only green run tests are scheduled. Engine 
certification and development tests, etc. are not 
presently incorporated into test schedules. 

prnhlpm Formulatio n for the Present SvStettL 

The problem subset being addressed by the ex¬ 
pert system prototype may be summarized by the 
following processes: 

1) Given an inventory of parts (e.g., turbine 
blades, etc.), assemble a virtual fleet of major com¬ 
ponents (e.g., turbopumps). A virtual fleet is a set of 
pumps that could be built if desired; however, only 
pumps actually required by the flight schedule will 
be assembled. 

2) Assign components to shuttle flights. This 
step determines which components from the virtual 
fleet are actually required. 

3) Schedule green run ground tests on com¬ 
ponents required for flight. 

4) Schedule assembly benches on components 
required for ground tests. 

Hard and Sa lt Constraints. 

All of the above steps must be performed subject 
to many constraints. For example, the pumps must 
be assembled from parts which will result in allow¬ 
able lifetimes sufficient to accommodate the desired 
flight schedule. Resources such as test stands and 
assembly benches must be scheduled to meet flight 
deadlines, testing requirements, and time restric¬ 
tions. Also, engines must be assigned to flight 
vehicles in a manner which minimizes engine chan¬ 
ges over the proposed flight schedule. 

Constraints on the assembly and scheduling 
processes may be categorized as being either hard 
or soft. Hard constraints are inviolable; an example 
is a requirement of three months between the com¬ 
pletion of a test on a pump and the use of the pump 
in a flight engine. Soft constraints are preferences 


which may be relaxed, if necessary, to accommodate 
hard constraints. An example of a soft constraint is 
the preferred use of a particular test stand for cer¬ 
tain types of tests, or the use of as few assembly 
benches as possible during assembly of a major 
component. The approach used in the present work 
allows the incorporation of both hard and 3oft con¬ 
straints into the assembly and scheduling process. 
In addition, soft constraints may be given relative 
degrees of importance. 

General Approach 

All of the major functions of CALMES have one 
aspect in common: they all require assignments of 
one object or entity to another. For instance, the 
part-component assignment problem assigns parts 
from an inventory to instances of components. Com¬ 
ponent-test stand assignment requires the assign¬ 
ment of components to available test slots on 
particular test stands. The four major functions of 
CALMES are each cast as instances of the formal¬ 
ized assignment problem, a familiar problem in 
operations research. 

The assignment problem is illustrated in Figs. 
2 and 3. Fig. 2 a) depicts a bipartite graph, i.e., a 
graph which consists of two distinct sets of nodes. 
The two sets of nodes are connected by weighted 
links. The links between node B and the nodes of the 
other side of the graphs are illustrated in Fig. 2 a). 
Each link is associated with a weight, which in the 
present implementation is a positive real number. 
T.inlrg between the other nodes are omitted for 
clarity. 

The system of weights can be conveniently rep¬ 
resented in matrix form, as depicted in Fig. 2 b). An 
optimal match is obtained when 1) each node on one 
side of the graph is matched with a single node on 
the other side of the graph, and 2) the sum of the 
weights (the cost) is a minimum. The circled weights 
in Fig. 2 b) illustrate the optimal matching for this 
graph. Note that each link assigned to the optimal 
match occupies a unique row and column of the 
weight matrix. 

optimal Matching. 

The process of finding an optimal match is il¬ 
lustrated in Fig. 3. The first step of the matching 
process is to select, as an initial matching, all pair¬ 
ings represented by the smallest weights in the 
columns of the matrix. This process results in the 
pairings A2, Cl, and D6, as depicted in Fig. 3 a). 
Nodes 3 and 4 cannot be matched to nodes on the 
other side of the graph, because their minimum 
weighted links conflict with matchings A2 and Cl, 
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a) BIPARTITE GRAPH 


1 2 3 4 5 



Figure 2. Optimal Matching Of Bipartite Graph 

respectively. If nodes 3 and 4 are matched arbitrari¬ 
ly with nodes B or E, the lowest sum of weights 
possible is 20. 

At this point, three nodes on both sides of the 
graph have been temporarily committed (i.e., they 
have been paired with nodes on the other side of the 
graph). The matching is now improved by introduc¬ 
ing the concept of augmenting paths. An augment¬ 
ing path through the graph consists of alternating 
uncommitted and committed segments, where the 
first and last segments are uncommitted. Fig. 4 
shows an augmenting path in Fig. 4 a). At first, there 
is only one committed segment (i.e., A2). If the 
committed segments are made uncommitted, and 
vice versa, the path will have two committed seg¬ 
ments and only one uncommitted segment, as in Fig. 
4 b). An augmenting path will always increase the 
number of pairings by one. In effect, augmenting 


paths "undo" initial, inferior pairings and replace 
them with pairings which are closer to the optimal 
match. After a single augmentation, depicted in Fig. 
3 b), the lowest possible cost has been reduced from 
20 to 16. After two augmentations, illustrated in Fig. 
3 c), the cost has been reduced further to 15. At this 
point, all possible pairings have been made, and the 
optimal match has been attained. 

The choice of augmenting paths is obviously of 
primary importance. Many methods have been 
studied; CALMES uses a modified simplex method 
to search for augmenting paths and to obtain global¬ 
ly optimal matchings. Discussion of the matching 
algorithm, sometimes referred to as the Hungarian 
method, is beyond the scope of this paper, however, 
the basic approach is fully documented in Refs. 1 
and 2. 

Results 

Assignment of Parts to Pumps 

The first function CALMES performs is the as¬ 
sembly of a fleet of pumps from an inventory of 
parts. In the present simplified implementation, 
only high pressure oxidizer turbopumps are con¬ 
sidered; in addition, only five different parts are 
assumed: turbine blades, impeller shaft, housing, 
turbine disc, and impeller. Any number of these 
parts can be in the inventory, each with associated 
allowable lifetimes and starts. The task of CALMES 
is to "build" as many pumps from the available 
inventory, murimiring either the allowable lifetime 
or the number of starts of the fully assembled fleet. 
The resultant set of pumps is "virtual", in that any 
of the pumps may possibly be actually built; how¬ 
ever, CALMES does not at this stage determine 
whether or not a particular pump should actually be 
assembled. 

The virtual assembly is initiated by optimally 
ma tching, for instance, the turbine blade inventory 
with the inventory of impellers. This is ac¬ 
complished by associati n g nodes on one side of the 
bipartite graph with instances of turbine blades in 
the inventory, and associating nodes on the other of 
the graph with instances of impellers. A matching 
between the two sides of the bipartite graph will in 
effect generate a set of partially assembled virtual 
pumps. 

The determination of the weights between the 
nodes is accomplished by applying appropriate 
heuristics. The resultant lifetime of any blade-im¬ 
peller matching will be the lower of the blade 
lifetime and the impeller lifetime. Therefore, it 
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a) INITIAL MATCHING (COST-20) 




C) TWO AUGMENTATIONS (COST -16) 


Pjgiifa 3. Development 01 Optimal Matching 

would not be productive to match a blade with a high 
lifetime with an impeller of a low lifetime, providing 
a blade with a lifetime close to that of the impeller 
is available in the inventory. Therefore, the weights 
between blades and impellers are made proportion¬ 
al to the absolute value of the difference of the blade 
and impeller lifetimes. As a result, a blade and 
impeller with equal lifetimes will most likely be 
matched; a match between blades and impellers 
with greatly differing lifetimes would be unlikely, 
but of course would not be impossible if such a 
matching resulted in a globally optimal matching. 

Once a virtual fleet of partially assembled 
pumps is obtained, the matching process is repeated 
between the partially assembled pumps and 
another parts inventory. In the present system, the 
fleet of partially assembled pumps will be matched 


o o—o o 

3 A 2 E 


a) BEFORE PATH AUGMENTATION 


o—o o—o 
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b) AFTER PATH AUGMENTATION 


Figure 4. Path Augmentation 

with turbine housings. Again, the weights between 
the partially assembled pumps and housings are 
made proportional to the difference in the respective 
allowable lifetimes. The second matching process 
results in another virtual fleet of partially as¬ 
sembled pumps; however, the pumps now have 
three parts: blades, impellers, and housings. 

Assignment of the partially assembled pumps 
with the two remaining parts is accomplished with 
two more matching steps. Weights are determined 
on the of allowable lifetimes of the partially 
assembled pump inventory and the inventories of 
turbine discs and impellers. The final result is a 
virt ual , completely assembled pump fleet. 

Abbreviated output from CALMES illustrating 
the part-component assignment process is 
reproduced in Pigs. 6-7. Pig. 5 shows the assignment 
of blade sets to individual pumps. The pumps are 
designated by alphabetic characters. Since both 
upper and lower cases are used, 62 pump instances 
can be designated. The first entry states that the 
second blade set (serial number 1000002, allowable 
life of 1723 seconds, 9 allowable starts) is assigned 
to pump T. Assignments for the remaining four 
parts are listed by CALMES in a similar manner, 
but are not reproduced here. 

At the bottom of Pig. 5 are blade sets which are 
not assigned to any pump. Parts with the highest 
lifetimes will generally not be assigned, because the 
only possible assignments would group them with 
other parts with lower lifetimes. Such an assign¬ 
ment would not increase the overall lifetime of the 
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PART-PUMP ASSIGNMENTS 
BLADES 


PART 

NAME 

S/N 

LIFE 

STARTS 

PUMP 

2 

BLADES 

1000002 

1723 

9 

T 

3 

BLADES 

1000003 

1020 

3 

c 

7 

BLADES 

1000007 

1951 

13 

M 

9 

BLADES 

1000009 

2003 

12 

I 

10 

BLADES 

1000010 

773 

7 

h 

81 

BLADES 

1000081 

509 

15 

P 

82 

BLADES 

1000082 

1975 

15 

K 

87 

BLADES 

1000087 

1791 

1 

S 

88 

BLADES 

1000088 

1549 

1 

u 

89 

BLADES 

1000089 

1929 

2 

N 


UNASSIGNED 

PARTS 




PART 

NAME 

S/N 

LIFE 

STARTS 

1 

BLADES 

1000001 

1317 

10 

4 

BLADES 

1000004 

988 

13 

5 

BLADES 

1000005 

2187 

3 

6 

BLADES 

1000006 

2298 

5 

8 

blades 

1000008 

2082 

11 

92 

BLADES 

1000092 

1494 

1 

93 

BLADES 

1000093 

858 

11 

94 

BLADES 

1000094 

2409 

10 

95 

BLADES 

1000095 

350 

8 


Figure 5. Assignment Of Parts (Blade Sets) To Pumps 

virtual fleet, and would thereby be a suboptimal 
assignment. Parts with low remaining lifetimes will 
not be assigned, because their assignment would 
lower the overall fleet lifetime. 

Parts are assigned to pumps in a manner which 
ensures an optimal pump fleet, in the sense that no 
rearrangement of parts will result in a mix of pumps 
with greater allowable lifetimes. The assignment 
algorithm allows pump fleets to be optimized on the 
basis of life, starts, or an arbitrary combination of 
desirable factors. 

Fig. 6 is a table summarizing the virtual pump 
fleet resulting from the parts assignments. Forty- 
seven pumps can be assembled from the parts in¬ 
ventory, with allowable lifetimes from 2324 seconds 
to 79 seconds. As stated above, this is the best 
possible fleet of pumps which can be obtained from 
the given parts inventory. Fig. 7 is a table which 
details which parts (by serial number) are assigned 
to each pump. 

As stated earlier, the results of Fig. 6-7 indicate 
a virtual fleet of pumps which can, if required, be 
assembled from the available parts. Which pumps 
should be actually assembled will depend on the 
requirements of the flight schedule. 


VIRTUAL PUMP FLEET 


>TANC2 

PUMP 

U/N 

LIFE 

STARTS 

DESIGNATION 

1 

HP OTP 

2000059 

2324 

1 

A 

2 

HP OTP 

2000027 

2311 

4 

B 

3 

HP OTP 

2000021 

2259 

3 

C 

4 

HPOTP 

2000064 

2229 

2 

D 

5 

HPOTP 

2000091 

2203 

1 

E 

43 

HPOTP 

2000044 

401 

3 

q 

44 

HPOTP 

2000009 

398 

2 

r 

45 

HPOTP 

2000046 

389 

3 

s 

46 

HPOTP 

2000031 

345 

3 

t 

47 

HPOTP 

2000094 

79 

1 

u 


Figure 6. Optimal Pump Inventory 


PUMP-PART ASSIGNMENTS 


PUMP 

BLADES 

HOUSING 

SHAFT 

TURBINE 

DISC 

MAIN 

A 

1000080 

2000059 

3000018 

4000042 

5000014 

B 

1000046 
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Assignmen t of Pumps to Flights. 

Once a virtual pump fleet has been assembled, 
the pumps must be assigned to engines which are 
scheduled to fly on the shuttle. At this point, 
knowledge of the proposed flight schedule is neces¬ 
sary. Under normal conditions, each flight requires 
the engines to fire approximately 660 seconds. To 
illustrate the flexibility of the present approach, a 
fictitious flight schedule with a range of engine 
firing requirements from 600 to 770 seconds is 
depicted in Fig. 8 a). (In real life, each flight would 
be assumed to require 770 seconds of engine firing 
time, to account for the possibility of worst-case 
flight abort scenarios). 

Pumps may, if remaining lifetime and starts 
permit, be used on multiple flights. Obviously, it is 
advantageous for pumps to be used on multiple 
flights which are contiguous, to minimize pump 
changes. Therefore, use of a pump on multiple 
flights is a soft constraint, i.e., it is a preferred 
feature in the assignments of pumps to flights, but 
it is not absolutely necessary for an assignment 
scenario to be acceptable. 
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Figure 8. Flight Schedule and Assignments 


The assignment of pumps to flights is initiated 
by grouping contiguous flights. Although the groups 
may be arbitrarily large, in the present implemen¬ 
tation the groups are restricted in size to five con¬ 
tiguous flights. Each group is assigned a node on one 
side of a bipartite graph. Nodes on the other side of 
the bipartite graph are associated with the virtual 
pump fleet generated in the previous step. Weights 
are assigned to links between the pump fleet and 
flight groups in the following manner: 

1. If a pump has a greater allowable lifetime and 
allowable starts than that represented within a 
flight group, the weight is equal to the difference 
between 1) the pump lifetime and 2) the total length 
of the flights in the flight group. 

2. If a pump has a smaller allowable lifetime 
than a flight group, or fewer number of allowable 
starts than are represented in the flight group, the 
weight is set to a large number. Since a pump must 
have a greater lifetime and number of starts than 
that represented by a flight group, an optimal as¬ 
signment which includes a link with a large-num- 
bered weight would represent the violation of a hard 
constraint. 

Using these heuristics, a pump will tend to be 
matched to a flight group with a lifetime require¬ 
ment nearly equal to (but less than) the lifetime of 
the pump. If a matching occurs which includes a link 
with a large-numbered weight, the hard constraint 


(2) has been violated. In an operational version of 
CALMES, some soft constraint could be relaxed 
until the assignment no longer contained the large- 
numbered matching. 

Generally, a flight schedule will allow several 
different arrangements of multiple flight groups. 
For instance, a 10-flight schedule will allow seven 
different arrangements of five flights. One arrange¬ 
ment is two groups of five flights (flights 1 through 
5 and 6 through 10). The other arrangements in¬ 
volve only a single group of five contiguous flights. 
The resulting groups will be flights 1 through 5, 2 
through 6, 3 through 7, etc. Each of these seven 
arrangements requires a separate matching process 
with the pump fleet. CALMES will retain the match¬ 
ing with the lowest global cost. 

Once all possible five-flight groups are matched, 
the process is repeated with four-flight groups. If no 
five-flight groups were matched, the 10-flight 
schedule in the example above would allow 13 ar¬ 
rangements of four-flight groups. Each arrange¬ 
ment of four-flight groups is optimally matched with 
the pump inventory; the best optimal match is 
retained. The process is repeated with flight groups 
of three, two, and finally single flights. At this point 
all possible matchings will have been made. The 
original pump inventory will have been assigned to 
single or multiple flights in an optimal manner, i.e., 
in a manner which minimizes the overall differences 
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Figure 10. Assembly Schedule Timeline 


between the lifetimes of the assigned pumps and the 
lifetime requirements of the flight groups. 

Fig. 8 depicts a) the input flight schedule 
(provided by the user), and b) the assignment of 
pumps to each scheduled flight. The input flight 
schedule is not realistic, as it schedules a shuttle 
flight each week from week 55 to week 66. However, 
it is included as an example of a "worst case 
scenario, since it is most demanding of ground test 
and component assembly resources. 

The assignment of pumps to flights is depicted 
graphically in a timeline format spanning 70 weeks 
(Fig. 8 b). The assignment graph indicates that 
pump "F" is assigned to engine E-2 for flights at 
weeks 56, 57 and 58, while pump T” is assigned to 
engine E-l for flights at weeks 55 and 56. The 
pump-flight assignment algorithm allows pumps to 
be assigned to flights in a preferential manner. In 


the example, single engines are to be assigned to as 
man y flights as possible, with early flights taking 
precedence over later ones. Other preference may 
be included conveniently if they are expressible in 
terms of weight assignments. 

Pump Tes t Scheduling 

Once pump-flight assignments are made, pump 
test scheduling is performed. This function is cast 
as an optimal match between nodes representing 1) 
dates on which test stands are available, and 2) 
dates on which pumps are required for flight. A hard 
constraint requiring three months between a pump 
test and its use on a flight must be accommodated 
by the test schedule. In addition, two soft con¬ 
straints are imposed, i.e., it is preferred that tests 
occur on test stand B-l, no earlier than four months 
before the flight. 
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A a in the pump-flight assignments, node match¬ 
ings which would represent a violation of a hard 
constraint are weighted with a large number. As a 
result, links between 1) a given flight and 2) test 
stand dates later than three months before the flight 
will be given large weights. The soft constraints are 
imposed by assigning small weights to links be¬ 
tween nodes representing test stand B-l and those 
representing pumps which have been previously 
assigned to flights. Slightly larger weights are as¬ 
signed to links between the pump nodes and nodes 
representing test stand A-2. In this manner the 
preference of using test stand B-l over test stand 
A-2 is imposed. An optimal match of the resulting 
bipartite graph will result in 1) a (hopefully) total 
avoidance of any match between nodes which would 
violate a hard constraint, and 2) a preferential 
match between nodes representing time slots on test 
stand B-l and the pump nodes. If an optimal match 
cannot be attained by assigning pump tests to test 
stand B-l, then test stand A-2 will be assigned some 
pumps for testing. If an optimal match can be ob- 
tained only by violating a bard constraint, a soft 
constraint must be relaxed until the hard constraint 
is satisfied. 

The results of scheduling pump tests for the 
given flight schedule are illustrated by the timelines 
depicted in Fig. 9. The flight assignment timeline is 
repeated for comparison. Note that some tests are 
scheduled for test stand A-2 so that no test is per- 
formed earlier than four months before the flight. In 
this example, all hard constraints are satisfied. 

Pump Ass embly Scheduling, 

The final step, scheduling of pump assemblies, 
is depicted in Fig. 10. As in the previous step, the 
scheduling algorithm allows preferences to be ac¬ 
commodated if possible. In this example, it is desired 
to schedule assemblies as late as possible, using as 
few assembly benches as possible. (Alternatively, a 
user could prefer to fill up all benches if possible, to 
only use certain benches during certain periods, or 
to specify variable assembly times, etc.) An obvious 
hard constraint is that pumps must be assembled 
before they are tested. 

The two sides of the associated bipartite graph 
will represent 1) available assembly benches and 
dates, and 2) pump testing dates determined in the 
previous step. As in the previous step, links between 
node pairs which would violate a hard constraint if 


the nodes were matched are given large weights. 
Node pairings which represent preferences, or soft 
constraints, are weighted with small weights; the 
magnitudes of these weights depend on the relative 
desirabilities of the associated node pairings. In the 
example, all required pumps can be assembled in 
time for their green run tests. Four assembly 
benches are required for the pump assembly. 

Summary 

A prototype of the SSME Component Assembly 
and Life Management Expert System (CALMES) 
has been developed and implemented at 
CASP/UTSI. The four major functions of CALMES, 
i.e., 1) the assignment of parts from an inventory to 
jnnta nrfts of high pressure oxidizer turbopumps, 2) 
the assignment of pumps to flights, 3) the schedul¬ 
ing of pump tests, and 4) scheduling of pump as¬ 
sembly benches are conceptualized as general 
optimal assignment problems. A graph-theoretical 
optimal matching algorithm, based on a modified 
simplex method, is applied to the major functions 
required by the SSME component assembly and 
scheduling processes. Heuristics pertaining to ideal 
and preferred scheduling and assembly constraints 
are incorporated into weighted links in the graphi¬ 
cal representations. The prototype successfully 
demonstrates the accommodation of hard and soft 
constraints in assembly and scheduling processes. 
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abstract 

An expert system, called LEADER, has been de¬ 
signed and implemented for automatic learning, de¬ 
tection, identification, verification and correction of 
anomalous propulsion system operations in real time. 
LEADER employs a set of sensors to monitor engine 
component performance, and to detect, identify and 
validate abnormalities with respect to varying engine 
dynamics and behavior. Two diagnostic approaches 
are adopted in the architecture of LEADER. In the 
first approach fault diagnosis is performed through 
learning and identifying engine behavior patterns. 
LEADER, utilizing this approach, generates few hy¬ 
potheses about the possible abnormalities. These hy¬ 
potheses are then validated based on the SSME de¬ 
sign and functional knowledge. The second approach 
directs the processing of engine sensory data and per¬ 
forms reasoning based on the SSME design, func¬ 
tional knowledge, and the deep-level knowledge, i.e., 
the first principles (physics and mechanics) of SSME 
subsystems and components. This paper describes 
LEADER’S architecture which integrates a design- 
based reasoning approach with neural network-based 
fault pattern matching techniques. The fault diag¬ 
nosis results obtained through the analyses of SSME 
ground test data are presented and discussed. 


INTRODUCTION 

Most of the present day knowledge based fault di¬ 
agnostic systems are not capable of performing fault 
diagnosis in real time. Many of them use a large 
number of production rules, that make the system 
response slow. Also, the interactions among the rules 
make addition and deletion of rules difficult. Further¬ 
more, a number of these diagnostic systems perform 
shallow reasoning that fails to diagnose faults vvhic 
require deep level knowledge. A number of these sys¬ 
tems also suffer from a major drawback, namely, the 
user is expected to understand the fault symptoms 
because of the lack of facilities to acquire knowledge 


automatically from the sensor data. Since these diag¬ 
nostic systems are generally unsuitable for real-time 
diagnostic applications, alternative approaches for de¬ 
veloping real-time diagnostic systems are needed. 

In recent years, a number of researchers have pro¬ 
posed diagnostic approaches for complex mechanical 
systems, some of them are based on causal reason¬ 
ing [1], qualitative reasoning with mechanism models 
[2*3], sensor-based reasoning [4], and pattern match¬ 
ing using neural networks [5]. These approaches 
have been applied to various domains such as digital- 
circuits, jet engine and other engineering domains. 
The expert svstem LEADER discussed in this paper is 
based on two’ different but complementing approaches 
that show a potential for real-time diagnosis oi taults 
occuring in a complex mechanical system such is a 
rocket engine. The two approaches are based on 
engine bahavior and engine design analyses, respec¬ 
tively. LEADER has been applied for jet engine and 
rocket engine fault diagnosis. The approaches and the 
techniques used in developing LEADER, along with 
the results of SSME fault diagnosis obtained through 
the analysis of ground test data are discussed in this 

paper. 

k£AD£R ARCHITECTURE -V. A PPROACH. 

The architecture of LEADER is divided into two 
diagnostic subsystems, and an integration subsystem. 
One of the diagnostic subsystem performs fault diag- 
nosis based on the analysis of SSME behavior. whe»£ 
eas the other diagnostic subsystem analyzes bo: 1L. 
design and functions to diagnose SSME faults. Each 
of the diagnostic subsystems performs diagnosis in¬ 
dependently of the other by sharing the Pressor 
time, or whenever possible by running in parallel. The 
two diagnostic subsystems allow LEADER to perform 
SSME fault diagnosis using two different but com¬ 
plementing approaches - one based on utilizing an 
patterns obtained from SSME behavior corresponding 
to the faults, and the other based on applying deep- 
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level reasoning to SSME design. Whenever the results 
from the two approaches concur, no further analysis 
is needed. Otherwise, suitable analysis followed by 
integration of the results from the two approaches is 
performed by the integration subsystem of LEADER. 
This is expected to result in an overall improvement 
in LEADER’S performance in terms of its reliability, 
speed, effeciency, and accuracy of results. As shown 
in figure 1, LEADER’S architecture is divided into the 
following three subsystems: 

A. Engine Behavior Analysis Based Diagnostic Sub- 
sytem (BABS). 

B. Engine Design Analysis Based Diagnostic Sub- 
system (DABS). 

C. Integration Subsystem (IS). 

The three subsystems, their modules, and the 
overall control strategy and approaches used to de¬ 
sign and develop LEADER are described below. 

A Fna-ine Behavior Analysis 
Rnsed Diagnostic S ubsystem 

This subsystem performs diagnosis by analyz¬ 
ing the SSME behavior data representing various ab¬ 
normal conditions. SSME behavior is analyzed by 
learning the discriminant behavior-patterns associ¬ 
ated with various SSME faults, and using them to di¬ 
agnose unknown faults by matching the learned fault- 
patterns with unknown fault patterns. The results 
of this subsystem indicate that this subsystem per¬ 
forms fault diagnosis fast and successfully, if the un¬ 
known fault indeed represents one of the fault cate¬ 
gories present in the training data. This subsystem is 
comprised of the following modules: 

(1) The Machine Learning Module. 

(2) The Diagnostic Module. 

(3) Hypothesis Validation Module. 

(1) The Machine Learning Module 

This module uses the available SSME test sensor 
data and identifies several parameter characteristics 
which are important in diagnosing SSME faults rep¬ 
resented by the test data. These characteristics in¬ 
clude sensor value slope variations and inter and m- 
tra sensor slope-trend relationships, slope-magnitude 
relationships, temporal relationships, magnitude rela¬ 
tionships and duration relationships. Figure 2 shows 
the slope and temporal relationships. For example, if 
a sensor shows a positive trend whereas another shows 
a negative trend during the occurence of a fault, the 
two sensors are considered to have a slope-relationship 


Tvp-1 as shown in figure 2 . Aiso, if one ot the sensors 
starts showing the trend earlier than the other sensor 
and both continue to show the same trend upto end 
time, they are related by temporal relationship Typ-3 
as shown in figure 2. 

The machine learning program associates each 
fault with a set of relevant sensor characteristics. In 
our implementation thirteen types of SSME faults and 
their associated discriminant descriptors were iden¬ 
tified bv the machine learning module, using SSME 
sensors data obtained during the twenty-seven ground 
tests. The approach and the learning techniques used 
in developing this module are similar to those dis¬ 
cussed and described in detail in our previous paper 
[6], and will not be discussed in detail here. 

(2) The Diagnostic Module 

After the diagnostic patterns are learned by the 
MLM in the form of characteristic and discriminant 
descriptors, the diagnostic module is responsible for 
using those patterns in diagnosing a fault represented 
bv the given sensors’ data. In our present implemen¬ 
tation, the diagnostic module uses a neural network 
based pattern matcher that has been trained with the 
behavior patterns obtained from the machine learn¬ 
ing module. The network is trained to associate a 
particular set of relational sensor behaviors with a 
particular fault. During the diagnosis phase, the neu¬ 
ral network performs matching of patterns associated 
with the faults it was trained with, to the patterns as¬ 
sociated with the unknown fault. The patterns used 
in training the neural networks are binary representa¬ 
tions of discriminant sensor-pair relationships found 
to be relevant to a particular fault. The binary rep¬ 
resentation used here should not be confused with 
the standard binary value representation of a number, 
and is used only to indicate the presence or absence o 
a relationship. Further explanation on these patterns 
is -riven later in this section. As a result of this pat¬ 
tern analysis and matching, the system derives several 
possibilities of faults in terms of hypotheses and their 
estimated likelihood of occurences which are aval able 
for further analysis and validation by the hypothesis 
validation module. 

A neural network consists ot a system of inc^r- 
connected processing elements. The networks used in 
the present effort consist of three layers. The hrst 
layer consists of input nodes, the outputs of which 
are connected to a second layer of nodes. The sec¬ 
ond laver is called the "hidden” layer, because it is 
only indirectly affected by the application of external 
stimuli to the input nodes. The outputs of the hidden 
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nodes are connected to the third layer, which con 
of the output nodes of the neural network. 

Each connection is associated with a weight The 
input nodes are activated by an external stimulus de¬ 
rived from the behavior relations between various sen- 

" r pairs. The activations of the input nodes are mul- 
?• i-Ll Kv the weights of the connections between the 
££»YtSU *. resulting ^ 

2 input, to the hidden layer. The ««*«'«»» « 
tht hidden laver are then multiplied by the weight. of 
the connections between the bidden —- 
the resulting values become the inpu 

Jayel The activations of the output layer represent 

Zoning must be applied in order to deternune the 
most likely fault. 

Table 1 shows a few examples of the learned de¬ 
scriptors provided by the machine learning module. 
Thise descriptors were used to train the neural net- 
lorS The first column of table 1 shows pairs of 
related sensors, the second column shows the faul 
associated with the sensor pairs and the t 
umn shows the binary representations of 
descriptors based on the five types ofsensorrelati 
ships, namely, slope-type, slope-ma^ud 
magnitude and duration relationships 
the presence of binary 0 or 1 indicates the presence 

or absence of a corresponding relatl °“ h ' P pF ^l 
unple, the first descriptor for sensor pair HPFTDI 
V1CCCDT has a pattern 0001; since the 
from the left is on, this indicates the presence of a 
slope relationship of Typ-4 between the two sensors. 

The diagnostic module has been 
training it with the patterns 

^frTpSin; about thirteen fault ^ 

wpes of faults, it identified them wit 
curafiy Table 2 shows the results of the diagn 

tic module for the five test cases excluded from the 
tic module d to the diagnostic mod- 

STthe' "respond,ng fault classes. The 
of table 2 Shows the fault mstancesus^i 
to test the diagnostic module performance, 
cases were obtained from actual SSME ground tests 
and represented SSME faults. The second and third 
columns of table 2 show the fault hypotheses and 


the corresponding relative likelihood of occurences, 
respectively. These results form the basis for a pr_ 
inarv diagnosis by the BABS subsystem. The fault- 
hypothesis with highest vote-count i» ^ch case shown 
in table 2 is considered to be the most likely fault cor 
responding to the sensor data of the test-case, ur 
ther vahdation of the preliminary results is performed 
as described in the next section. 

(3) Hypothesis Validation Module (HVM) 

This module is responsible for validating the fault 
hypotheses resulting from the diagnostic module in 
the light of SSME design. It determines if the events 
associated with a hypothesis and deser.bed m term 
of SSME design indeed occured, and are supp 
the sensor data of the fault currently being diagnosed 
^ the diagnostic module. If the events associated 
with a hypothesis are all found to have occured as de¬ 
duced by direct or inferred evidence, the P °" is ' 

L hypothesis is considered to be a valid diagnose, 
representing the fault currently being diagnosed. 

The strategy used by the hypothesis validation 

module involves knowing the events associated with 

the hypothesis being validated and using t e p 
^ss«of event-justification and event.con.tttd.c.mn 

^determine whether the events assoc,ated with a h>- 
oothesb^have occuted ot not. The processes of event- 
justification and event-contradiction use ‘he miknown 

fault sensor data and SSME design to infer the knowl 
Si or evidence that would confirm or deny the oc¬ 
curence of events associated with the hypothesis being 

validated. 

The process of event-justification proceeds as fol- 

An event-list associated with the hypoth«is b> 

. vA*t*A i« accessed. Events in the event-list 

^, a g^r description of the way the faulty 

edge and parameter-translation knowledg . 
nroaches for event-justification are used by the even 
proacnes lor e c j Hireet-iustiftcation ap- 

instification process, namely .a direct j 

pr^h. audio indirect-ju^i^appro^iThe 
direct-justification approach involves using 

unknown-fault sensor ^“^"“‘^“/hvpotbesis be- 

^ir«r:^t:i^to„tradic,- 


lows 
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ing evidence is obtained directly are placed under a 
contradicted-events list, and are processed later by 
the event-contradiction process. On the other hand, 
all those events for which no supporting or contradict¬ 
ing evidence is found are processed using the indirect- 
justification approach. The indirect-justification ap¬ 
proach uses SSME design knowledge and parameter- 
translation rules to find other ways of finding support¬ 
ing or contradicting evidence. For example if an event 
suggests that temperatures in the proximity of a com¬ 
ponent should go up, but does not explicitly mention 
the sensors, then the components in the proximity of 
the component are located from the SSME design and 
their temperature sensors are obtained and checked 
for a positive trend in the unknown fault sensor-data. 
An example illustrates parameter translation: when 
an event expects the temperature of hot-gas in a com¬ 
ponent to increase, the pressure of hot-gas should be 
checked if the temperature sensor does not indicate 
a rise in temperature. If the pressure of the hot- 
gas seems to be rising, it should be considered as a 
supporting evidence since a temperature rise can be 
translated to a pressure rise. 

The Event-contradiction processing involves 
finding out whether the events contradicted belong 
to the upstream, downstream, proximity, or within 
a component, and determining the nature of contra¬ 
diction. Contradictions of events expected to occur 
within a component (for a given fault-hypothesis to 
be true) are not tolerated, although the contradicted- 
events belonging to upstream or proximity may be 
considered more tolerable than the contradictions of 
events belonging to downstream or within a compo¬ 
nent indicated faulty by the hypothesis being vali¬ 
dated. As a result of determination and short analy¬ 
sis of the locations of contradicted events, the event- 
contradiction process gives a list of intolerable con¬ 
tradictions and a separate list of contradictions that 
can be overlooked. 

If none of the hypotheses is found to be valid 
based on the validation process mentioned above, the 
occurence of a new type of fault may be indicated. In 
cases of new types of faults the deep-level reasoning 
based^ubsystem is most likely to succeed in perform¬ 
ing fault diagnosis using any available results from 
The Hypothesis Validation Module. 

Figure 3 illustrates some of the results obtained 
from the HVM module for a test case in which the 
module was given the sensor data of the fault cur¬ 
rently being diagnosed, and three fault-hypotheses 
to be validated. The three fault hypotheses repre¬ 


sented Main Combustion Chamber Manifold (MCC- 
MANIFOLD-FAULT) fault. Liquid Oxygen Injector 
(MCC-LOX-INJ-FAULT) fault, and High Pressure 
Fuel Pump (HPFP-FAULT) fault. As shown in figure 
3 the occurences of various expected events within, 
upstream, downstream, and/or proximity of faulty 
components represented by hypotheses MCC-LOX- 
INJ-FAULT and MCC-MANIFOLD-FAULT are con¬ 
tradicted, thus invalidating the occurences of these 
two hypotheses for the sensor data of the fault cur¬ 
rently being diagnosed. The HPFP-FAULT hypothe¬ 
sis does not have any contradicted events in any of the 
four categories shown in figure 3. These results from 
the HVM module are correct because the unknown 
fault sensor data was obtained from ground test data 
representing the HPFP fault. 

R F.ngina Desig n Analysis Based 
Diagnosti c Subsystem (DABS) 

This subsystem performs diagnosis using design 
based reasoning. In order to cope with the temporal 
engine dymanics, the concept of Time Frame (TF) is 
used The TF concept is particularly important be¬ 
cause abnormality and trend variations are indicated 
by different sensors in a unique temporal sequence for 
each fault. A TF is that particular instance of time 
when one or more sensors start indicating an abnor¬ 
mal parameter condition. Each TF is given a time 
frame number which is incremented by one for the 
successor time frame and so on. Also associated with 
each TF is the group of sensors that start showing 
abnormal engine parameters and trends at that TF. 
Every time one or more sensors start indicating an 
abnormal condition of the engine, a new TF is cre¬ 
ated and given a time frame number higher than the 
previous TF number by one. Figure 4 shows the time 
frames created for simulating real time sensor indica¬ 
tions towards abnormal engine parameters and their 
trends for a High Pressure Fuel Pump Inlet block¬ 
age fault. The time frame numbers are indicated by 
arrows, and the sensors indicating abnormal parame¬ 
ters and their trends are mentioned in a list along the 
trend portion of the corresponding sensor data curve. 
In order to perform diagnosis in a real-time environ¬ 
ment, information obtained from the sensors indicat¬ 
ing abnormal engine parameters and their trends is 
mapped into various TFs before presenting it to the 
deep-level reasoning module. 

The DABS subsystem consists of the following 
modules: 

(1) SSME Design and Functional Module (DIM). 

(2) Deep Level Reasoning Module (DLRM). 
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(3) Integration Module (IM). 

(1) SSME Design k Functional Module (DFM) 

This module is responsible for providing the de¬ 
sign and functional knowledge pertaining to SSME 
subsystems, components, or subcomponents when¬ 
ever requested by the other modules of LEADER. The 
SSME design and functional module incorporates an 
SSME design and functional knowledge base (DFKB). 
The major SSME subsystems such as the oxidizer- 
subsystem, the fuel subsystem, the main combustion 
subsvstem, their major components and subcompo¬ 
nents along with their functional and physical re¬ 
lationships, interactions, and associated sensors are 
represented in DFKB. Object representation is used 
to represent each subsystem and its major compo¬ 
nents and subcomponents. For example, the oxidizer- 
subsystem is defined as an object which, in addi¬ 
tion to the subsvstem-level knowledge, points to com¬ 
ponent objects such as the High Pressure Oxidizer 
Turbine (HPOT). Component-level objects, besides 
storing component-level and some of the parent level 
knowledge, also point to their subcomponents such 
as the main oxidizer pump. Thus, a structural hi¬ 
erarchy of knowledge is established. Figure 5 shows 
a schematic of a Space Shuttle Main Engine design. 
The major components have been indicated. Figure 
6 shows the component-object representation of High 
Pressure Oxidizer Pump (HPOP). 

The various object attributes include pointers to 
the parent subsystem, children subcomponents, up¬ 
stream, downstream, proximity, and high interaction 
components. Some of the object attributes shown in 
figure 6 indicate functions of HPOP in terms of the 
change in temperature and pressure of liquid oxygen 
at the input and output of HPOP. Other attributes 
describe the types of fault-primitives that can occur 
in HPOP, the associated sensors and other relevant 
information. 

DFKB also includes a number of causal-symptom 
rules on the mechanism involved with SSME com¬ 
ponent design, performance, and interaction. The 
causjj-symptom rules have been divided into five 
"roups depending on the nature of the physical re¬ 
lationships described by them, as discussed below. 

The physical relationships between a basic com¬ 
ponent such as a turbine and a substance such as hot- 
gas which it carries have been represented in the form 
of causal-symptom relationship rules involving the 
heuristics and the physics that describes the cause(s) 


for a certain sensor indication towards an abnormal 
engine parameter and its trend. The causes may be 
described in terms of other engine parameters and 
their trends, or in terms of fault-primitives such as 
a leak, blockage, burn, structural-failure occuring in 
one or’more components or in upstream, downstream 
or within the reference component(s). An example of 
such a rule is: 

SYMPTOM (Rising Turbine Speed) 

CAUSE (Rising Turbine Hot-gas Temperature) 

This rule suggests that the cause of an increase 
in turbine speed may be attributed to a rise in the 
temperature of hot-gas flowing into the turbine. 

The physics behind thermal/hydraulic parameter 
trends shown by a substance such as fuel or hot-gas 
has also been described in the form of causal-symptom 
laws/rules. For example, whenever the hot-gas tem¬ 
perature increases beyond the one desired in SSME, 
a basic cause is the fuel and liquid oxygen ratio im¬ 
balance in the hot-gas, with more liquid oxygen and 
comparatively less fuel in the hot-gas mixture. Ail 
example of such a rule is: 

SYMPTOM (Rising Hot-gas Temperature) 

CAUSE (OR (Increasing Liquid-Oxygen Flow)) 

(Dropping Fuel Flow) 

This rule suggests that a trend towards a rise in 
hot-gas temperature in SSME may be attributed to 
either an increase in the flow of liquid-oxygen or a 
drop in fuel supply/flow in the engine. 

Another set of causal-symptom rules/laws based 
on the mechanics and design of SSME and its com¬ 
ponents, describes the interaction between any two 
components used in SSME design. For example, High 
Pressure Fuel Turbine and High Pressure Fuel Pump 
are two components of types Turbine and Pump, re¬ 
spectively, used in SSME design. This allows deduc¬ 
ing the influences that components with high interac¬ 
tion in a design like SSME would have on each other. 
These laws also allow inferring the possible causes for 
a component’s parameter trend in terms of other com¬ 
ponents interacting highly with the symptom param¬ 
eter. For example, a decrease in load on a pump such 
as the High Pressure Fuel Pump can be one of the 
possible causes for the increase of speed of its driving 
turbine, the High Pressure Fuel Turbine. An example 
of a causal-symptom rule representing this group is as 
follows: 

SYMPTOM (Dropping Turbine Speed) 

CAUSE (Increasing Pump Load) 
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This rule suggests that a possible cause for a 
trend towards a drop in turbine speed may be an in¬ 
creasing load on the driven pump. 

In order to be able to infer the possible causes or 
predicted-effects of a parameter trend in terms of its 
upstream or downstream component, heuristic rules 
describing the causes and predictions in terms of up¬ 
stream and/or downstream component connections in 
SSME design are used. These heuristics are based on 
the type of a sensor such as a temperature, pressure, 
speed/rpm, or a flow measurement sensor, as well its 
trend. For example, an increase in pressure of a liquid 
such as fuel would indicate that there is a blockage in 
the downstream component(s). These heuristics are 
more of a general nature and depend on the type of 
the given parameter. An example of such a causal- 
symptom rule is as follows: 

SYMPTOM (OR (Dropping Liquid Flow)) 
(Dropping Liquid Pressure) 

CAUSE (OR (Leaking Liquid Upstream)) 
(Blocked Liquid Upstream) 

This rule suggests that the possible causes for a 
symptom like a drop in flow or pressure of a liquid 
in SSME may be the occurence of a liquid leak or 
blockage upstream of the location where the symp- 
toms have been indicated. 

Finally, another set of domain-specific heuristics 
has been used, which indicates the cause(s) for a gen¬ 
eral scenario presented by the parameters showing 
certain trends. An example of such a heuristic is 
shown by a causal-symptom rule below. 

SYMPTOM (AND (Rising Temperature Sensors)) 
(Dropping Pressure Sensors) 

CAUSE (OR (Leaking Fuel)) 

(Blocked Fuel) 

This rule suggests that if most of the tempera¬ 
ture sensors in SSME show rising trends, whereas the 
majority of pressure sensors indicate a droping trend, 
the likely cause of the indicated abnormality is fuel 
leak or blockage in SSME. 

(2) Deep-Level Reasoning Module (DLRM) 

This module is responsible for performing the de¬ 
sign based forward reasoning to infer the causes, lo¬ 
cation and associated events of a given unknown fault 
by utilizing the unknown fault sensor data, SSME de¬ 
sign and functional knowledge, and domain heuristics. 
The deep level reasoning module is comprised of three 


problem-solvers, namely the subsystem level problem 
solver (SLPS), the component level problem solver 
(CLPS), and the subcomponent level problem solver 
(SCLPS). It uses a hypotheses blackboard to write 
the intermediate results, such as most likely faulty 
subsystem(s), faulty components, associated events, 
and supporting sensors. A hierarchy of known fault 
hypotheses and associated subsystems or components 
is also maintained and used by the various problem 
solvers of this module. 

Our approach in developing the three problem 
solvers is focussed more towards identifying the root 
causes as well as the affected subsystems, compo¬ 
nents, and subcomponents, that are responsible for 
the abnormal behavior of SSME, and predicting fu¬ 
ture faults. This approach is adopted because aaer 
a fault occurs, its symptoms and affects propagate 
throughout the system as time passes. Thus, if the 
cause(s) of a fault is known, the fault propagation may 
be contained or stopped by taking corrective mea¬ 
sures to eliminate as much as possible the causes of 
the fault. The three problem solvers mentioned above 
work interactively to identify and diagnose new faults. 
Each of these problem-solvers is discussed in detail 
below. 

The purpose of the subsystem-level problem- 
solver is to identify the subsystems appearing to be at 
fault, and prioritize them in the order in which they 
should be further investigated to identify the specific 
components in which the fault has occurred. As soon 
as the information associated with a TF is presented 
to the DLRM, the SLPS is triggered. It uses certain 
heuristics on sensor data. For example, in general 
the speed, vibrational, and flow measuring sensors are 
more reliable in pointing towards faults than pressure 
and temperature measuring sensors. Heuristics like 
these allow the SLPS to attach weights or varying 
degrees of importance to the various sensors belong¬ 
ing to a TF. From the SSME design and functional 
module, the location of each sensor associated with 
the given TF is obtained in terms of the subsystem 
they belong to. The total weight of all the sensors 
of a SSME subsystem indicating an abnormal engine 
condition gives that subsystem’s fault-likelihood num¬ 
ber. The various subsystems, corresponding fault- 
likelihood numbers, and associated sensors and TFs, 
are written on the subsystem-level of the hypothesis- 
blackboard. The subsystem with the highest fa un¬ 
likelihood number is then investigated further by the 
component-level problem-solver. 

The major objective of the component-level 
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problem-solver is to identify the components at fault 
and the root-causes that led to the SSME fault(s). 

It generates a fault-mechanism tree to achieve this 
objective, and also attempts to predict the compo¬ 
nents affected and conditions that will result due to 
the present fault(s). A fault-meachanism tree (FMT) 
is employed. The FMT is a tree-structure that de¬ 
scribes the mechanism of fault(s) in terms of the ba¬ 
sic underlying causes, effects, and predictions. The 
FMT consists of causal nodes at top, pointing down¬ 
ward to the corresponding symptom nodes, and the 
symptom nodes pointing towards the predictions or 
future symptom nodes. A symptom-node can be con¬ 
sidered as a causal-node relative to the node pointed 
to at the lower-level. Apart from the SSME de¬ 
sign and functional knowledge, the causal-symptom 
rules of physics, mechanics, and heuristics behind the 
various SSME components and parameters are uti¬ 
lized to generate the FMT. Figures 7a, 7b, <c illus¬ 
trate various stages of development of a FMT by the 
Deep>-Level Reasoning Module when presented with 
the sensor data representing "Blockage of High Pres¬ 
sure Fuel Pump Inlet" fault. Each node represents a 
mechanism-event, the TF when the mechanism-event 
most likely occured, and the likelihood referred to 
as Confidence-Level (CL) that the node lies on the 
specific path of the FMT which explains the causes 
and consequences of the fault currently being diag¬ 
nosed with highest confidence-level values at each 
node of the path. Such a path is referred to as the 
best path of the FMT. Various nodes of the FMI 
are created during the process of reasoning. Fig¬ 
ure 7a shows the FMT developed during the TF 1. 
Each node of the FMT is shown by a block, with 

the causal/symptom/prediction event represented in 

the form of a list. The CL associated with each node 
and the TF during which the event(s) associated with 
the node may have occurred is represented by the 
CL value and the TF is indicated with the node. In 
figure 7a, the events associated with the top nodes 
may be responsible for triggering the events associ¬ 
ated with the lower level nodes. Using pointers from 
parent(s) to children, a hierachy is maintained, with 
events associated with the parent nodes giving birth 
to the events associated with corresponding children 
nodes. The top right node in figure 7a suggests that 
a drop in fuel flow to HPFP may be a possible cause 
for a decrease in HPFP pump load as indicated by 
the corresponding child node on the right. A fault- 
event prediction made during the TF 1 is shown by 
the node at the bottom of figure 7a. The predicted 
fault-event indicates an increase in speed of the HPFT 
turbine. The confidence-level of each node is calcu¬ 
lated by adding the CLs of its children nodes and in¬ 


crementing the resulting CL value by one, the initial 
CL value of a node. The FMTs in figures 7b, and 7c 
can be interpreted similarly. The node with no point¬ 
ers coming into or out of it in figure 7b does not sup¬ 
port any of the previously occured events, or no cause 
has been obtained corresponding to the event associ¬ 
ated with this node. Also, this node does not have 
any child, which indicates that it has not been found 
to cause any other fault-events. Figure 7c shows (in 
thick lines) the best path of the FMT. The best path 
in figure 7c may be formalized into a hypothesis in¬ 
dicating that the HPFP fault is due to a drop in the 
flow of fuel to HPFP. 


Triggering of the CLPS performs the following 
steps: 

- It reads the subsystem-level part of the 
hypothesis-blackboard to obtain the subsystem with 
the highest likelihood of fault, and the sensors associ¬ 
ated with the latest TF. If no such sensor is available, 
then it either goes for the next likely faulty subsystem 
or simply returns, depending on the relative likelihood 
of the various subsystems to be at fault. 


- Once the sensors obtained during the latest 
TF and supporting the highest priority subsystem 
are read, the location of the sensors in terms of 
the components of this subsystem are obtained from 
the design knowledge. Also, the substance(s) associ¬ 
ated with each of those component-names is obtained 
from the SSME design and functional module. The 
component-type such as turbine, pump, injector is 
obtained for each component-name. 


- Knowing the components, the corresponding 
omponent-types, the associated active sensors and 
rends, and wherever applicable the associated sub- 
tance with each component, one or more lists follow- 
ng certain formats called templates are created. The- 
ies templates mostly represent symptoms of the fau 1 
:urrently being diagnosed, and are matched with the 
symptom side of appropriate causal-symptom rules 
jf the DFKB to heuristically generate the possible 
:auses and future predictions from the sensors and 
trends of the fault currently being diagnosed tem¬ 
plates representing the symptoms of the fault cur¬ 
rently being diagnosed are called symptom-templates. 
The five types of symptom templates defined by us 
and understood by the DLRM are described later ... 
this section. During this step, the CLPS creates a 
possible symptom-templates for each set comprising 
of a related sensor, its trend, indicated componen - 
type, and substance, as applicable follow,ng the for¬ 
mat of the five symptom-template. 
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- Templates created from related information 
comprising an abnormal parameter and its trend, 
component-type, and substance (as applicable), are 
matched with the symptoms encoded in the causal- 
svmptom rules, preferably of a high success rate cat¬ 
egory known for that template-type. In this way a 
number of possible causes for each symptom-template 
is generated, and pointers set from the causes to¬ 
wards the matching symptoms-templates. When¬ 
ever a symptom-template fails to match any causal- 
svmptom rule, an attempt is made to match it with 
the causal side of the causal-symptom rules. When¬ 
ever a match is found with the causal side of a rule, 
the corresponding symptom-part is considered as a 
prediction, whose other (new) possible causes are then 
searched by matching this prediction template with 
the symptom-sides of the causal-symptom rules. 

- The above step is carried out for each set of 
related information comprising an abnormal param¬ 
eter and its sensor, trend, component, component- 
tvpe, substance, for the sensors indicating abnormal 
behavior of SSME during the latest TF. 

- As a result of the above steps a fault-mechanism 
tree/hierarchy is created, which is used to find the 
most likely search path leading from the cause(s) of 
fault(s) to the symptoms and onto the predictions 
wherever the chain points to the symptoms expected 
to occur after the latest TF. 

- The above steps are repeated for relevant com¬ 
ponents every time another TF information is pre¬ 
sented to the CLPS, and any additional causes and 
symptoms generated are added appropriately to the 
fault-mechanism tree created so far. 

- The next step is to find the best path in the 
fault-mechanism tree. Search for the best path is per¬ 
formed by choosing the path with highest confidence- 
level value at each node of the path. The best path 
of the fault-mechanism is selected and formalized into 
a component-level fault hypothesis. Component-level 
results from the CLPS are written on the component- 
level part of the hypothesis-blackboard. 

Iih. our implementation five types of symptom- 
templates have been defined, each of them may have 
a higher rate of success in matching with a certain 
category/set of causal-symptom rules and less with 
other categories of causal-symptom rules. The five 
types of symptom-templates are as follows. 

Template-typ-1: (Trend Component-typ 
Substance Parameter-typ) 


Template-typ-2: (Trend Component-typ 
Parameter-typ) 

Template-typ-3: (Trend Substance 
Parameter-typ) 

Template-typ-4: (Majority-trend 
Parameter-typ) 

Template-typ-5: (Trend Parameter-typ) 

Template types 1, 2 have a higher rate of match¬ 
ing with causal-symptom rules on the interaction be¬ 
tween two components, or between a component and 
the substance contained in it. Template type 3 has 
more successful matching with causal-symptom rules 
on substances (liquid/mixture). Template type 4 has 
a higher rate of successful matching with the general 
domain heuristics, and Template type 5 has more suc¬ 
cess with the rules on sensor-behavior heuristics. 

After diagnosing the component-level faults, the 
DLRM triggers the SCLPS in order to investigate 
the subscomponent-level faults corresponding to the 
already identified component-level faults. Knowing 
the components at fault, and the description of the 
component-level fault(s) such as causes, symptoms, 
and associated plausible fault-primitives, the SCLPS 
looks into the design of each of the components. 
From the design of each faulty component, infor¬ 
mation on corresponding subcomponents and fault- 
primitives possible with those sub-components is ob¬ 
tained. Knowing the nature of a fault that oc¬ 
curred inside a component, its subcomponent with 
the best match between the nature of component 
fault, and the subcomponent function, location, de¬ 
sign,and plausible fault-primitives associated with its 
subcomponents is chosen to be the faulty subcompo¬ 
nent. The described approach for the SCLPS is being 
implemented in LEADER. Results from the SCLPS 
are written on the subcomponent-level part of the 
hypothesis-blackboard. 

The overall fault-hypotheses are formalized and 
a brief explanation is then generated from the best 
path(s) of resulting FMT and the results on the three 
levels of the hypothesis blackboard. 

r Integration Sub system US! 

The diagnostic results based on the SSME behav¬ 
ior and design analyses are integrated by the Integra¬ 
tion Subsystem and the resulting hypothesis is consid¬ 
ered final. IS checks if the results obtained from both 
the design and behavior analyses based approaches 
concur, indicating the same hypothesis. If they con- 
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cur, then the hypothesis is considered final: other¬ 
wise, a brief analysis of the certainty of results from 
each approach is carried out, so that the subsystem 
either combines both the hypotheses or validates one 
of them to give the final results. 


RESULTS 

The results of fault diagnosis obtained from 
LEADER illustrate the diagnostic power which 
results from integrating two complementary ap¬ 
proaches, one based on the analysis of engine design 
and the other on the analysis of the engine behav¬ 
ior. We have tested LEADER with a limited amount 
of rocket engine data and expect to test it further 
for numerous types of SSME faults as more data on 
SSME becomes available. The approach discussed for 
the DLRM has been tested for some of the fuel sub¬ 
system components’ faults and the main-combustion 
subsystem components’ faults. The FMTs developed 
for each fault case has been correctly analyzing the 
sensor data of the test case, but in most cases the 
speed and reliability of the component-level fault di¬ 
agnosis depends on the accuracy and reliability of the 
causal-symptom rules in the DFKB. We expect to im¬ 
plement the SCLPS and the integration module using 
the approach discussed here. 

Although we are still far from complete and ex¬ 
haustive testing of LEADER, the results obtained so 
far with the available data show that the approaches 
discussed in this paper have a potential for perform¬ 
ing real-time fault diagnosis in a complex mechanical 
system such as SSME. 
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Fig. 7a - Fault Mechanism Tree developed 
during Time-Frame 1 shown in 
figure 4 



Fig. 7b - Fault Mechanism TVee resulting 
after Time-Frame 2 shown in 
figure 4 



Fig. 7c - Fault Mechanism Tree resulting 
after Time-Frame 3 shown in 
figure 4 
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MICRO GRAVITY FLUID MANAGEMENT 

• Two-Phase Flow in Helical and Spiral Coils 

• The Subscale Orbital Fluid Transfer Experiment 


• Low Gravity Transfer Line Chilldown 




Project Title: Two-Phase Flow in Helical and Spiral Coils 

Principal Investigator: Edward G. Keshock 

Other Investigators: Mia L. Bush, Adel Omrani and An Yan 


Project Report 


Coiled tube heat exchangers involving two-phase 
flows are used in a variety of application areas, ex¬ 
tending from the aerospace industry to petrochem¬ 
ical, refrigeration land power generation industries. 
The optimal design in each situation requires a funda¬ 
mental understanding of the heat, mass and momen¬ 
tum transfer characteristic of the flowing two-phase 
mixture. However, two-phase flows in lengths of hori¬ 
zontal or vertical straight channels with heat transfer 
are often quite difficult in themselves to understand 
sufficiently well to permit accurate system designs. 

In the case of possible applications to systems 
operating in a space environment, coiled flow chan¬ 
nels may be advantageously employed, since the cen¬ 
trifugal forces experienced by the two-phase fluid acts 
to promote contact with the channel wall, thereby 
tending to insure good contact with the wall, and 
an enhanced heat transfer. Because of such possi¬ 
ble improved heat transfer performances it becomes 
important to identify the conditions under which heat 
transfer rates are optimized. Correspondingly, it is of 
importance to know the “price” that must be paid 
for such optimal performance. From a fundamental, 
physical viewpoint, these questions are resolved to a 
very large extent by (1) prediction/identification of 
the two-phase liquid-vapor flow patter, and (2) the 
accompanying pressure gradient to achieve or mam- 
tain specific flow patterns. 

In the case of horizontal or vertical straight flow 
channels a variety of flow pattern maps have been 
generated and proposed for application in a 1-g en¬ 
vironment [1-4]. Similarly a variety of methods for 
predicting two-phase pressure drops have been devel¬ 
oped (see [5], for example). Relatively few studies 
have been made of helical coil flows, however [6-9]. 
The most recent significant study of helical two-phase 
flows is that of Uddin [10] who observed a variety 
of flow patterns under various flow conditions, and 
generated a map identifying the various flow regime 
boundaries (Figure l). The transition are considered 
to be dependent upon the interactive effects of (1) the 


gravitational force of the test environment, (2) the in¬ 
duced centrifugal force field, and (3) shear stresses in¬ 
duced by secondary flow within the vapor phase por- 
tion of the flow. 

In general, two-phase flow in coiled geometries 
can be expected to be much more complex than in 
straight channels because of the additional factors of 
an induced centrifugal force field and the resulting in¬ 
duced secondary flows. Ones intuition would lead one 
to expect that the liquid phase would produce a phase 
separation, with the liquid phase being centrifuged to 
the outer wall. However, the increased complexity 
of helical flows becomes apparent from the studies of 
Bannerjee et al. [6], who observed that under some 
conditions the liquid phase flowed along the tube wail 
nearest the coil axis - a phenomenon labeled film in- 
version.” 

A second factor tending to complicate coil flows 
is the inducement of secondary flows, as depicted 
schematically in Figure 2. The secondary flow de¬ 
picted in Figure 2(a) would be expected even in a 
single-phase flow. The pattern depicted m Figure 
2(b), proposed by Uddin, indicates an annular liquid 
film surrounding a vapor-core flow in which a rota¬ 
tional secondary flow pattern is set up. It should be 
noted that this latter flow pattern has not been exper¬ 
imentally verified, but has only been hypothesized. 

Sundstrand Corporation developed the basic 
concepts for two-phase helical flows and the present 
project is an extension of experimental work con¬ 
ducted at Sundstrand. 

The present study has the following general oh- 
jectives: 

1. Observe two-phase flow patterns of air-water and 
R-113 working fluids over a range of flow condi¬ 
tions, for helical and spiral coil geometries, of 
circular and rectangular cross-section. 

2. Compare observed flow patterns with predictions 
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of existing flow maps. 

3. Study criteria for flow regime transitions for pos¬ 
sible modifications of existing flow pattern maps. 

4. Measure associated pressure drops across the 
coiled test sections over the range of flow con¬ 
ditions specified. 

The experimental measurements and observa¬ 
tions will be conducted using an apparatus donated 
by the Sundstrand Corporation, and now located and 
installed in the Two-Phase Flow Laboratory of the 
Mechanical and Aerospace Engineering Department, 
UTK. A schematic diagram of the modified apparatus 
is shown in Figure 3. 

Flow regime observations will be made using a 
high shutter speed video camcorder. Additionally, 
several selected operating points will be observed us- 
ing a high-speed motion picture camera. 

Pressure drop measurements will be made for 
each operating condition, and compared with values 
calculated using several methods from the literature 
intended for use in straight channel flows. 

The system is now operational and the first set 
of data has been only recently obtained for a heli¬ 
cal coils glass test section of 0.325 inch I.D., and 3 
inch diameter helix axis, with a pitch of 1.” Future 
tests are planned for a 0.325” I.D. spiral coils section, 
and three rectangular cross-section of three different 
aspect ratios. 
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Project Title: 


The Subscale Orbital Fluid Transfer Experiment 


Principal Investigators: Dr. J. S. Meserole and Dr. Frank G. Collins 

Other Investigators: Mr. Ogden Jones, Dr. Basil Antar, Dr. Reinhard Menzel 

and Mr. Perry Gray 


Project Report 


PROJECT O BJECTIVES 

CASP is a subcontractor to Boeing Aerospace 
Corporation to provide support for the concept defi¬ 
nition and design of a subscale orbital fluid transfer 
experiment. 

fVSQ ACC OMPLISHMENTS 

The Subscale Orbital Fluid Transfer Experiment 
(SOFTE) contract has three phases: concept defini¬ 
tion, preliminary design and final design. Most of 
the work during this year concentrated upon the first 
phase. 

SOFTE is an experiment that will look at the 
fluid mechanics of the process of transfer of a satu¬ 
rated fluid between two containers. The experiment 
will be placed in two GAS can containers; the tanks 
will be in one container and the power and electron¬ 
ics will be in a second container. Since GAS cans are 
being used, the experiment will be autonomous. 

The fluid supply tank will contain a liquid acqui¬ 
sition device (LAD) to assist in fluid removal. The 
LAD will be in four bands, each containing wire 
screens. The receiver tank will contain a elastomeric 
diaphragm that can be inflated to push the fluid back 
to the supply tank. Vaporized fluid will be used to 
inflate the elastomeric diaphragm. The experiment 
will be designed to investigate the no-vent fill process 
of a cryogenic fuel tank, LAD breakdown, LAD refill 
and vented fill. 

The work during the present year consisted of 
examining concepts for visual observation of the fluid 
transfer process, examination of methods for accu¬ 
rately metering the amount of fluid transferred be¬ 
tween the two tanks, examination of possible test flu¬ 
ids, and consideration of the materials to use for the 
elastomeric diaphragm. 


The objective of the visual observation is to lo¬ 
cate the fluid-vapor interfaces and. if possible, quan¬ 
tify the amount of vapor and the area of the interface. 

It is proposed to use video cameras to view the over¬ 
all process in each tank and to place borescopes or 
other devices through the tank walls to obtain de¬ 
tailed. undisiorted views inside the tanks of critical 
portions of the transfer process. Further work will 
continue to find an economical means for providing 
this detailed view, which clearly would increase the 
data obtained from the experiment. 

Fluid gauging of cryogenic fluids at 0-g remains 
one of the greatest challenges for orbital resupply ot 
cryogenic liquids. In the present experiment the flow 
meter is the primary means of gauging the fill level o 
the receiver tank and should be as accurate as pos¬ 
sible In addition, the meter must be rugged enough 
to survive lift-off, have a small pressure drop and be 
of small size. It is desirable to measure mass flow di¬ 
rectly because of the expected density change of the 
fluid caused by the temperature change during the 
course of the experiment. No commercial meter met 
these requirements and it was determined that the 
best compromise would be a turbine meter. The tur¬ 
bine meter can accurately measure volume flow if the 
fluid passing through the meter is known to always 
be in a liquid state. To measure mass flow, addi¬ 
tional pressure and temperature measurements must 
be made. The overall expected accuracy and reliabil¬ 
ity of this system will be examined during the next 

phase. 

Numerous fluids having a vapor pressure in the 
range 0.5 to 0.75 atmosphere were examined. The 
base-line fluid was Freon 113, which is compatible 
with the tank materials although it has a vapor pres¬ 
sure that is lower than desired. No safe alternative 
could be found; however, the duPont Freon alterna- 
tives are still under consideration. 

During the next phase the preliminary design ot 



the optical system and the elastomeric diaphragm will 
be completed. Abo, the overall accuracy and reliabil¬ 
ity that can be achieved with the turbine flow meter 
will be determined. Thb information will be incor¬ 
porated into the overall experiment design, which is 
being conducted by Boeing Aerospace Corporation. 
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Low Gravity Transfer Line Chilldown 

B. N. Antar and Jere S. Meserole 

F. G. Collins, A. Hedayatpour and T. McGee 


Project Report 


eij pnificant Acc o mplishments m the East Yga i 

A code has been developed that solves for the 
transfer line chilldown time and flow and heat trans¬ 
fer characteristics in one-g environment. The code 
solves the transient, one dimensional, space average 
mass, momentum and energy conservation equations 
for liquid-vapor two-phase flow in tubes. The phys¬ 
ical configuration solved is that appropriate for bot¬ 
tom coolant injection in a vertically supported heated 
tube. Four distinct regions are considered consecu¬ 
tively: fully liquid, inverted annular, dispersed, and 
fully vapor flow. The conservation equations for both 
the liquid and the vapor are solved in each region 
separately. Also, in each region the mass and energy 
transport between each phase as well as the energy 
and momentum transport between the tube wall an 
the fluid are accounted for. A finite wall thickness is 
also considered. 

The model described above was solved numeri¬ 
cally through a mixed finite difference scheme with 
forward time marching. The inverted annular regime 
was resolved using a semi-implicit finite differencing 
while the dispersed regime was solved explicitly. Also, 
a staggered mesh was used in which the velocity was 
resolved at mesh boundaries while all other field vari¬ 
ables were resolved at the mesh centroids. Different 
mesh sizes were used depending on the region o so¬ 
lution. A coarse mesh was used in the dispersed flow 
region while a much finer mesh was used m both the 
inverted annular flow region and the tube walls. 

The code was written and verified against the re¬ 
sults of Kawaji and Banergee (J. Heat Transfer vol. 
110 1988). In that work a model was developed tor 
a vertical tube chilldown with bottom flooding with 
water. They validated their mode! with an experi¬ 
ment specifically constructed for that purpose. The 
results from our code compared favorably with the 
results of Kawaji and Banergee for the case of a U 
ft long and 0.5 in I.D. heated tube at 1000 °F with 
water injected at 30 °F. Two inlet water speeds were 


used 5 in/sec and 3 in/sec. The results of both.cases 
compared well with both the experiments of Kawaji 
and Banergee and their calculations. Figures 1 and 
2 show wall temperature histories at four locations 
along the tube axis for both flow rates. Figures 3 and 
4 show the void fraction for the two flow rate cases at 
three axial locations. 


jbcus nf Current Research Activities 

The code validations will be completed by com- 
taring other field variables (e.g. vapor velocity, heat 
luxes, etc.) with the experiments and numerical anal¬ 
yses of Kawaji and Banergee. Next, the code will be 
tested for crvogen chilldown characteristics by simu¬ 
lating an a liquid nitrogen chilldown experiment. In 
that experiment the cooling fluid will be liquid nitro¬ 
gen while the test section will be a O.o m O.D. by 
12 ft long copper tube at ambient temperature. The 
tube will be supported vertically with liquid nitrogen 
injected at the bottom of the tube while the top is 
open to the atmosphere. A test stand for this mode 
is presently under construction. This experiment will 
be run and numerically simulated for different liquid 




' 1_ax ae 


Plans for th e Next Year 

For the next year the plans are to validate the 
de with liquid nitrogen, bottom injection experi- 
ents. Also, the code will be modified to handle the 
se for which the tube will exit into a receiver tank at 
iy prescribed temperature and pressure. The code 
ill also be modified for low gravity conditions. For 
is end a series of low gravity aircraft experiments 
ill be performed to gather all the necessary data 
:eded for upgrading the code. The first '°w gravity 
cperiment is scheduled for flight in June, 090 and 
is hoped it will soon be followed by the next two 
ights. In the first flight experiment the working fluid 
ill be freon 113 and the test section will be O.o 
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in I.D. and 4 ft long quartz tube heated to approxi¬ 
mately 300 °C. In the following two flight experiments 
liquid nitrogen will be used for the working fluid. 






ELECTRIC PROPULSION 


• Electrostatic Plasma Accelerator (EPA) 

• Ion Propulsion 

• Magnetic Annular Arc (MAARC) 
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Electrostatic Plasma Accelerator (EPA) 


Project Title: 

Principal Investigators: John R. Brophy and Graeme Aston 


Project Report 


jjsiTBnnucTiQli 

The lifetime of a typical commercial, communica¬ 
tions satellite is limited by the mass of station keeping 
propellant which is carried on-board. For a typical 
communications satellite at geosynchronous altitude 
with a ten year useful life, it is generally true that each 
one second increase in specific impulse results m a 
station keeping propellant mass reduction of approx¬ 
imately one kilogram. A variety of electric propul¬ 
sion thrusters are being developed which promise spe¬ 
cific impulse levels hundreds of seconds greater than 
those possible with chemical thrusters. For a fixed 
launch vehicle delivery mass to geosynchronous orbit, 
the propellant mass savings enabled by high specific 
impulse thrusters can be used to either significantly 
increase the satellite station keeping lifetime or to in¬ 
crease the number of on-board transponders. As a 
result, the economic benefits to commercial users of 
higher specific impulse electric propulsion thrusters 
can be very large. 

The application of electric propulsion to com¬ 
munications satellites, however, has been limited to 
the use of hydrazine thrusters with electric heaters 
for thrust and specific impulse augmentation. These 
electrothermal thrusters operate at specific impulse 
levels of approximately 300 s with heater powers of 
about 500 W. Low power arejets (1-3 kW) are cur¬ 
rently being investigated as a way to increase specific 
impulse levels to approximately 500 s. Ion propulsion 
svsterns can easily produce specific impulses of 3000 s 
or greater, but have yet to be applied to communica¬ 
tions satellites. The reasons most often given for not 
using, ion propulsion systems are their high leve o 
overall complexity, low thrust with long bum times, 
and the difficulty of integrating the propulsion system 
into existing commercial spacecraft busses. 

The Electrostatic Plasma Accelerator (EPA) is a 
thruster concept which promises specific impulse lev¬ 
els between low power arejets and those of the ion 
engine while retaining the relative simplicity of the 
arejet. The EPA thruster produces thrust through 
the electrostatic acceleration of a moderately dense 


plasma. No accelerating electrodes are used and the 
specific impulse is a direct function of the applied dis¬ 
charge voltage and the propellant atomic mass. 


PHENOMENOLOGIC AL MODEL 


The EPA thruster consists of a hollow cathode 
surrounded by a cylindrical discharge chamber shell 
which is completely open at the end opposite the cath¬ 
ode. The discharge chamber walls are maintained at 
anode potential, which may be between 28 to 60 volts 
positive of the cathode. A magnetic field is applied 
to the discharge chamber through the use of external 
permanent magnets to improve the plasma genera¬ 
tion efficiency as well as the efficiency with which the 
plasma is exhausted from the thruster. 


During thruster operation a plasma is created 
y electron bombardment of the propellant gas in a 
lasma arc discharge which is sustained by the ap- 
,lied voltage between the anode and the cathode. The 
,lasma created in the discharge chamber will charac- 
eristically assume a potential slightly higher (a few 
rolts) than the applied voltage. Electrons, by virtue 
,f their much smaller mass, are far more mobile than 
he discharge chamber ions and will attempt to leave 
he discharge chamber plasma. An initial start-up 
ransient will occur in which more electrons will leave 
,he thruster than ions. This will cause the discharge 
chamber plasma to acquire a potential more positive 
;han the ambient space plasma potential to prevent 
the further unbalanced loss of electrons. As the dis- 
:harge plasma potential becomes more positive, ions 
will begin to be accelerated from the thruster. le 
potential difference between the discharge chamber 
plasma and ambient space plasma will automatically 
adjust itself so that the loss rate of ions and electrons 
are equal; thus a neutral plasma is exhausted from 
the thruster. The net ion velocity (and ultimately 
the specific impulse) is determined by the magnitude 


The magnitude of the potential difference be¬ 
tween the discharge chamber plasma and the ambi¬ 
ent space plasma will be approximately equal to the 



applied discharge voltage provided the potential of 
the cathode can be maintained within a few volts 
of the ambient space plasma. This will be the case 
provided sufficient electrons can escape from the dis¬ 
charge chamber to neutralize the electrostatically ac¬ 
celerated ions. In order for an electron to «icape from 
the discharge chamber it must have sufficient energy 
to overcome a potential barrier approximately equal 
to the applied voltage. Only electrons m the high en¬ 
ergy tail of the Maxwellian energy distribution have 
sufficient energy to overcome this potential difference. 
Furthermore, there will be sufficient numbers of such 
electrons only if the electron temperature is above 
some threshold level (which is a function of the po¬ 
tential barrier to be overcome). 

The number of electrons in a Maxwellian energy 
distribution, with energies above a certain level, can 
be readily calculated from statistical thermodynam¬ 
ics. For an applied voltage of 32 volts, the electron 
temperature must be between 4 and 5 eV for there 
to be sufficient electrons to neutralize the ions leav¬ 
ing the thruster. At 60 volts, the required electron 
temperature is between 7 and 8 eV. Electron temper¬ 
atures of these magnitudes can easily be produced m 
an electron bombardment type of discharge chamber. 


TT Thrust Vfecha qfcm Verification 

The thrust producing mechanism for the EPA 
thruster was first discovered by chance in 1979 in 
the course of the extended SERT II (Space Electric 
Rocket Test II) test program. 1 In these tests, one of 
the SERT II mercury ion engines was operated m a 
discharge only mode without the high voltage applied 
to the accelerator system. While operating m this 
mode it was discovered that the thruster was pro 
ducing a measured thrust level of 0.8 mN. Other tests 
indicated that the thrust level decreased as the dis¬ 
charge voltage was reduced (which is consistent w,th 
the phenomological model discussed above). Based 
on the measured thrust of 0.8 mN, Ref. 1 calculated 
a specific impulse of 300 s at an input power of 80 W, 
resulting in a thrust-to power ratio of 10 mN/kW for 
operation in this mode. The 300 s specific impulse 
level (for mercury propellant) is much higher than 
can be accounted for by simple gas d y namlc ex P| n ' 
sion. These tests demonstrate the validity of the EPA 
thruster thrust producing mechanism, even though 
the SERT II ion thruster configuration was far from 
optimum for this mode of operation. 


task OBJECTIVES 

For the EPA thruster to be a viable propul¬ 
sion concept, overall thruster efficiencies approaching 
30approaching 1000 s must be achieved. To achieve 
these efficiency levels, a discharge chamber which can 
produce ions for an energy cost on the order of 120 
eV/ion while ionizing at least 90mechamsms have 
been identified including: the power required to ionize 
the propellant (equivalent to the frozen flow loss in an 
electrothermal or chemical engine), the loss of union¬ 
ized and unaccelerated propellant from the thruster, 
cosine losses resulting from beam divergence, and the 
production of a non-uniform exhaust velocity profile 
due to multiply charged particles. The major effi¬ 
ciency loss mechanisms are believed to be the power 
required to ionize the propellant, and the loss of 
unionized propellant. 

Consideration of the characteristics of the col- 
lisional phenomena which accompany the electron 
bombardment ionization process, together with the 
electron emission processes required to maintain the 
arc discharge, result in a minimum energy expendi¬ 
ture of approximately 60 eV to produce each ion in 
the discharge chamber plasma. It is expected that 
not all of the ions created in the discharge charn 
ber will be accelerated into the exhaust. Some frac¬ 
tion of the ions will recombine with electrons at the 
walls of the discharge chamber. This ion recombina¬ 
tion process will increase the overall average energy 
expenditure per exhaust ion above 60 eV/ion. Ad¬ 
vances in plasma discharge chamber design and un¬ 
derstanding have progressed to the point where is- 
charge chambers for plasma contactors can now ionize 
90133 eV/ion. 2 Thus, the performance objectives for 
the EPA thruster appear to be within reach. 

The specific objectives of the EPA thruster proof 
of concept program are: 

1. Operate the Test Bed Ion (TBI) Engine* as 
an EPA thruster and parametrically vary ie 
thruster geometry and operating conditions to 
quantify the electrostatic plasma acceleration et- 

fect. 

2 Based on the experimental results of (1), develop 
' an analytical model for EPA thruster design and 
performance prediction. 

* This device was developed under a prior 
NASA/LeRC funded program at the Jet Propu sion 
Laboratory and is currently on loan to CAS1 



3. Design and build a radiation cooled EPA thruster 
based upon the output from (1) and (2) and per- 
formance optimize this thruster. 

PRESENT PRO GRAM STATUS 

A brief description of the present program status 
is given below. The Test Bed Ion Engine discharge 
chamber has been received from JPL. A support stand 
has been fabricated and the engine has been mounted 
to the support stand and installed in an oil diffusion 
pumped vacuum facility. A novel hollow cathode con¬ 
figuration was designed for use with this discharge 
chamber and has been fabricated. This cathode uses 
a shaped porous tungsten insert which functions both 
to provide a reservoir of low work function material 
and also as a choked flow orifice to separate the rela¬ 
tively high pressure hollow cathode plasma from the 
much more dilute discharge chamber plasma. 

Eleven ferroresonant DC power supplies have 
been installed in an instrument rack to serve as the 
discharge power supply for the EPA thruster. These 
power supplies can be combined to give a maximum 
current of 275 A at 28 V. A variety of plasma probes 
are being assembled for probing the EPA thruster 
exhaust, including: a Faraday/Langmuir probe rake 
combination and a retarding potential analyzer. The 
Faraday/Langmuir probe rake will be used axial and 
radial position downstream of the thruster. The re¬ 
tarding potential analyzer will be used to verify the 
beam ion energy and any possible distribution in ion 
energies. 


REFERENCES 

1. Kerslake, W. R., and Ignaczak, L. R., ”SERT 
II 1980 Extended Flight Thruster Experiments,” 
AIAA 81-0665, April 1981. 

2. Beattie, J. R-, Williamson, W. S., and Matossian, 
J. N., ”Plasma Contactor Neutralizer System for 
ATAS Missions,” IEPC-88-030, October 1989. 


SUMMARY 

The goal of the present program is to demon¬ 
strate feasibility of the EPA thruster concept. The 
program includes both experimental and theoretical 
investigations of the EPA acceleration mechanism and 
discharge chamber performance. The theoretical in¬ 
vestigations will include the development of a dis¬ 
charge chamber model which describes the relation- 
ships;between the engine size, plasma properties and 
overall performance. One use of the model will be 
to determine the discharge chamber operating condi¬ 
tions which produce sufficiently high electron temper¬ 
atures for ion beam neutralization. The experimental 
investigations will be used to demonstrate ion beam 
acceleration and neutralization, and discharge cham¬ 
ber performance. 





Ion Propulsion 
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INTRODUCTION 

An ion engine is a plasma thruster which pro¬ 
duces thrust by extracting ions from the plasma and 
accelerating them to high velocity with an electro¬ 
static field. The ions are then neutralized and leave 
the engine as high velocity neutral particles. The 
advantages of ion engines are high specific impulse 
and efficiency and their ability to operate over a wide 
range of input powers. In comparison with other elec¬ 
tric thrusters, the ion engine has higher efficiency and 
specific impulse than thermal electric devices such 
as the arcjet, microwave, radiofrequency and laser 
heated thrusters and can operate at much lower cur¬ 
rent levels than the MPD thruster. However, the 
thrust level for an ion engine may be lower than a 
thermal electric thruster of the same operating power, 
consistent with its higher specific impulse, and there¬ 
fore ion engines are best suited for missions which can 
tolerate longer duration propulsive phases. 

The critical issue for the ion engine is lifetime, 
since the prospective missions may require operation 
for several thousands of hours. The critical compo¬ 
nents of the ion engine, with respect to engine life¬ 
time, are the screen and accelerating grid structures. 
Typically, these are large metal screens that must sup¬ 
port a large voltage difference and maintain a small 
gap between them. Metallic whisker growth, distor¬ 
tion and vibration can lead to arcing, and over a long 
period of time ion sputtering will erode the grid struc¬ 
tures and change their geometry. In order to study 
the effects of long time operation of the grid struc- 
ture, we are developing computer codes based on the 
particle-in-cell (PIC) technique and laser induced flu¬ 
orescence (LIF) diagnostic techniques to study the 
physical processes which control the performance and 
lifetime of the grid structures. 


EESEARCH OBJECTIVES 

In order to verify the results of the code and to 


aid in the development of physical models used in 
the code, we are developing an experimental capabil¬ 
ity for the operation and diagnostics of ion engines. 
The test chamber is approximately 3 m in diameter 
by 9 m in length and will be pumped cryogemcally 
using approximately 1 kW of 20 K gaseous helium. 

A three axis computer controlled translation system 
will permit movement of the engine to facilitate pre¬ 
cise spatial diagnostics of the exhaust flow. Mr. Jim 
Sovey of NASA Lewis Research Center has agreed to 
provide us with an ion engine for diagnostic testing. 

In additional to conventional electric probe diag¬ 
nostics, we are developing a LIF technique to measure 
the velocity, density and temperature of the neutral 
and, possibly, ion species in the exhaust flow. The di¬ 
agnostic technique uses the single frequency line from 
a ring-dye laser to excite fluorescence from the species 
of interest. As the laser frequency is scanned across 
the doppler broadened absorption line, the emitted 
radiation is measured using an optical detector. The 
strength of this optical signal is proportional to the 
velocity distribution function of the absorbing species, 
resolved along the laser beam direction. Calculation 
of the zeroth, first and second moments of this mea¬ 
sured distribution provides a measurement of the den¬ 
sity, mean velocity and temperature of the absorbing 

species. 

The basic function of the grid structure is to ex¬ 
tract ions from the plasma maintained in the interior 
of the engine and to accelerate them through the grids 
and out of the engine. In this extraction and accelera¬ 
tion process some of the ions will strike the surfaces of 
the grids and remove metal atoms from the structure 
by sputtering. In addition, some collisions between 
fast ions and neutral particles will occur which result 
in charge exchange that produces slow ions near the 
grid surfaces. These slow ions will be accelerated to 
the grid surfaces by the strong electric field and will 
cause additional sputtering. Collisions between ions 
or electrons and neutral particles, particularly in the 
rare gas propellants, can also produce metastable 



neutral particles which may also produce sputtering 
when they impinge on the grid surfaces. We are de¬ 
veloping a detailed computer code in order to study 
the effect of grid geometry and the imposed potential 
field on these life-limiting processes and the efficiency 
of the acceleration process. 


Particle-in-cell code 

To model ion-thruster performance, we are devel¬ 
oping a particle-in-cell code [1], in which a single pair 
of screen-grid and accelerator-grid holes spans an ax¬ 
ially symmetric computational domain, as illustrated 
in figure 1. Thus, the (generally nonuniform) com¬ 
putational grid extends from the plasma inside the 
ion engine out toward the vacuum of space. Under 
typical operating conditions, the electrical potential 
depends not only on the fixed voltages on the metal 
grids, but also on the extracted ion flux. The particle- 
in-cell (PIC) method is particularly suited to model 
such complex dynamical behavior. Its main features 
are illustrated in figure 2: 

During each iteration in time, one first solves 
Poisson’s equation, whose inhomogeneous part is 
given by the instantaneous charge density. Then, one 
obtains the electric field for each grid point. The 
next two steps are at the core of the PIC-concept: 
Each particle, representing a “cloud” of charge den¬ 
sity, and generally located within some cell of the 
computational grid, is accelerated and displaced ac¬ 
cording to the interpolated values of the electric field 
at the grid points. Also, new particles are added into 
the computational domain by random sampling from 
the plasma phase-space, and particles that impinge 
on the electrodes or exit through any of the other 
boundaries are removed from the computation, io 
close the cycle, one calculates the new charge density 
on the grid points from the local particle positions. 
Interpolations from grid coordinates to particle co¬ 
ordinates and vice versa are performed using linear 
volume-weighting, as illustrated to the right of the 
flow diagram in figure 2. When the system reaches 
steady-state, or at any other point during the compu¬ 
tation'' such quantities as thrust, electrical efficiency, 
and electrode impingement currents can be obtained. 

Many sophistications can be made to this model. 
Already, electrons are taken into account, assuming 
instantaneous thermalization with the ion flow. In 
the near future, neutral atoms will be added to the 
code, whose charge exchange interactions with the 
ions and electrons will be modeled using Monte-Carlo 
techniques. Also, metastable neutrals may be taken 


into account, since these are particularly accessible 
by laser diagnostic techniques. Ultimately, we expect 
the code to predict accurate values of the specific im¬ 
pulse of the device, and reasonable estimates of grid 
lifetimes, so that optimization studies of grid geome¬ 
try and different grid structures (for example, floating 
metallic or dielectric screen grids) can be performed. 

PRFT TMTNARY RESULTS 

To illustrate some of the capabilities of the PIC 
code, a representative calculation for a single-hole ion 
engine operated with neon gas is shown in figure 3. 
This calculation corresponds roughly to the geome¬ 
try shown in figure 1, with a screen grid aperture of 
2 mm, and voltages as indicated in figure 1. Initial 
conditions are for a plasma temperature 2500 K, ion 
density 10 16 m -3 , and a perveance of 10 A/V -. 
Some 15,000 particles were included in the calcula¬ 
tion. In figure 3(a), contours of the electrical poten¬ 
tial are shown. Note that, even downstream from the 
accelerator grid, these contours are bent around the 
centerline, due to the ion density. Contours of the ion 
density are shown in figure 3(b). The impingement on 
the electrodes, on the order of twenty percent in fig¬ 
ure 3, and the divergence of the ion plume are highly 
sensitive to the geometry and operating conditions of 
the ion-thruster. Finally, figure 3(c) shows a plot of 
average velocities of the ions, as obtained from the 
PIC code. 

Not only will the PIC code be run for different 
geometries and operating conditions of the thruster, 
but also comparisons with experimental results are 
planned, notably with work performed at Colorado 
State University by Paul Wilbur and coworkers (see, 
for example, ref. (2)), and, once our own facilities 
become operational, with detailed measurements per- 
formed here at the Space Institute. 

[1] C. K. Birdsall and A. B. Langdon, Plasma 
Physics via Computer Simulation (McGraw-Hill, 
New York, 1985). 

[2] G. Aston, Ion Extraction from a Plasma, Ph D. 
thesis (Colorado State University, Fort Collins, 
1980). 
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Project Title: Magnetic Annual Arc (MAARC) 


Principal Investigator: Herbert Thomas 


Project Report 


Electromagnetic propulsion systems offer high 
specific impulses (> 1500) and high thrust-density as 
compared to conventional chemical rocket propulsion 
systems. Understanding and simulating the physi¬ 
cal processes in electromagnetic propulsion devices 
(as well as testing with refined diagnostic equipment) 
hold the key to improved performance and longevity 
of these devices. This in turn is an essential require¬ 
ment for a practical application in future space en- 
deavors. 

The physical processes involved in electromag¬ 
netic accelerators, especially in Arc-Jets and MPD- 
thrusters, are quite complicated and deny an accurate 
description by analytical means. Therefore it is neces¬ 
sary to develop computer-codes to simulate the main 
effects occurring in electromagnetic accelerators. De¬ 
pending on the density levels at which a particular 
device operates, tow different approaches are neces¬ 
sary. For Low density devices a statistical descrip¬ 
tion of the flow is necessary, while for higher density 
where collisional processes dominate the continuum 
equations are required. The Magnetic Annular Arc 
(MAARC) model being developed assumes a collision 
dominated plasma and the conservation equations for 
mass, momentum and energy including terms for the 
electromagnetic effects are as follows. 


dp . dpu ^ dpv _ Q 
at **" dz * dy 

^ t fjp*' = (/* BU * + £/’- 

d Jgl * ±(p»„) t- + p) = </* 3)„ - §- = r z , + ±r„ 

a y * l {puw} * L (pvw ) = (/* B), t -/» - 

dt “ dz o') 

a l f! - £(.(* -,). - = ikVT) +9 

at ' dz' ' B-J 


Maxwell’s equations in combination with Ohm s 
law are used to determine the electric and magnetic 
fields and the current distribution: 


at 

c~ dt 


While Maxwell’s equations can be used to solve 
electromagnetic problems having a relatively low elec¬ 
trical conductivity, numerical instability problems 
arise if the electrical conductivity is high. The numer¬ 
ical instability prohibits the use of Maxwell’s equa¬ 
tions in their original form for electric thrusters. How¬ 
ever, it is possible to combine Maxwell’s equations 
and Ohm’s law to obtain the general induction equa- 

tion. 

£! +7x £il-7x(ix5)rV* -o 

The resulting induction equation is an elliptic, 
2nd order non-linear differential equation. It is solved 
bv a successive over-relaxation method (SOR), which 
enable the simultaneous solution of a system of tri¬ 
diagonal equations, sweeping along lines of constant 
F and T) (the coordinates in computational domain). 
The effect of an applied magnetic field are calcu¬ 
lated once and are part of the right-hand-side. The 
MAARC geometry is depicted in Figure 1. 

The boundary condition for the insulator results 
from the parallel-current condition, fixing the value of 
the magnetic field times the radius along the insula¬ 
tor. The conditions imposed on the electrodes consist 

of: 

a) the total current passing through the electrodes 
(input variable) 

b) The magnetic field value is zero for downstream 
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c) The electric field is perpendicular to the elec- 
trodes. 

Calculations of a test-case showed convergent be¬ 
havior. Because of the differencing used in those cal¬ 
culation, the accuracy of the magnetic field strengt 
is only approximately one per mille. This inaccuracy 
leads to errors of the order of 1 to 10 for the current 
density and the electric field, which are derived from 
the magnetic field strength. It is therefore necessary 
to reformulate the Induction equation in a strong con¬ 
servative form to reduce the inaccuracy of the mag¬ 
netic field strength, such that the error of the current 
density and the electric field becomes tolerable. This 
is especially important, since the boundary conditions 
are derived from the current densities and the electric 
field values. 

Reformulation of the Induction equation into a 
strong conservative form is now in process. 



Figure 1. MAARC Geometry 
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LASER MATERIALS PROCESSING 


• Effect on Pulse Duty Cycle on Inconel 718 Laser Welds 
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Project Title: Effect of Pulse Duty Cycle on Inconel 71S Laser Welds 

Principal Investigators: M. H. McCay, T. D. McCay and N. B. Dahotre 
Other Investigators: C. M. Sharp, A. Sedghinasab and S. Gopinathan 

Project Report 


ABSTRACT 

Crack sensitive Inconel 718 was laser pulse 
welded using a 3.0 kW C0 2 laser. Weld shape, struc¬ 
ture and porosity were recorded as a function of the 
pulse duty cycle. Within the matrix studied, the 
welds were found to be optimized at a high (17 ms 
on 7 ms off) duty cycle. These welds were superior 
in appearance and lack of porosity to both low duty 
cycle and CW welds. 


BACKGROUND 

Since its development thirty years ago (1), In¬ 
conel 718 has grown in use until it accounts for 35% 
of all superalloy production (2). The material is rou¬ 
tinely employed in gas turbine disks because of its 
resistance to creep and stress rupture at elevated tem¬ 
peratures. Additional uses which take advantage of 
its environmental cracking resistance have developed 
in nuclear power plant and cryogenic facility piping. 

Inconel 718 is known as a nickel-iron base alloy. 
It contains chromium for oxidation resistance, molyb¬ 
denum and tungsten as solid-solution strengthened, 
titanium, aluminum and columbium to combine wit 
nickel as precipitation strengthened, and carbon an 
other elements for various complex effects (3). In con¬ 
trast to alloys that are strengthened by the ordered 
face centered cubic (fee) gamma prime phase, Inconel 
718 is strengthened by the ordered body centered 
tetragonal (bet) gamma double prime. The nominal 
composition is given in Table I. 

When Inconel 718 is welded, metallurgical eval¬ 
uation often reveals microfissures in the heat affected 
zone. In general, these are thought to be due to liqua¬ 
tion, a phenomena that is influenced by m.crostruc- 
ture and constituent inhomogeneities (4). The mi- 
crofissures create problems both in hardware fabri¬ 
cation and repair and, therefore, the welding of In¬ 
conel 718 has continued to be a challenge for both 


researchers and manufacturers. 

Laser welding with its inherently small heat af- 
fected zone has been considered for Inconel 718. A 
study (5) on the welding and cutting of Inconel 718 
with a 15kW CW (continuous wave) C0 2 laser in¬ 
dicated that it could be an effective joining method 
for up to one half inch thicknesses. The CW laser 
welds produced isolated instances of microfissuring 
and some porosity. Several additional studies have 
been made of pulse laser welding on Inconel 718 us¬ 
ing YAG lasers (6,7). In these instances the material 
cracked in the nailhead in both individual and over¬ 
lapping welds. The cracks were intergranular and at¬ 
tributed to solidification cracking. 

The purpose of this investigation is to study the 
solidification morphology and weld structure includ¬ 
ing porosity, of Inconel 718 as the laser welding pro¬ 
cess proceeds from CW to individual non-overlapping 
pubes. 

flYPFRTMENTA Ir PROCEDURE 

WELDING - Bead-on-plate laser welding of the 
Inconel 718 was accomplished with a Rofin Sinar 
RS3000 laser. The RS3000 is a commercially avail¬ 
able fast axial flow, RF (27.12 MHz) excited, C0 2 
(10.6 /un wavelength). For welding the laser may be 
operated in two different configurations: 1) Quasi C W 
- output is pulse width modulated at 5 KHz to give 
control of average power from 0.2 to 3.7 kW; 2) Gated 
Pulse - output can be gated from a CW condition up 
to 25 kHz; peak power remaining fixed at 3. ( kW . L he 
cavity can be configured to produce beams with sev¬ 
eral different spatial modes including TEMoo. TEM io> 
and TEM20- 

One-eighth inch (1/8”) thick wrought alloy 
(ASTM grain size #6) was welded using the standar 
RS3000 welding head (150 mm f/6 optics, circular po¬ 
larization). The welding head nozzle permits slue! - 
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', n .r o-ases, both coaxial flow and plasma control, to 
be”used to moderate the welding process and protect 
the molten metal during solidification. Helium, at a 
flow rate of approximately 60 liters per minute, was 
selected as the shielding gas. The laser beam was di¬ 
rected perpendicular to the workpiece. Figure 1 is a 
schematic of the welding nozzle and the sample orien¬ 
tation. The beam was operated in the TEMio mode 
(ring with a central peak). Gated pulsing was em¬ 
ployed using the parameters shown in Table II. T 
measured peak power, i.e. delivered at the work piece, 
was 3 2 kW. This is below the 3.7 kW output due to 
the optical system losses in the beam delivery system. 
The energy per pulse was kept constant by ^mtain- 
ing the same pulse on-time and only varying the pulse 
off-time. A pulse on-time of 17 ms was necessary £ 
provide complete penetration of the sample. Pulse 
off- time was varied from 1 ms to 21 ms mincrements 
of 2 ms. The samples are designated STA , 

3 STA 21, indicating the pulse off-time used for 
the particular weld. Laser power output calibration 
was done prior to and after welding with a Scientech 
Model 390 water cooled calorimeter. Figure 
the delivered average power as a function o t e se¬ 
lected pulse off-times. For comparison, a CW bead 
on plate weld was made under similar gas ow an 
focus conditions, using 2.5 kW at a weld speed of 
1 5 m/min. Several individual (isolated) pulses were 
also made for characterization, a function of pulse 
off-time. 


DIAGNOSTICS - Elasma Spectroscopy - Th « ra - 
diation emitted by the metallic plasma formed above 
the weld pool was focussed onto the entrance slit of 
a 1/4 meter spectrometer. The transfer optics con¬ 
sisted of a 6:1 telescope which reduced the size o 
the plasma image to about 1 mm at the plane o 
entrance slit. A grating with 1200 lines/mm was chc> 
sen for obtaining high resolution spectra, providing 
spectral window of approximately 60 nrm Five sepa¬ 
rate spectral windows of constant extent ( nm) we 
used to examine the many spectral lines of interest, 
most of which are in the 250 - 500 nm range. An av¬ 
erage of 100 temporal (20 msec each) scans over the 
length'of the bead was obtained by using an optical 
multichannel analyzer (OMA). This enhanced the sig¬ 
nal to noise ratio and negated the transients within 
the plasma which result from flow field, me ting an 
solidification process effects. Each spectra ™ o 
was calibrated using suitable emission line calibration 
tamps. The spectral data was later transferred to a 
personal computer for post-facto reduction and anal- 

vsis. 


Since the nature of the decaying (during the 
beam off-time) plasma may effect the coupling of the 
next laser pulse to the workpiece, the plasrna rate o 
decay is of significance in pulsed laser welding. It 
is currently impossible to model such a process due 
to unknowns such as the appropriate recombination 
coefficients, collision cross sections and other atomic 
properties. This is chiefly due to the plasma contain¬ 
ing most of the constituents of the substance being 
welded as well as any shield gas. In order to esti¬ 
mate the rate of decay of the plasmas an experiment 
was set up wherein the radiation was monitored us¬ 
ing a nanosecond response silicon p-i-n photodiode. 
The photodiode output was stored on an oscilloscope. 
This data allowed the time constant of the plasma de- 
cay to be estimated. 

Video -A videocamera with a 4.2 neutral density 
filter recorded the plasma image on videotape during 
the weld passes. The video frames were subsequently 
digitized and analyzed using a VICOM image proces- 
sor. 

q a rnnle Preheat - The effect of P^heating/cool- 
ing on the microstructure of the Inconel -18 weldment 
JL also studied. One of the samples was heated to 
730°C in an electric resistance furnace. After soakin 0 
for one hour, it was immediately subjected to laser 
pulse welding using a pulse off-time of 21ms. A sec¬ 
ond sample was immersed in liquid Nitrogen or 
hour and then subjected to laser pulse welding using 
the above conditions. A third specimen of Inconel 
was pulse welded at room temperature, also using the 
same conditions. 

M^allographv - The welds were metallographi- 
cally examined in the XY, YZ, and XZ planes (See 
Figure 1) using an optical microscope an a Vl 
camera connected to a VICOM digital .mage process¬ 
ing system. Top and bottom weld bead width and the 
variation in bead width (maximum width - minimum 
width) wort measured in the YZ plane. Gas pore dl- 

ameter gas pore spacing and gas pore area fraction 
ameier, 6“ " r G area f rac - 

were measured in the XZ plane, o . , ., 

lion is the total centerline pore area divided by 
total centerline weld area. Since the weld shape and 
solidification structure were expected to var. 
dimensions, successive sections -re examined by pol¬ 
ishing in 0.1 mm increments in the YZ an , P 
for samples with pulse off-t.mes of , ms, 11 ms, lo ms, 

and 21 ms. 
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RESULTS 

Figures 3 and 4 show typical top and bottom 
views, (XY planes) respectively, of the pulsed welds. 
The periodic changes in width correlate with the laser 
pulse rate. Since the workpiece travel speed was 
maintained constant for all of the runs, the physical 

distance between pulses increases with increasing off- 
time. This effect can also be seen in the longitudinal 
views along the weld centerline (XZ plane) in Figure 
5. As the pulse off-time increases, the melt zones from 
the individual pulses begin to separate. 

Figures 6 and 7 show the minimum and maxi¬ 
mum top and bottom (respectively) bead widths as 
the pulse off-time is increased. For both the top and 
bottom beads there is an overall general increase in 
variation of the bead width (Figure 8) which reflects 
the periodic correlation of width with laser pulse rate. 
In addition, pores begin to appear more frequently 
within the weld center as the pulse off-times increase. 
Figure 9 demonstrates the increase in poring, showing 
the pore area fraction as a function of pulse off-time. 

Using an average radius for the large pores of 
0.025 cm and a surface energy value (8) of 1770 
dynes/cm, the pressure across the wall of the gas bub¬ 
ble forming the pore can be calculated from the fol- 
lowing equation: 

P = 2 cr/r 

where a is the surface energy and r is the radius of the 
gas bubble. These values yield a pressure of 106 mm 
Hg. The temperature of nickel at that vapor pressure 
(8) is 2370°C. For the other Inconel 718 constituents 
one obtains temperatures of 2140°C for chromium, 
4200°C for molybdenum, and 2370°C for iron. Ex¬ 
cept for molybdenum, the anticipated weld pool tem¬ 
peratures are thus well above those required to va¬ 
porize the Inconel constituents and produce these gas 
bubbles. The dendrite arm spacings adjacent to the 
pores are significantly smaller than those further dis¬ 
tant, indicating that the cooling rates are greater in 
the vifeinity of the pores. This is consistant with va¬ 
porization of the metal to form the pores. 

These same gases would be expected to make 
up the weld keyhole and the plasma above the sur¬ 
face. Figure 10 presents the radiating constituents 
within the plasma which is comprised primarily of 
neutral atoms, singly ionized metal atoms and elec¬ 
trons. Nickel, iron, chromium and niobium are defi¬ 


nitely identified as being present in the plasma. 

Figures 11 and 12 are mid-section views in the 
XY plane for STA 7 and STA 21, respectively. The 
weld pool in STA 21 appears.to have had sufficient 
time to almost completely solidify between pulses. 
The 100X views show the growth structure and the 
appearance and location of the pores. 

Figure 13 is a series of photographs taken from 
the video image of the plasma obtained during weld¬ 
ing. The plasma ignites, grows, encounters the de- 
struct jet, and dissipates with each laser pulse. This 
can also be seen in the photodiode data (Figure 14). 
The central (plateau-shaped) trace in the figure is the 
excitation pulse for the single laser pulse indicated by 
the arrow at the bottom of the scope display. The 
plasma ignites within 2 milliseconds behind the beam 
“on” trigger and dissipates within 2 milliseconds after 
the trigger is “off”. 

Successive cross-section (YZ plane) polishes can 
be used to reveal the change of weld pool geometry 
throughout a laser pulse. Figure 15 shows schematic 
cross-sections throughout one pulse for two different 
duty cycles. Figure 16 is an optical micrograph show¬ 
ing a typical cross-section for CW and a single (1< 
msec) laser pulse. The geometry for a series of pulses 
varied with pulse off-time. The short off-time sample 
had the most uniform depth weld pool geometry with 
significantly less variation. The weld middle width of 
STA 7 was 0.8740 ± 0.0320 mm and that of STA 21 
was 0.7312 ± 0.1295 mm over the duration of a pulse. 

The width to depth ratios for the top nailhead 
are presented in Figure 17. As the pulse off-times 
increase, so do both the ratio and its maximum and 
minimum values. The nailhead top width (W T ) to 
neck width (W iV ) is 1 353 ± 0.077 mm for STA i and 
decreases to 1.73 ± 0.244 for STA 21 where 0.0 ti mm 
and 0 244 mm represents the standard deviation in 
the measured widths. The W T /W jV ratio was largest 
for the CW sample, being 2.269 ± 0.088. This ra¬ 
tio for a single pulse (17 msec) was 1.995. Since 
convection exaggerates the nailhead, its effect must 
be lessening as pulse off-time increases. This is con¬ 
firmed by a similar reduction in the amount of band¬ 
ing/perturbations in the solidifation process. The 
number of visible bands was highest for STA 7 (six) 
and lowest for STA 21 (one). 

The short off-time weld (STA 7) is fairly uni¬ 
form with columnar grains growing toward the weld 
centerline. There is a distinguishable nailhead the 
long off-time weld (STA 21) shows considerable van- 
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ation in size and shape. Columnar solidification pro¬ 
ceeds perpendicular to the laser beam. These grains 
appear equiaxed on cross-sections removed from the 
pulse centerline. The nailhead is small and in some 
sections almost non-existant. 

CONCLUDING REMARKS 

Figure 18 is a schematic of the welding process 
as it occurs for pulses with two different off-times. 
The laser beam ignites a metallic plasma which in¬ 
fluences the thermal coupling to the alloy. Although 
the metallic plasma extinguishes between the STA 7 
pulses, the molten zone is continuous, producing a 
large nailhead with a small width to depth ratio. So¬ 
lidification occurs at a consistent angle with respect 
to the weld direction. There is significant convective 
motion as evidenced by the size of the nailhead, the 
solidification banding, and the minimal porosity. 

As the pulse off-time increases (e.g. STA 21), the 
size, shape and continuity of the molten zone begins 
to vary. The nailhead area decreases, but the width 
to depth ratio increases. The plasma extinguishes 
between pulses and the entire molten region nearly 
solidifies. As a result, the solidification direction be¬ 
comes perpendicular to the beam and convection is 
reduced. Pores caused by vaporization of the alloy 
constituents within the keyhole are trapped by the 
more rapid solidification and the lessened convection. 

The pulsed welds at short off-times were signifi¬ 
cantly superior to both the long off-time pulsed welds 
and CW welds. A more definitive study to opti¬ 
mize the laser weld parameters is now underway. The 
cracking behavior of the pulsed welds will be reported 
at the completion of that work. 
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Table I. Composition of Inconel 718 


ELEMENT 

WT% 

Ni 

53.0 

Cr 

18.6 

Mo t 

3.10 

Nb 

5.00 

Al 

0.40 

Ti 

0.90 

Fe 

18.5 

Mn 

0.20 

Si 

0.30 

C 

0.04 


CG-U98 


boam 




Fig. 1 - Schematic of Welding Nozzle and Sample 


Table II. Inconel 718 Pulse Weld Parameters 


Variable 

UNITS 

Weld Speed 

1.5 m/m in 

Focus 

1 mm Below XY Plane 

Pulse On Time 

17 ms 

Pulse Off Time 

1-21 ms (by 2 ms) 

Peak Power 

3.2 kW 

Gas 

Helium 

Flow Rate 

60 £/min 


CG-UTo 
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Fig 7 * Maximum and minimum bottom bead 
widths as a function of puke off-time. 


Pulse Off-Time (ms) 

Fig. 8 - Variation in top and bottom bead widths 
as a function of pulseoff-time. 



PULSE OFF-TIME (ms) 

Fig. 9 - Gas pore area fraction as a function of 
laser pulse off-time. 



RELATIVE INTENSITY CT RELATIVE INTENSITY 0) RELATIVE INTENSITY 



WAVELENGTH (A) 



Fig. 10 - The emission spectrum of Inconel 718 
in helium shield gas for (a) the 
spectral range of 2500 to 5000 A°, 

(b) the spectral range of 3500 to 
4000 A°, and (c) the spectral 
range of 4000 to 4500 A°. 
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13 - Sequential photographs of plasma 
shapeduring a single pulse for (a) 17 
ms on / 7 ms off and (b) 17 ms on / 
28 ms off. 
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Oscilloscope trace of photodiode 
output during (a) single pulse and 
(b) multiple pulses (21 ms 
off-time). 


Fig. 15 - Typical cross section grain structures 
throughout one pulse for (a) b iA 

7 and (b) STA 21. 
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Fig- 17 - The top nailhead width to depth 

ratio as a function of pulse off-time. 
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Fig. 18 - Schematic of pulsed welding process 
with a representative temperature 
profile for (a) short (7 ms) and 
(b) long (21 ms) duration pulse 
off-times. 
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